% bakat euler test


f=@(t,y) -8*t*y;

N=80;
h=0.1;

y=1;
t=0;
yv=[];
tv=[];
for i=1:N
    y=(1/(1+8*h*(t+h)))*y;
    t=t+h;
    yv=[yv;y];
    tv=[tv;t];
end

clf;
plot(tv,yv,'*');
hold;
plot(tv,exp(-4*tv.*tv),'r');