Метод Пикара

Форум пользователей пакета Maple

Модератор: Admin

Кли
Сообщения: 145
Зарегистрирован: Пт ноя 04, 2016 4:54 pm

Метод Пикара

Сообщение Кли » Вс мар 17, 2019 9:47 pm

Здравствуйте!!! Помогите решить дифф.уравнение методом Пикара, мне не удаётся в цикле вычислить интегралы

Код: Выделить всё

restart;
with(plots):
n:=10:
a:=1; b:=2;
f:=unapply(((3*x^2*exp(-x)-(x+1)*y)/x),x,y);
G:=simplify(dsolve({diff(y(x),x)=f(x,y(x)),y(a)=1}));
N:=5: h:=(b-a)/N;
for i from 0 to N do
x[i]:=a+i*h:
end do;
y[0]:=1:
for k from 1 to N do
y[k]:=evalf(y[0]+int(f(x,y[k-1]),x)-subs(x=a,int(f(x,y[k-1]),x))):
end do;
s1:=plot(rhs(G),x=a..b,color=red):
s2:=plot([seq(y(x)],x=a..b}):
display(s1,s2);