%歐拉向后公式
function E=Eula2(f,x0,y0,xN,N)
%f一階微分方程的函數(shù)
%x0，y0為初始條件
%xN為取值范圍的一個(gè)端點(diǎn)
%N為區(qū)間個(gè)數(shù)
x=zeros(1,N+1);y=zeros(1,N+1);
x(1)=x0;y(1)=y0;
h=(xN-x0)/N;
for n=1:N
    x(n+1)=x(n)+h;
    %用迭代法求一y（n+1）
   z0=y(n)+h*feval(f,x(n),y(n));
   for k=1:3
       z1=y(n)+h*feval(f,x(n+1),z0);
       if abs(z1-z0)<1e-3
           break;
       end
       z0=z1;
   end
   y(n+1)=z1;
end
E=[x',y'];