Bonjour à tous !
Dans le cadre de mes TIPE j'aurais besoin de résoudre un système d'équations différentielles non linéaires du second ordre... avec Maple.
Est-ce possible ?
Cela fait plusieurs heures que je cherche comment faire, mais déjà avec deux équations j'obtiens sans cesse le message d'erreur suivant :
"Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations"
Mon but serait de résoudre ce système à trois équations (Z,X, et teta sont les inconnues) :
http://physique-ricochet.societeg.co...uationsMvt.jpg
Voilà ce que j'ai fait avec les deux premières équations seulement (téta étant considérée comme une constante): cela ne fonctionne déjà pas...
(le début est l'initialisation des constantes)
> R := 0.25e-1; h := 0.275e-1; ro := 1; M := .1; C := 1; beta := evalf[10](.2*Pi*1/180); j0 := evalf[10](M*R^2); j1 := evalf[10]((1/4)*M*(R^2+(1/3)*h^2));
>
> sim1 := arccos(1-s/R)-(1-s/R)*sqrt(1-(1-s/R)^2);
> sim := subs(s = abs(z(t))/sin(evalf[10](13*Pi*1/180)), sim1); eqn := [diff(diff(x(t), t), t) = -ro*R*sim*((diff(x(t), t))^2+(diff(z(t), t))^2)*C*sin(beta+evalf[10](13*Pi*1/180))*sin(evalf[10](13*Pi*1/180))/(2*M), diff(diff(z(t), t), t) = ro*R*sim*((diff(x(t), t))^2+(diff(z(t), t))^2)*C*sin(beta+evalf[10](13*Pi*1/180))*cos(evalf[10](13*Pi*1/180))/(2*M)-g*R/R((diff(x(t), t))^2+(diff(z(t), t))^2)^2];
init := x(0) = 0, (D(x))(0) = 2, z(0) = 2.25, (D(z))(0) = .5;
sol := dsolve({eqn, init}, numeric);
sol(1);
Pourriez-vous m'aider, me conseiller ou me donner une syntaxe qui marche svp ?
(Je me suis aidée de l'aide de Maple... sans succès).
Merci beaucoup d'avance.
Remarque : j'utilise Maple 11.
-----