Plotting ODEs on maple

Question

I'm a beginner on maple and i don't know how to plot an ODEs system like this.

Here's the code

g := 9.81;
m := 0.25e-2;
h := 0.73e-3;
xi := -.5;
yi := -.5;
zi := .3;
vxi := 8;
vyi := 2;
vzi := 3;
wx := (0*2)*Pi;
wy := (13*2)*Pi;
wz := (0*2)*Pi;
eqx := m*(diff(x(t), t, t)) = -h*sqrt((diff(x(t), t))^2+(diff(y(t), t))^2+   (diff(z(t), t))^2)*(diff(x(t), t))+a*(wy*(diff(z(t), t))-wz*(diff(y(t), t)));
eqy := m*(diff(y(t), t, t)) = -h*sqrt((diff(x(t), t))^2+(diff(y(t), t))^2+(diff(z(t), t))^2)*(diff(y(t), t))+a*(wz*(diff(x(t), t))-wx*(diff(z(t), t)));
eqz := m*(diff(z(t), t, t)) = -m*g-h*sqrt((diff(x(t), t))^2+(diff(y(t), t))^2+(diff(z(t), t))^2)*(diff(z(t), t))+a*(wx*(diff(y*t, t))-wy*(diff(x(t), t)));
ei := x(0) = xi, y(0) = yi, z(0) = zi, (D(x))(0) = vxi, (D(y))(0) = vyi, (D(z))(0) = vzi;
F := dsolve({ei, eqx, eqy, eqz}, {x(t), y(t), z(t)}, numeric, output = listprocedure);

Thanks in advance.

Carl Love · Accepted Answer · 2016-06-22 13:13:26Z

0

The command is

plots:-odeplot(F, [x,y,z](t), 0..100);

answered Jun 22, 2016 at 13:13

Carl Love

1,4718 silver badges7 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

gxxx Over a year ago

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

Carl Love Over a year ago

That error message is unrelated to the plotting command. The message is about mathematical properties of the system itself. I'll work on it later tonight.

Collectives™ on Stack Overflow

Plotting ODEs on maple

1 Answer 1

2 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

2 Comments

Your Answer

Sign up or log in

Post as a guest

Related