Maple example

Dette er et eksempel på numerisk løsning af differentialligninger ved Euler og Runge-Kutta metoder ved brug af maple. Eksemplet er opgave 1.30 i A&B s.81. Hvis man kopier alt mellem de to stiplede linier til en fil man f.eks. kalder maple.txt så kan man give kommandoen: xmaple maple.txt og vælge `maple text' i den `pop up box' der kommer. Under EDIT i maple programmet vælger man execute worksheet. Nu kan man så efter behag ændre f(t,x) og alle parametrene.

foreuler= Eulers metode

rk2, rk4= Runge-Kutta til 2. og 4. orden.

------------------------------------

> f:=(t,x)->1-t-(1/2)*x^2;

> Xe:=dsolve({diff(x(t),t)=f(x(t),t), x(0)=0}, x(t), type=numeric,

> method=classical[foreuler], stepsize=0.1);

> Xe(0.2);

> Xrk2:=dsolve({diff(x(t),t)=f(x(t),t), x(0)=0}, x(t), type=numeric,

> method=classical[rk2], stepsize=0.1);

> Xrk2(0.2);

> Xrk4:=dsolve({diff(x(t),t)=f(x(t),t), x(0)=0}, x(t), type=numeric,

> method=classical[rk4], stepsize=0.1);

> Xrk4(0.2);

----------------------------------------


Jan Philip Solovej
Last modified: Fri Mar 6 15:02:15 MET 1998