function [Xout,Yout]=MyEulerPro(fun,x0,xt,y0,PointNumber)
%MyEulerPro 用改進的歐拉法解微分方程
if nargin<5 | PointNumber<=0 %如果函數僅輸入4個參數值，則PointNumer默認值為100
    PointNumer=100;
end
if nargin<4 %y0默認值為0
    y0=0;
end
h=(xt-x0)/PointNumber;%計算所取的兩離散點之間的距離
x=x0+[0:PointNumber]'*h;%表示出離散的自變量x
y(1,:)=y0(:)';
for i=1:PointNumber    %迭代計算過程
    f1=h*feval(fun,x(i),y(i,:));
    f1=f1(:)';
    f2=h*feval(fun,x(i+1),y(i,:)+f1);
    f2=f2(:)';
    y(i+1,:)=y(i,:)+1/2*(f1+f2);
end
Xout=x;
Yout=y;