blob: 317fc72461fcefa33ff87c43bd063ad401304906 (
plain)
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
|
% Copyright The Numerical Algorithms Group Limited 1992-94. All rights reserved.
% !! DO NOT MODIFY THIS FILE BY HAND !! Created by ht.awk.
\newcommand{\LieExponentialsXmpTitle}{LieExponentials}
\newcommand{\LieExponentialsXmpNumber}{9.42}
%
% =====================================================================
\begin{page}{LieExponentialsXmpPage}{9.42 LieExponentials}
% =====================================================================
\beginscroll
\xtc{
}{
\spadpaste{ a: Symbol := 'a \bound{a}}
}
\xtc{
}{
\spadpaste{ b: Symbol := 'b \bound{b}}
}
Declarations of domains
\xtc{
}{
\spadpaste{ coef := Fraction(Integer) \bound{coef}}
}
\xtc{
}{
\spadpaste{ group := LieExponentials(Symbol, coef, 3) \free{coef} \bound{group}}
}
\xtc{
}{
\spadpaste{ lpoly := LiePolynomial(Symbol, coef) \free{coef} \bound{lpoly}}
}
\xtc{
}{
\spadpaste{ poly := XPBWPolynomial(Symbol, coef) \free{coef} \bound{poly}}
}
Calculations
\xtc{
}{
\spadpaste{ ea := exp(a::lpoly)$group \free{a} \free{lpoly} \free{group} \bound{ea}}
}
\xtc{
}{
\spadpaste{ eb := exp(b::lpoly)$group \free{b} \free{lpoly} \free{group} \bound{eb}}
}
\xtc{
}{
\spadpaste{ g: group := ea*eb \free{ea} \free{eb} \bound{g}}
}
\xtc{
}{
\spadpaste{ g :: poly \free{g} \free{poly}}
}
\xtc{
}{
\spadpaste{ log(g)$group \free{g} \free{group}}
}
\xtc{
}{
\spadpaste{ g1: group := inv(g) \free{g} \free{group} \bound{g1}}
}
\xtc{
}{
\spadpaste{ g*g1 \free{g} \free{g1}}
}
\endscroll
\autobuttons
\end{page}
%
|
