Mathematik für Physiker - Numerische Physik: Modellierung

512 ANHANG B. MATLAB: THE BASICS
[t,y]=ode23(funk,tspan,ystart);
elseif val == 3
[t,y]=ode113(funk,tspan,ystart);
elseif val == 4
[t,y]=ode15s(funk,tspan,ystart);
elseif val == 5
[t,y]=ode23s(funk,tspan,ystart);
elseif val == 6
[t,y]=ode23t(funk,tspan,ystart);
elseif val == 7
[t,y]=ode23tb(funk,tspan,ystart);
elseif val == 8
[t,y]=eulervorwaerts(funk,xstart,xend,dx,ystart);
elseif val == 9
[t,y]=eulerrueckwaerts(funk,xstart,xend,dx,ystart);
elseif val == 10
[t,y]=leapfrog(funk,xstart,xend,dx,ystart);
elseif val == 11
[t,y]=rungekutta(funk,xstart,xend,dx,ystart);
end
plot(t,y,’-ok’)
xlabel(’x’,’Fontsize’,14);
ylabel(’y’,’Fontsize’,14);
13. März 2007 c○ M.-B. Kallenrode