function [x,y]=naeuler(dyfun,xspan,y0,h)
%用途: Euler法解常微分方程 y'=f(x,y),y(x0)=y0
%格式: [x,y]=naeuler(dyfun,xspan,y0,h),dyfun為函數f(x,y)，xspan為求解區間[x0,xN]
%      y0為初值y(x0)，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';