Question: Animating a double pendulum/Numerical Ode

HI everyone, im currently using maple 15 and im trying to animate at least the mass of my pendulum. My problem is that the  coordinate of each mass would be the solution of two numerical ODE and I can't animate this. I am able to plot the solution of my ODE (wich correpond to the two angle in the pendulum), but when it comes to animating... its an another piece of work. So here is what I tried: the code would correspond to the animation two masses... no rope yet.

 restart; with(DEtools, odeadvisor); with(plots); with(plottools);

Systeme := [2*(diff(T1(t), t, t))+cos(T1(t)-T2(t))*(diff(T2(t), t, t))+sin(T1(t)-T2(t))*(diff(T2(t), t))^2+19.6*sin(T1(t)) = 0, diff(T2(t), t, t)+cos(T1(t)-T2(t))*(diff(T1(t), t, t))-sin(T1(t)-T2(t))*(diff(T1(t), t))+9.8*sin(T2(t)) = 0, T1(0) = 1, (D(T1))(0) = 0, T2(0) = 1, (D(T2))(0) = 1];

 sol := dsolve(Systeme, type = numeric);

 plots[odeplot](sol, [T1(t), T2(t)], 0 .. 20); (It works, I see weird graph!)

M1 := proc (x, y) plots[display](plots[pointplot]([x, y], color = red, symbol = solidcircle, symbolsize = 25)) end proc;

M2 := proc (x, y) plots[display](plots[pointplot]([x, y], color = red, symbol = solidcircle, symbolsize = 25)) end proc;

A1 := plots[animate](M1, [cos(T1(t)), sin(T1(t))], t = 0 .. 8);

A2 := plots[animate](M2, [cos(T1(t))+cos(T2(t)), sin(T1(t))+sin(T2(t))], t = 0 .. 8);

plots[display]([A1, A2]);

Error, (in plots/animate) points cannot be converted to floating-point values

Error, (in plots/animate) points cannot be converted to floating-point values

Error, (in plots:-display) expecting plot structures but received: [A1, A2]

Error, (in plots/animate) points cannot be converted to floating-point values

Any ideas? Thanks for your help:)

Please Wait...