Метод Рунге-Кутта

@parkito · Регистрация: 22.03.2010

Здравствуйте, помогите пожалуйста решить проблему.
Решаю простой дифур

Matlab M

dy=-2*x*y

Создаю m-файл

Matlab M

1
2
3

function dy=fu(x,y)
dy=-2*y*x^2;
end

Запускаю

Matlab M

1	[T,Y]=ode45(@fu,[0,2],1

Получаю ошибку

Matlab M

1 2	Error in ode45 (line 222) rh = norm(f0 ./ max(abs(y),threshold),inf) / (0.8 * rtol^pow);

Что я делаю неправильно?

@S_el · 26.02.2015, 02:22

Сообщение от parkito

Что я делаю неправильно?

Ошибки я не вижу.попробуйте запустить мой код:

Matlab M

1
2
3

A=dsolve('Dy = -2*x^2*y','y(0)==1','x')
[x,Y]=ode45(@fu,[0,2],1);
plot(x,Y,x,subs(A))

как видите численное решение полностью совпадает с аналитическим.

@parkito · 26.02.2015, 03:09 **[ТС]**

S_el, Такая же ситуация

Зосима · 26.02.2015, 17:10

parkito, странно, у меня работает

@parkito 14 / 13 / 10 Регистрация: 22.03.2010 Сообщений: 695
	26.02.2015, 03:09 [ТС]	3
	S_el, Такая же ситуация 0

Зосима 5142 / 3480 / 356 Регистрация: 02.04.2012 Сообщений: 6,387 Записей в блоге: 16
	26.02.2015, 17:10	4
	parkito, странно, у меня работает 0

Опции темы

Решение