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
|
/*
* multiplikation.maxima
*
* (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
*/
Malpha: matrix(
[ 0, 0, 0, 0, 0, 0, -m0 ],
[ 1, 0, 0, 0, 0, 0, -m1 ],
[ 0, 1, 0, 0, 0, 0, -m2 ],
[ 0, 0, 1, 0, 0, 0, -m3 ],
[ 0, 0, 0, 1, 0, 0, -m4 ],
[ 0, 0, 0, 0, 1, 0, -m5 ],
[ 0, 0, 0, 0, 0, 1, -m6 ]
);
Malpha2: expand(Malpha . Malpha);
Malpha3: expand(Malpha . Malpha2);
Malpha4: expand(Malpha . Malpha3);
Malpha5: expand(Malpha . Malpha4);
Malpha6: expand(Malpha . Malpha5);
Malpha7: expand(Malpha . Malpha6);
Malpha8: expand(Malpha . Malpha7);
p: m0 * identfor(Malpha)
+ m1 * Malpha
+ m2 * Malpha2
+ m3 * Malpha3
+ m4 * Malpha4
+ m5 * Malpha5
+ m6 * Malpha6
+ Malpha7;
expand(p);
m(X) := m0 + m1*X + m2*X^2 + m3*X^3 + m4*X^4 + m5*X^5 + m6*X^6 + X^7;
invert(Malpha);
|
