% euler.m 
% MATLAB上提供的“改進的歐拉方法”
function s=euler(fun,x0,xn,y0,n)
% 用euler法計算常微分方程初值問題
% x0和y0分別為初值條件y(x0)=y0中的x0,y0
% xn為x取值區間的最后一個節點的橫坐標值，n為把區間分成的等份數
if nargin<5
    error
    return
end
h=(xn-x0)/n;
for i=1:n
    yp=y0+h*feval(fun,x0,y0);
    x0=x0+h;
    yc=y0+h*feval(fun,x0,yp);
    y0=(yp+yc)/2;
end
s=y0;
return

% 首先建立鏈接向fun的函數
function y=ff2(x,y)
y=y-2*x/y;