2012-06-18 3 views
1

Это было бы его: у '(т) = у (т)/тт^2/г^2 * т у (1) = 1Решение этого дифференциального уравнения с помощью Matlab?

я попробовал: функция hazi3b()

[T, Y] = ode45 (@bfugveny, [1 12], 1); участка (T, Y, 'дх')

конца

и:

функция д = bfugveny (т, у)

д = у (т)/т - т^2/(у^2 * т);

конец

+0

Что пошло не так? –

+0

Попытка доступа к y (1.22); индекс должен быть положительным целым или логическим. Ошибка в bfugveny (строка 3) dy = y (t)/t - t^2/(y^2 * t); Ошибка в ode45 (строка 262) f (:, 2) = feval (odeFcn, t + hA (1), y + f * hB (:, 1), odeArgs {:}); Ошибка в hazi3b (строка 2) [T, Y] = ode45 (@bfugveny, [1 12], 1); –

ответ

2

Вам не нужно писать y(t) в формуле.

y перешел в ваш оракул, уже есть предположение для y -оценено-на-времени-t.

Так что попробуйте

dy = y/t - t^2/(y^2*t); 
+0

Спасибо, это сработало :) –

+0

Так что принимайте ответ .... благодарность хорошая, но прием хорош для вас и человека, который вам помог. Проголосуйте тоже. – duffymo