% EulerDemo.m
% Illustration des Eulerverfahrens anhand eines
% einfachen Beispieles: d/dt x = 1+x^2+t^3, x(1) = -4
clc
clear
help EulerDemo
xstart=-4
tstart = 1
tend = 2
t = tstart;
x=xstart;
tau = input('Schrittweite eingeben: ')
maxstep = (tend-tstart)/tau;
for istep=1:maxstep
    xi(istep) = x;
    ti(istep) = t;
    f  = 1 + x^2 + t^3;
    x = x + tau * f;
    t = t + tau;
    fprintf('%d %f  %f  \n', istep, t, x)
end

plot(ti,xi)