%Lösning av y'=sin(ty),y(0)=6 0<=t<=2 med Eulers metod
%y(2) tabelleras
tslut=2;
n=5;h=tslut/n;%minska h för noggrannare lösning
antal=6;%antal halveringar av h
tabell=[];
for k=1:antal
    t=0;y=6;
    h=h/2;n=2*n;
    for i=1:n
        y=y+h*sin(t*y);   
        t=t+h;
    end
    fel=y-5.8639;
    tabell=[tabell;[n y fel]];
end
format short e
disp('  antal steg    y med Euler     fel')
tabell
pause
%Richardsonextrapolera 
korr=[0;diff(tabell(:,2))];
yext=tabell(:,2)+korr;
disp('  antal steg    y med Euler     fel       korr         y-ext')
tabell=[tabell korr yext]


