function E=N85chapt9ex1_1(a,b,N,ya)
%學(xué)號(hào) 20042374  姓名 鐘健升   作業(yè)號(hào)  85
%p264 第1(1)題
%f是微分方程右端表達(dá)式
%g是微分方程解的表達(dá)式
%a,b是區(qū)間[a,b]中的a,b
%N是等分區(qū)間的個(gè)數(shù)
syms x y
f=x+y; g=-x-1+2*exp(x);
h=(b-a)/N;
y=zeros(1,N+1);
x=zeros(1,N+1);
z=zeros(1,N+1);
m=zeros(1,N+1);
y(1)=ya;
z(1)=ya;
m(1)=0;
x=a:h:b;
for i=1:N
    y(i+1)=y(i)+h*(x(i)+y(i));
    z(i+1)=-x(i+1)-1+2*exp(x(i+1));
    m(i+1)=z(i+1)-y(i+1);
end
E=[x',y',z',m'];
plot(x,y,'r.',x,z,'b*')