%calleuler
clear;clf;
a=0;b=10;N=100;y0=1;h=(b-a)/N;
[X,Y]=euler(a,b,y0,N,'f1');
plot(X,Y,'k-o',X,exact(X),X,-2*h*ones(N+1,1),'k-^',0,-2,'.','markerfacecolor','k')   % 
%plot(X,Y,'k-o',X,exact(X),X,h*(2-2*exact(X)),'k-^',0,-2,'.','markerfacecolor','k')   %  %-2*h*ones(N+1,1),
leg1=[' Euler with h = ',num2str(h)];
leg2=[' Exact '];
legend(leg1,leg2,4)
xlabel('x')
ylabel(' y_{exact}(x) and the Euler Solution Y(x) ')