restart: lstd:={ } ;`(* Menge der Differentialgleichungen *)`: tabf:={ } ;`(* Liste der zu berechnenden Funktionen *)`: lsg:=dsolve(lstd,tabf); fx:=subs(lsg,x(t)); fy:=subs(lsg,y(t)); plot([fx,fy,t=0..2*Pi],scaling=constrained);
Beispiel:
Randbedingungen:
restart: lstd:={D(x)(t)=2*y(t),D(y)(t)=-0.5*x(t),x(0)=1,D(y)(1)=2}; tabf:={x(t),y(t)}; lsg:=dsolve(lstd,tabf); fx:=subs(lsg,x(t)); fy:=subs(lsg,y(t)); plot([fx,fy,t=0..2*Pi],scaling=constrained);