% Framåt Euler exempel


f=@(t,y) y;


Nv=5:100;
a=0; b=1;

errv=[];
hv=[];
for k=1:length(Nv)
    N=Nv(k);
    h=(b-a)/N;
    hv=[hv;h];
    y=1;
    t=0;

    %yv=[];
    %tv=[];
    for i=1:N
        y=y+h*f(t,y); 
        t=t+h;
        %yv=[yv;y]
        %tv=[tv;t];
    end
    errv=[errv;abs(y-exp(1))];
end
clf;
loglog(hv,errv,'*');
hold on;
loglog(hv,hv,'r--');
ylabel('Felet vid t=1');
xlabel('h = steglängd')

%clf;
%plot(tv,yv,'*');
%hold on;
%plot(tv,exp(tv),'r')

