%改進的euler法解一維微分方程的通用子函數
%函數的輸入變量為：T為x的取值上限，n為把[0，T]區間等份的數目，u0為方程的初始值
function u=feuler(T,n,u0)
h=T/n;
u=[];
u(1)=u0;
for i=1:n-1
    u(i+1)=u(i)+h*f(h*i,u(i));
    u(i+1)=u(i)+(h/2)*(f(h*(i+1),u(i+1))+f(h*i,u(i)));
end