Oui, effectivement, désolé d'avoir était aussi flou.
Code : Tout sélectionner
purge(t_S, t_I1, t_I2, t_I12, t_IPrime1, t_IPrime2, t_I2_1, t_I1_2,
g_1,g_2,g_3,g_4,g_5,g_6,g_7,g_8,g_9,g_10,g_11,g_12,g_13,g_14,g_15,g_16,g_17,g_18,g_19,g_20,g_21,g_22,g_23,g_24,g_25,g_26,g_27,g_28,g_29,g_30,g_31,g_32,g_33,g_34,g_35,g_36,g_37,g_38,g_39,g_40,g_41,g_42,g_43,g_44,g_45,g_46,g_47,g_48,
S, I1, I2, I12, IPrime1, IPrime2, I2_1, I1_2,
q1_S, q1_I1, q1_I2, q1_I12, q1_IPrime1, q1_IPrime2, q1_I2_1, q1_I1_2,
q2_S, q2_I1, q2_I2, q2_I12, q2_IPrime1, q2_IPrime2, q2_I2_1, q2_I1_2,
q12_S, q12_I1, q12_I2, q12_I12, q12_IPrime1, q12_IPrime2, q12_I2_1, q12_I1_2,
a1_S, a1_I1, a1_I2, a1_I12, a1_IPrime1, a1_IPrime2, a1_I2_1, a1_I1_2,
a2_S, a2_I1, a2_I2, a2_I12, a2_IPrime1, a2_IPrime2, a2_I2_1, a2_I1_2,
a12_S, a12_I1, a12_I2, a12_I12, a12_IPrime1, a12_IPrime2, a12_I2_1, a12_I1_2,
t_S_S, t_S_I1, t_S_I2, t_S_I12, t_S_IPrime1, t_S_IPrime2, t_S_I2_1, t_S_I1_2,
t_I1_S, t_I1_I1, t_I1_I2, t_I1_I12, t_I1_IPrime1, t_I1_IPrime2, t_I1_I2_1, t_I1_I1_2,
t_I2_S, t_I2_I1, t_I2_I2, t_I2_I12, t_I2_IPrime1, t_I2_IPrime2, t_I2_I2_1, t_I2_I1_2,
t_I12_S, t_I12_I1, t_I12_I2, t_I12_I12, t_I12_IPrime1, t_I12_IPrime2, t_I12_I2_1, t_I12_I1_2,
t_IPrime1_S, t_IPrime1_I1, t_IPrime1_I2, t_IPrime1_I12, t_IPrime1_IPrime1, t_IPrime1_IPrime2, t_IPrime1_I2_1, t_IPrime1_I1_2,
t_IPrime2_S, t_IPrime2_I1, t_IPrime2_I2, t_IPrime2_I12, t_IPrime2_IPrime1, t_IPrime2_IPrime2, t_IPrime2_I2_1, t_IPrime2_I1_2,
t_I2_1_S, t_I2_1_I1, t_I2_1_I2, t_I2_1_I12, t_I2_1_IPrime1, t_I2_1_IPrime2, t_I2_1_I2_1, t_I2_1_I1_2,
t_I1_2_S, t_I1_2_I1, t_I1_2_I2, t_I1_2_I12, t_I1_2_IPrime1, t_I1_2_IPrime2, t_I1_2_I2_1, t_I1_2_I1_2,
d1_S, d1_I1, d1_I2, d1_I12, d1_IPrime1, d1_IPrime2, d1_I2_1, d1_I1_2,
d2_S, d2_I1, d2_I2, d2_I12, d2_IPrime1, d2_IPrime2, d2_I2_1, d2_I1_2,
d12_S, d12_I1, d12_I2, d12_I12, d12_IPrime1, d12_IPrime2, d12_I2_1, d12_I1_2,
l1, l2, l12,
b1,b2,b12, theta, gamma, G, A, TWH, Dec, L, B
);
F1:= t_S*q1_S + t_I1*q1_I1 + t_I2* q1_I2 + t_I12*q1_I12 + t_IPrime1*q1_IPrime1 + t_IPrime2*q1_IPrime2 + t_I2_1*q1_I2_1 + t_I1_2*q1_I1_2;
F2:= t_S*q2_S + t_I1*q2_I1 + t_I2* q2_I2 + t_I12*q2_I12 + t_IPrime1*q2_IPrime1 + t_IPrime2*q2_IPrime2 + t_I2_1*q2_I2_1 + t_I1_2*q2_I1_2;
F12:= t_S*q12_S + t_I1*q12_I1 + t_I2* q12_I2 + t_I12*q12_I12 + t_IPrime1*q12_IPrime1 + t_IPrime2*q12_IPrime2 + t_I2_1*q12_I2_1 + t_I1_2*q12_I1_2;
Psi_1 := (t_S * g_1 * t_I1 ) * S * B_1 * F1 * a1_S;
Psi_2 := (t_I1 * g_2 * t_I2 ) * I1 * d1_I1;
Psi_3 := (t_I1 * g_3 * t_I2 ) * I1 * B_2 * F2 * a2_I1;
Psi_4 := (t_I1 * g_4 * t_I2 ) * I1 * t_I1_I2;
Psi_5 := (t_I2 * g_5 * t_S ) * I2 * d2_I2;
Psi_6 := (t_S * g_6 * t_I2 ) * S * B_2 * F2 * a2_S;
Psi_7 := (t_I2 * g_7 * t_I1 ) * I2 * d2_I2;
Psi_8 := (t_I2 * g_8 * t_I1 ) * I2 * B_1 * F1 * a1_I2;
Psi_9 := (t_I2 * g_9 * t_I1 ) * I2 * t_I2_I1;
Psi_10 := (t_I1 * g_10 * t_S ) * I1 * d1_I1;
Psi_11 := (t_I1 * g_11 * t_I12 ) * I1 * B_2 * F2 * a2_I1;
Psi_12 := (t_S * g_12 * t_I12 ) * S * B_12 * F12 * a12_S;
Psi_13 := (t_I2 * g_13 * t_I12 ) * I2 * B_1 * F1 * a2_I2;
Psi_14 := (t_I12 * g_14 * t_I1 ) * I12 * d2_I12;
Psi_15 := (t_I12 * g_15 * t_S ) * I12 * d12_I12;
Psi_16 := (t_I12 * g_16 * t_I2 ) * I12 * d1_I12;
Psi_17 := (t_I12 * g_17 * t_IPrime1 ) * I12 * d2_I12;
Psi_18 := (t_IPrime1 * g_18 * t_S ) * IPrime1 * d1_IPrime1;
Psi_19 := (t_I12 * g_19 * t_IPrime2 ) * I12 * d1_I12;
Psi_20 := (t_IPrime2 * g_20 * t_S ) * IPrime2 * d2_IPrime2;
Psi_21 := (t_S * g_21 * t_I1_2 ) * S * B_12 * F12 * a12_S;
Psi_22 := (t_I1 * g_22 * t_I1_2 ) * I1 * B_2 * F2 * a2_I1_2;
Psi_23 := (t_S * g_23 * t_I2_1 ) * S * B_12 * F12 * a12_S;
Psi_24 := (t_I2 * g_24 * t_I2_1 ) * I2 * B_2 * F2 * a1_I2;
Psi_25 := (t_I1_2 * g_25 * t_S ) * I1_2 * d12_I1_2;
Psi_26 := (t_I1_2 * g_26 * t_I1 ) * I1_2 * d2_I1_2;
Psi_27 := (t_I1_2 * g_27 * t_I2 ) * I1_2 * d1_I1_2;
Psi_28 := (t_I2_1 * g_28 * t_S ) * I2_1 * d12_I2_1;
Psi_29 := (t_I2_1 * g_29* t_I2 ) * I2_1 * d1_I2_1;
Psi_30 := (t_I2_1 * g_30 * t_I1 ) * I2_1 * d2_I2_1;;
Psi_31 := (t_I1_2 * g_31 * t_I12 ) * I1_2 * d1_I1_2;
Psi_32 := (t_I1_2 * g_32 * t_I12 ) * I12 * B_2 * F2 * a2_I1_2;
Psi_33 := (t_I1_2 * g_33 * t_I12 ) * I1 * t_I1_2_I12;
Psi_34 := (t_I12 * g_31 * t_I2_1 ) * I12 * d1_I12;
Psi_35 := (t_I12 * g_35 * t_I2_1 ) * I12 * B_2 * F2 * a2_I12;
Psi_36 := (t_I12 * g_36 * t_I2_1 ) * I12 * t_I12_I2_1;
Psi_37 := (t_I1_2 * g_37 * t_I2_2 ) * I1_2 * d1_I1_2;
Psi_38 := (t_I1_2 * g_38 * t_I2_1 ) * I1_2 * B_2 * F2 * a2_I1_2;
Psi_39 := (t_I1_2 * g_39 * t_I2_1 ) * I1_2 * t_I1_2_I2_1;
Psi_40 := (t_I2_1 * g_40 * t_I12 ) * I2_1 * d2_I2_1;
Psi_41 := (t_I2_1 * g_41 * t_I12 ) * I2_1 * B_1 * F1 * a1_I2_1;
Psi_42 := (t_I2_1 * g_42 * t_I12 ) * I2_1 * t_I2_1_I12;
Psi_43 := (t_I12 * g_40 * t_I1_2 ) * I12 * d2_I12;
Psi_44 := (t_I12 * g_44 * t_I1_2 ) * I12 * B_1 * F1 * a1_I12;
Psi_45 := (t_I12 * g_45 * t_I1_2 ) * I12 * t_I12_I1_2;
Psi_46 := (t_I2_1 * g_46 * t_I1_2 ) * I2_1 * d2_I2_1;
Psi_47 := (t_I2_1 * g_47 * t_I1_2 ) * I2_1 * B_1 * F1 * a1_I2_1;
Psi_48 := (t_I2_1 * g_48 * t_I1_2 ) * I2_1 * t_I2_1_I1_2;
afficher("Psi");
Vbase:=[- (Psi_7 + Psi_9 + Psi_14 + Psi_26 + Psi_30 - Psi_2 - Psi_3 - Psi_4 - Psi_10 - Psi_11 - Psi_22 ),
- (Psi_2 + Psi_4 +Psi_16 + Psi_27 + Psi_29 - Psi_5 - Psi_7 - Psi_8 - Psi_9 - Psi_13 - Psi_24),
- (Psi_31 + Psi_33 + Psi_40 + Psi_42 - Psi_14 - Psi_15 - Psi_16 - Psi_17 - Psi_19 - Psi_34 - Psi_35 - Psi_36 - Psi_43 - Psi_44 - Psi_45),
Psi_18,
Psi_20,
-(Psi_34 + Psi_36 + Psi_37 + Psi_39 - Psi_28 - Psi_29 - Psi_30 - Psi_40 - Psi_41 - Psi_42 - Psi_46 - Psi_47 - Psi_48),
-(Psi_43 +Psi_45 + Psi_46 + Psi_48 - Psi_25 - Psi_26 - Psi_27 - Psi_31 - Psi_32 - Psi_33 - Psi_37 - Psi_38 - Psi_39)];
afficher("VBase");
V:=[[diff(Vbase[1],I1),diff(Vbase[1],I2),diff(Vbase[1],I12),diff(Vbase[1],IPrime1),diff(Vbase[1],IPrime2),diff(Vbase[1],I2_1),diff(Vbase[1],I1_2)],
[diff(Vbase[2],I1),diff(Vbase[2],I2),diff(Vbase[2],I12),diff(Vbase[2],IPrime1),diff(Vbase[2],IPrime2),diff(Vbase[2],I2_1),diff(Vbase[2],I1_2)],
[diff(Vbase[3],I1),diff(Vbase[3],I2),diff(Vbase[3],I12),diff(Vbase[3],IPrime1),diff(Vbase[3],IPrime2),diff(Vbase[3],I2_1),diff(Vbase[3],I1_2)],
[diff(Vbase[4],I1),diff(Vbase[4],I2),diff(Vbase[4],I12),diff(Vbase[4],IPrime1),diff(Vbase[4],IPrime2),diff(Vbase[4],I2_1),diff(Vbase[4],I1_2)],
[diff(Vbase[5],I1),diff(Vbase[5],I2),diff(Vbase[5],I12),diff(Vbase[5],IPrime1),diff(Vbase[5],IPrime2),diff(Vbase[5],I2_1),diff(Vbase[5],I1_2)],
[diff(Vbase[6],I1),diff(Vbase[6],I2),diff(Vbase[6],I12),diff(Vbase[6],IPrime1),diff(Vbase[6],IPrime2),diff(Vbase[6],I2_1),diff(Vbase[6],I1_2)],
[diff(Vbase[7],I1),diff(Vbase[7],I2),diff(Vbase[7],I12),diff(Vbase[7],IPrime1),diff(Vbase[7],IPrime2),diff(Vbase[7],I2_1),diff(Vbase[7],I1_2)]];
afficher("V");
invV:= inverse(V);
afficher("invV");
R:=[[r,0,0,0,0,0,0],
[0,r,0,0,0,0,0],
[0,0,r,0,0,0,0],
[0,0,0,r,0,0,0],
[0,0,0,0,r,0,0],
[0,0,0,0,0,r,0],
[0,0,0,0,0,0,r]];
afficher("R");
res := simplifier(developper(solve(det(F*invV-R)=0,r)));
afficher(res);
La fonction "afficher" sert uniquement à suivre l'avancée des calculs.
Le script s'arrête entre afficher("V") et afficher("invV").
J'en déduis donc que le problème vient de la ligne invV:= inverse(V).
J'avoue ne pas comprendre pourquoi après une journée de calculs, je n'ai pas eu de résultat.
Merci de votre aide.