MaplePrimes - answers and comments on Question, Spacecurve of Numerical Solution to DE

example

acer — Tue, 02 Oct 2012 06:19:07 Z

restart:

solprocs:=dsolve({diff(x(t),t)=x(t)/diff(z(t),t),
                  diff(x(t),t)=y(t)^3-z(t), diff(y(t),t)=z(t)^4-y(t)^3,
                  x(0)=5,y(0)=4,z(0)=1},numeric,
                  output=listprocedure):

X:=eval(x(t),solprocs):
Y:=eval(y(t),solprocs):
Z:=eval(z(t),solprocs):

plots:-spacecurve([X,Y,Z,0..1],axes=box,labels=[x,y,z]);

acer

Not inventing the wheel

Markiyan Hirnyk — Tue, 02 Oct 2012 08:21:09 Z

Also the example from Maple Help to ?DEplot3d may be of interest.

For acer's example

Markiyan Hirnyk — Tue, 02 Oct 2012 19:33:41 Z

Thanks

jschulzb — Tue, 02 Oct 2012 22:47:10 Z

thanks guys! I ended up using the DE method, but since I wish to animate it later I'm gunna try the spacecurve method now,

animation

acer — Wed, 03 Oct 2012 09:49:07 Z

@jschulzb Note that DEplot and DEplot3d can do animations too.

That may well be an effective and reasonably efficient way to show evolution in time from a/some initial conditions (ICs).

If, on the other hand, you intend on animating with respect to some other parameter that appears in an IVP's DEs or in its ICs and if you need it quicker than DEplot offers then you might take a glance at this.