function E=euler(f1521,a,b,ya,N)
%  f1521為該問題中的函數f(x,y)。
%  a,b分別為取值范圍的左右端點。
%  ya為給定初值y(a)。
%  N為迭代步數。
%  h為步長。
%  輸出值為對應每個節點的近似值。
h=(b-a)/N;
T=zeros(1,N+1);
Y=zeros(1,N+1);
T=a:h:b;
Y(1)=ya;
for j=1:N
    Y(j+1)=Y(j)+h*feval('f1521',T(j),Y(j));
end
T=[T' Y']