eulerwykres

%wykres

figurę

plot(t, y,'-') xlabel ('czas')

title('Rozwiązanie metodą Eulera')

%błąd

bladl=['Euler zwykła :    ' num2str(100*abs(ye(end)-yl(end))/ye(end))

blad2=['Euler ulepszona:    ' num2str(100*abs(ye(end)-y2(end))/ye(end))

blad3=['Euler wsteczna :    ' num2str(100*abs(ye(end)-y3(end))/ye(end))

blad={bladł;blad2;blad3} ; text(0.05,1.05,blad)