Blame view

Giac_maj/giac-1.4.9/examples/lewisw/lwN1 9.34 KB
6663b6c9   adorian   projet complet av...
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
  res := -l1*l2*p11*p22*q1^2*q4^2+l2*p11*p22*q1^2*q4^2+l1*p11*p22*q1^2*q4^2-p11*p22*q1^2*\
  q4^2+l1*l2*p12*p21*q1^2*q4^2-l2*p12*p21*q1^2*q4^2-l1*p12*p21*q1^2*q4^2+p12*p21*\
  q1^2*q4^2+l1*n22*p11*p22*q1*q4^2-n22*p11*p22*q1*q4^2+l2*n11*p11*p22*q1*q4^2-n11*\
  p11*p22*q1*q4^2-l2*p11*p22*q1*q4^2-l1*p11*p22*q1*q4^2+2*p11*p22*q1*q4^2-l1*n22*\
  p12*p21*q1*q4^2+n22*p12*p21*q1*q4^2-l2*n11*p12*p21*q1*q4^2+n11*p12*p21*q1*q4^2+\
  l2*p12*p21*q1*q4^2+l1*p12*p21*q1*q4^2-2*p12*p21*q1*q4^2+g^2*p11*p22*q4^2-n11*n22*\
  p11*p22*q4^2+n22*p11*p22*q4^2+n11*p11*p22*q4^2-p11*p22*q4^2-g^2*p12*p21*q4^2+\
  n11*n22*p12*p21*q4^2-n22*p12*p21*q4^2-n11*p12*p21*q4^2+p12*p21*q4^2+l1*l2*n11*\
  p22*q1*q2*q3*q4-l2*n11*p22*q1*q2*q3*q4-l1*n11*p22*q1*q2*q3*q4+n11*p22*q1*q2*q3*\
  q4-l1*l2*g*p21*q1*q2*q3*q4+l2*g*p21*q1*q2*q3*q4+l1*g*p21*q1*q2*q3*q4-g*p21*q1*\
  q2*q3*q4-l1*l2*g*p12*q1*q2*q3*q4+l2*g*p12*q1*q2*q3*q4+l1*g*p12*q1*q2*q3*q4-g*p12*\
  q1*q2*q3*q4+l1*l2*n22*p11*q1*q2*q3*q4-l2*n22*p11*q1*q2*q3*q4-l1*n22*p11*q1*q2*q3*\
  q4+n22*p11*q1*q2*q3*q4+l1*g^2*p22*q2*q3*q4-g^2*p22*q2*q3*q4-l1*n11*n22*p22*q2*\
  q3*q4+n11*n22*p22*q2*q3*q4+l1*n11*p22*q2*q3*q4-n11*p22*q2*q3*q4-l2*g*p21*q2*q3*\
  q4+g*p21*q2*q3*q4-l1*g*p12*q2*q3*q4+g*p12*q2*q3*q4+l2*g^2*p11*q2*q3*q4-g^2*p11*\
  q2*q3*q4-l2*n11*n22*p11*q2*q3*q4+n11*n22*p11*q2*q3*q4+l2*n22*p11*q2*q3*q4-n22*\
  p11*q2*q3*q4-l1*l2*n11*p22*q1*q3*q4+l1*n11*p22*q1*q3*q4+l1*l2*p22*q1*q3*q4-l1*\
  p22*q1*q3*q4+l1*l2*g*p21*q1*q3*q4-l1*g*p21*q1*q3*q4+l1*l2*g*p12*q1*q3*q4-l2*g*\
  p12*q1*q3*q4-l1*l2*n22*p11*q1*q3*q4+l2*n22*p11*q1*q3*q4+l1*l2*p11*q1*q3*q4-l2*\
  p11*q1*q3*q4-l1*g^2*p22*q3*q4+l1*n11*n22*p22*q3*q4-l1*n22*p22*q3*q4-l1*n11*p22*\
  q3*q4+l1*p22*q3*q4-l2*g^2*p11*q3*q4+l2*n11*n22*p11*q3*q4-l2*n22*p11*q3*q4-l2*\
  n11*p11*q3*q4+l2*p11*q3*q4-l1*n22*p11*p22*q1*q2*q4+n22*p11*p22*q1*q2*q4-l2*n11*\
  p11*p22*q1*q2*q4+n11*p11*p22*q1*q2*q4+l2*n11*p22*q1*q2*q4-n11*p22*q1*q2*q4+l1*\
  n22*p12*p21*q1*q2*q4-n22*p12*p21*q1*q2*q4+l2*n11*p12*p21*q1*q2*q4-n11*p12*p21*q1*\
  q2*q4-l1*g*p21*q1*q2*q4+g*p21*q1*q2*q4-l2*g*p12*q1*q2*q4+g*p12*q1*q2*q4+l1*n22*\
  p11*q1*q2*q4-n22*p11*q1*q2*q4-2*g^2*p11*p22*q2*q4+2*n11*n22*p11*p22*q2*q4-n22*p11*\
  p22*q2*q4-n11*p11*p22*q2*q4+g^2*p22*q2*q4-n11*n22*p22*q2*q4+n11*p22*q2*q4+2*g^2*\
  p12*p21*q2*q4-2*n11*n22*p12*p21*q2*q4+n22*p12*p21*q2*q4+n11*p12*p21*q2*q4-g*p21*\
  q2*q4-g*p12*q2*q4+g^2*p11*q2*q4-n11*n22*p11*q2*q4+n22*p11*q2*q4+2*l1*l2*p11*p22*\
  q1^2*q4-l2*p11*p22*q1^2*q4-l1*p11*p22*q1^2*q4-l1*l2*p22*q1^2*q4+l1*p22*q1^2*q4-\
  2*l1*l2*p12*p21*q1^2*q4+l2*p12*p21*q1^2*q4+l1*p12*p21*q1^2*q4-l1*l2*p11*q1^2*q4+\
  l2*p11*q1^2*q4-l1*n22*p11*p22*q1*q4-l2*n11*p11*p22*q1*q4+l2*p11*p22*q1*q4+l1*p11*\
  p22*q1*q4+l1*n22*p22*q1*q4-l1*p22*q1*q4+l1*n22*p12*p21*q1*q4+l2*n11*p12*p21*q1*\
  q4-l2*p12*p21*q1*q4-l1*p12*p21*q1*q4+l1*g*p21*q1*q4+l2*g*p12*q1*q4+l2*n11*p11*\
  q1*q4-l2*p11*q1*q4+l1*l2*g^2*q2^2*q3^2-l2*g^2*q2^2*q3^2-l1*g^2*q2^2*q3^2+g^2*\
  q2^2*q3^2-l1*l2*n11*n22*q2^2*q3^2+l2*n11*n22*q2^2*q3^2+l1*n11*n22*q2^2*q3^2-n11*\
  n22*q2^2*q3^2-2*l1*l2*g^2*q2*q3^2+l2*g^2*q2*q3^2+l1*g^2*q2*q3^2+2*l1*l2*n11*n22*\
  q2*q3^2-l2*n11*n22*q2*q3^2-l1*n11*n22*q2*q3^2-l1*l2*n22*q2*q3^2+l1*n22*q2*q3^2-\
  l1*l2*n11*q2*q3^2+l2*n11*q2*q3^2+l1*l2*g^2*q3^2-l1*l2*n11*n22*q3^2+l1*l2*n22*q3^2+\
  l1*l2*n11*q3^2-l1*l2*q3^2-l1*g^2*p22*q2^2*q3+g^2*p22*q2^2*q3+l1*n11*n22*p22*q2^2*\
  q3-n11*n22*p22*q2^2*q3-l2*g^2*p11*q2^2*q3+g^2*p11*q2^2*q3+l2*n11*n22*p11*q2^2*\
  q3-n11*n22*p11*q2^2*q3+l2*g^2*q2^2*q3+l1*g^2*q2^2*q3-2*g^2*q2^2*q3-l2*n11*n22*\
  q2^2*q3-l1*n11*n22*q2^2*q3+2*n11*n22*q2^2*q3-l1*l2*n11*p22*q1*q2*q3+l2*n11*p22*\
  q1*q2*q3+l1*l2*g*p21*q1*q2*q3-l2*g*p21*q1*q2*q3+l1*l2*g*p12*q1*q2*q3-l1*g*p12*\
  q1*q2*q3-l1*l2*n22*p11*q1*q2*q3+l1*n22*p11*q1*q2*q3+l1*l2*n22*q1*q2*q3-l1*n22*\
  q1*q2*q3+l1*l2*n11*q1*q2*q3-l2*n11*q1*q2*q3+l1*g^2*p22*q2*q3-l1*n11*n22*p22*q2*\
  q3+l1*n22*p22*q2*q3+l2*g*p21*q2*q3+l1*g*p12*q2*q3+l2*g^2*p11*q2*q3-l2*n11*n22*\
  p11*q2*q3+l2*n11*p11*q2*q3-l2*g^2*q2*q3-l1*g^2*q2*q3+l2*n11*n22*q2*q3+l1*n11*\
  n22*q2*q3-l1*n22*q2*q3-l2*n11*q2*q3+l1*l2*n11*p22*q1*q3-l1*l2*p22*q1*q3-l1*l2*\
  g*p21*q1*q3-l1*l2*g*p12*q1*q3+l1*l2*n22*p11*q1*q3-l1*l2*p11*q1*q3-l1*l2*n22*q1*\
  q3-l1*l2*n11*q1*q3+2*l1*l2*q1*q3+g^2*p11*p22*q2^2-n11*n22*p11*p22*q2^2-g^2*p22*\
  q2^2+n11*n22*p22*q2^2-g^2*p12*p21*q2^2+n11*n22*p12*p21*q2^2-g^2*p11*q2^2+n11*\
  n22*p11*q2^2+g^2*q2^2-n11*n22*q2^2+l1*n22*p11*p22*q1*q2+l2*n11*p11*p22*q1*q2-\
  l1*n22*p22*q1*q2-l2*n11*p22*q1*q2-l1*n22*p12*p21*q1*q2-l2*n11*p12*p21*q1*q2-l1*\
  n22*p11*q1*q2-l2*n11*p11*q1*q2+l1*n22*q1*q2+l2*n11*q1*q2-l1*l2*p11*p22*q1^2+l1*\
  l2*p22*q1^2+l1*l2*p12*p21*q1^2+l1*l2*p11*q1^2-l1*l2*q1^2;
  
   s1 := (l2*g*a22^3-l1*g*a12*a21*a22^2-l2*n22*a21*a22^2+l2*n11*a21*a22^2+2*l2*g*a12*a22^2-\
  l1*g*a12*a22^2+l2*g^2*a22^2-l2*n11*n22*a22^2+l2*n11*a22^2+l1*n22*a12*a21^2*a22-\
  l1*n11*a12*a21^2*a22-l2*g*a21^2*a22-l1*g*a12^2*a21*a22-l2*g^2*a12*a21*a22-l1*g^2*\
  a12*a21*a22+l2*n11*n22*a12*a21*a22+l1*n11*n22*a12*a21*a22-2*l2*n22*a12*a21*a22+\
  l1*n22*a12*a21*a22+l2*n11*a12*a21*a22-2*l1*n11*a12*a21*a22-l2*g*a21*a22+l2*g*a12^2*\
  a22-l1*g*a12^2*a22+l2*g^2*a12*a22-l1*g^2*a12*a22-l2*n11*n22*a12*a22+l1*n11*n22*\
  a12*a22+l2*n11*a12*a22-l1*n11*a12*a22+l1*g*a12*a21^3+l1*g^2*a12^2*a21^2-l1*n11*\
  n22*a12^2*a21^2+l1*n22*a12^2*a21^2-l2*g*a12*a21^2+2*l1*g*a12*a21^2-l2*g^2*a12^2*\
  a21+l1*g^2*a12^2*a21+l2*n11*n22*a12^2*a21-l1*n11*n22*a12^2*a21-l2*n22*a12^2*a21+\
  l1*n22*a12^2*a21-l2*g*a12*a21+l1*g*a12*a21)/(l2*g*a12*a21*a22^2-l1*g*a12*a21*a22^2-\
  l2*n22*a21*a22^2+l1*n22*a21*a22^2+l1*l2*g*a12*a22^2-l1*g*a12*a22^2-l1*l2*n22*a22^2+\
  l1*n22*a22^2+l2*n11*a12*a21^2*a22-l1*n11*a12*a21^2*a22-l2*g*a21^2*a22+l1*g*a21^2*\
  a22-l1*l2*g*a12^2*a21*a22+2*l2*g*a12^2*a21*a22-l1*g*a12^2*a21*a22+l1*l2*n22*a12*\
  a21*a22-2*l2*n22*a12*a21*a22+l1*n22*a12*a21*a22+l1*l2*n11*a12*a21*a22+l2*n11*a12*\
  a21*a22-2*l1*n11*a12*a21*a22-l1*l2*g*a21*a22-l2*g*a21*a22+2*l1*g*a21*a22+l1*l2*\
  g*a12^2*a22-l1*g*a12^2*a22-l1*l2*n22*a12*a22+l1*n22*a12*a22+l1*l2*n11*a12*a22-\
  l1*n11*a12*a22-l1*l2*g*a22+l1*g*a22-l1*l2*n11*a12^2*a21^2+l2*n11*a12^2*a21^2+\
  l1*l2*g*a12*a21^2-l2*g*a12*a21^2-l1*l2*g*a12^3*a21+l2*g*a12^3*a21+l1*l2*n22*a12^2*\
  a21-l2*n22*a12^2*a21-l1*l2*n11*a12^2*a21+l2*n11*a12^2*a21+l1*l2*g*a12*a21-l2*\
  g*a12*a21);
  
   s2 := (l2*g*a12*a22^2-l2*n22*a22^2+l2*a22^2-l1*g*a12^2*a21*a22+l1*n22*a12*a21*a22+\
  l2*n11*a12*a21*a22-l2*a12*a21*a22-l1*a12*a21*a22-l2*g*a21*a22+l2*g*a12^2*a22-\
  l1*g*a12^2*a22-l2*n22*a12*a22+l1*n22*a12*a22+l2*a12*a22-l1*a12*a22-l1*n11*a12^2*\
  a21^2+l1*a12^2*a21^2+l1*g*a12*a21^2+l2*n11*a12^2*a21-l1*n11*a12^2*a21-l2*a12^2*\
  a21+l1*a12^2*a21-l2*g*a12*a21+l1*g*a12*a21)/(l2*g*a12*a22^2-g*a12*a22^2-l2*n22*\
  a22^2+n22*a22^2-l1*g*a12^2*a21*a22+g*a12^2*a21*a22+l1*n22*a12*a21*a22-n22*a12*\
  a21*a22+l2*n11*a12*a21*a22-n11*a12*a21*a22-l2*g*a21*a22+g*a21*a22+l2*g*a12^2*\
  a22-g*a12^2*a22-l2*n22*a12*a22+n22*a12*a22+l2*n11*a12*a22-n11*a12*a22-l2*g*a22+\
  g*a22-l1*n11*a12^2*a21^2+n11*a12^2*a21^2+l1*g*a12*a21^2-g*a12*a21^2-l1*g*a12^3*\
  a21+g*a12^3*a21+l1*n22*a12^2*a21-n22*a12^2*a21-l1*n11*a12^2*a21+n11*a12^2*a21+\
  l1*g*a12*a21-g*a12*a21);
  
   s3 := (l2*p21*a22^3-l1*p21*a12*a21*a22^2-l2*p22*a21*a22^2+l2*p11*a21*a22^2+2*l2*p21*\
  a12*a22^2-l1*p21*a12*a22^2-l2*p11*p22*a22^2+l2*p12*p21*a22^2+l2*p11*a22^2+l1*\
  p22*a12*a21^2*a22-l1*p11*a12*a21^2*a22-l2*p12*a21^2*a22-l1*p21*a12^2*a21*a22+l2*\
  p11*p22*a12*a21*a22+l1*p11*p22*a12*a21*a22-2*l2*p22*a12*a21*a22+l1*p22*a12*a21*a22-\
  l2*p12*p21*a12*a21*a22-l1*p12*p21*a12*a21*a22+l2*p11*a12*a21*a22-2*l1*p11*a12*a21*\
  a22-l2*p12*a21*a22+l2*p21*a12^2*a22-l1*p21*a12^2*a22-l2*p11*p22*a12*a22+l1*p11*\
  p22*a12*a22+l2*p12*p21*a12*a22-l1*p12*p21*a12*a22+l2*p11*a12*a22-l1*p11*a12*a22+\
  l1*p12*a12*a21^3-l1*p11*p22*a12^2*a21^2+l1*p22*a12^2*a21^2+l1*p12*p21*a12^2*a21^2-\
  l2*p12*a12*a21^2+2*l1*p12*a12*a21^2+l2*p11*p22*a12^2*a21-l1*p11*p22*a12^2*a21-l2*\
  p22*a12^2*a21+l1*p22*a12^2*a21-l2*p12*p21*a12^2*a21+l1*p12*p21*a12^2*a21-l2*p12*\
  a12*a21+l1*p12*a12*a21)/(l2*p21*a12*a21*a22^2-l1*p21*a12*a21*a22^2-l2*p22*a21*\
  a22^2+l1*p22*a21*a22^2+l1*l2*p21*a12*a22^2-l1*p21*a12*a22^2-l1*l2*p22*a22^2+l1*\
  p22*a22^2+l2*p11*a12*a21^2*a22-l1*p11*a12*a21^2*a22-l2*p12*a21^2*a22+l1*p12*a21^2*\
  a22-l1*l2*p21*a12^2*a21*a22+2*l2*p21*a12^2*a21*a22-l1*p21*a12^2*a21*a22+l1*l2*p22*\
  a12*a21*a22-2*l2*p22*a12*a21*a22+l1*p22*a12*a21*a22+l1*l2*p11*a12*a21*a22+l2*p11*\
  a12*a21*a22-2*l1*p11*a12*a21*a22-l1*l2*p12*a21*a22-l2*p12*a21*a22+2*l1*p12*a21*a22+\
  l1*l2*p21*a12^2*a22-l1*p21*a12^2*a22-l1*l2*p22*a12*a22+l1*p22*a12*a22+l1*l2*p11*\
  a12*a22-l1*p11*a12*a22-l1*l2*p12*a22+l1*p12*a22-l1*l2*p11*a12^2*a21^2+l2*p11*\
  a12^2*a21^2+l1*l2*p12*a12*a21^2-l2*p12*a12*a21^2-l1*l2*p21*a12^3*a21+l2*p21*a12^3*\
  a21+l1*l2*p22*a12^2*a21-l2*p22*a12^2*a21-l1*l2*p11*a12^2*a21+l2*p11*a12^2*a21+\
  l1*l2*p12*a12*a21-l2*p12*a12*a21);
  
   s4 := (l2*p21*a12*a22^2-l2*p22*a22^2+l2*a22^2-l1*p21*a12^2*a21*a22+l1*p22*a12*a21*a22+\
  l2*p11*a12*a21*a22-l2*a12*a21*a22-l1*a12*a21*a22-l2*p12*a21*a22+l2*p21*a12^2*\
  a22-l1*p21*a12^2*a22-l2*p22*a12*a22+l1*p22*a12*a22+l2*a12*a22-l1*a12*a22-l1*\
  p11*a12^2*a21^2+l1*a12^2*a21^2+l1*p12*a12*a21^2+l2*p11*a12^2*a21-l1*p11*a12^2*\
  a21-l2*a12^2*a21+l1*a12^2*a21-l2*p12*a12*a21+l1*p12*a12*a21)/(l2*p21*a12*a22^2-\
  p21*a12*a22^2-l2*p22*a22^2+p22*a22^2-l1*p21*a12^2*a21*a22+p21*a12^2*a21*a22+l1*\
  p22*a12*a21*a22-p22*a12*a21*a22+l2*p11*a12*a21*a22-p11*a12*a21*a22-l2*p12*a21*\
  a22+p12*a21*a22+l2*p21*a12^2*a22-p21*a12^2*a22-l2*p22*a12*a22+p22*a12*a22+l2*\
  p11*a12*a22-p11*a12*a22-l2*p12*a22+p12*a22-l1*p11*a12^2*a21^2+p11*a12^2*a21^2+\
  l1*p12*a12*a21^2-p12*a12*a21^2-l1*p21*a12^3*a21+p21*a12^3*a21+l1*p22*a12^2*a21-\
  p22*a12^2*a21-l1*p11*a12^2*a21+p11*a12^2*a21+l1*p12*a12*a21-p12*a12*a21);
  
  time(res:=peval(res, [q1,q2,q3,q4],[s1,s2,s3,s4]));
  res;