function E=N85chapt9ex5MendEuler(a,b,N,ya)
%學號 20042374  姓名 鐘健升   作業(yè)號  85
%p264 第5題改進Euler方法
%f是微分方程右端表達式
%g是微分方程解的表達式
%a,b是區(qū)間[a,b]中的a,b
%N是等分區(qū)間的個數
syms x
f=exp(x^2);
h=(b-a)/N;
y=zeros(1,N+1);
x=zeros(1,N+1);
y(1)=ya;
x=a:h:b;
for i=1:N
    y1=y(i)+h*exp(x(i)^2);
    y2=y(i)+h*exp(x(i+1)^2);
    y(i+1)=(y1+y2)/2;
end
E=[x',y']