[幻紫风]

对初值问题

y'=1/(1+x*x)-2*y*y,0<x<1;
y(0)=0;
试用Euler法取步长h=0.1和h=0.2计算其近似解，用Matlab实现求解。
源代码:
function[x,y]=Euler1(dyfun,xspan,y0,h)
x=xspan(1):h:xspan(2);
y(1)=y0;
for n=1:length(x)-1
y(n+1)=y(n)+h*feval(dyfun,x(n),y(n));
end
x=x';y=y';