function [x,y,yanal]=problem1(h,xf)
i=1;
t(i)=0;
x(i)=0;
y(i)=1;
while x(i)<xf
y(i+1)=y(i)+h*(x(i)+y(i));
x(i+1)=x(i)+h*y(i);
t(i+1)=t(i)+h;
yanal(i+1)=2*exp(x(i))-x(i)-1;
chi=((yanal(i)-y(i))/(y(i)))*100
i=i+1
end
endfunction
На самом деле я пытаюсь решить дифференциальное уравнение (dy / dt) = t + y, где y (0) = 1; мне нужно найти y (1) для другого размера шага, например, 0,4 0,1 0,05, и если аналитическое решение равно y (t) = 2 * exp (t) -t-1, то найти разницу в% между y и yanaltical.