Метод Эйлера для решения ОДУ

@cfif · Регистрация: 30.11.2015

Author24 — интернет-сервис помощи студентам

дано дифференциальное уравнение вида dy/dx=1+y/x+exp(y/x)
Программу написал, но почему-то не рисует график. Может я ошибся где-нибудь. Помогите пожалуйста

Matlab M

function test_eiler_pr
clc;
f=@(x,y) 1+y/x+exp(y/x);
a=1;b=1.9;
n=20000;
h=(b-a)/n;
y(1)=0;
for k=1:n
    x(k+1)=a+(k-1)*h;
    y(k+1)=y(k)+h*f(x(k),y(k));
end
%lx=length(x);disp('lx');disp(lx);
%ly=length(y);disp('ly');disp(ly);
plot(x,y,'r--');

@VSI · 05.11.2016, 18:11

cfif, Решение ОДУ методом Эйлера в Матлабе

@Centurio · 05.11.2016, 18:28

cfif, у вас не задано x(1), в результате оно принимается равным 0, происходит деление на 0.

@cfif · 06.11.2016, 14:55 **[ТС]**

спасибо, разобрался

@VSI Модератор $Эксперт по математике/физике$ 5240 / 4027 / 1385 Регистрация: 30.07.2012 Сообщений: 12,301
	05.11.2016, 18:11	2
	cfif, Решение ОДУ методом Эйлера в Матлабе 0

@Centurio Модератор 1703 / 1557 / 520 Регистрация: 13.09.2015 Сообщений: 5,379
	05.11.2016, 18:28	3
	cfif, у вас не задано x(1), в результате оно принимается равным 0, происходит деление на 0. 1

@cfif 0 / 0 / 0 Регистрация: 30.11.2015 Сообщений: 33
	06.11.2016, 14:55 [ТС]	4
	спасибо, разобрался 0