(4) O
(5) O
(1)
O O
(9) O
O
(2)
(7)
O
(8) (6) (3)
(10) O
O
O
O O
restart:
read "qFPS.mpl":
with(qFPS):
Sitzung 01 (q-Ableitungen und q-Shifts)
q-Ableitung bzw. q-Shift von Potenzen qdiff(x^k,x,q);
qkK1 xkK1 qK1 qshift(x^k,x,q);
q x k
Zweite q-Ableitung bzw. doppelter q-Shift des q-Pochhammersymbols qdiff(qpochhammer(x,q,k),[x$2],q);
K qkK1 KqkCq qpochhammer x,q,k qK1 2 K1Cx K1Cq x qshift(qpochhammer(x,q,k),[x$2],q);
K1Cq x qk K1Cx qk qpochhammer x,q,k K1Cq x K1Cx
q-Ableitung bzw. q-Shift der kleinen q-Exponentialfunktion qdiff(qexp(x,q),x,q);
Kqexp x,q qK1 qshift(qexp(x,q),x,q);
1Kx qexp x,q qdiff(qexp((1-q)*x,q),x,q);
qexp 1Kq x,q Inverser q-Shift
qshift(qpochhammer(x,q,k),x,1/q);
KqCx qpochhammer x,q,k KqCx qk
Kompliziertere, höhere q-Ableitungen
qdiff(sinq(x,q)*f(x,y)+cosq(y*x,q)*g(y),[x,y,y],q);
Ky2 q2 x qK1 qC1 cosq y x,q Dqy,y g y Cy q2 y2 q3 x2K2 q2 y2 x2Cy2 q x2 K1 Dqy,y g y sinq y x,q C qK1 x cosq x,q Dqx,y,y f x,y Csinq x, q Dqx,y,y f x,y Ccosq x,q Dqy,y f x,y Kx qC1 cosq y x,q g y
Cy x2 q2 sinq y x,q g y C qC1 y2 q3 x2Kq2 y2 x2K1 Dqy g y sinq y x,q Ky x qC1 q2CqK1 cosq y x,q Dqy g y
Sitzung 02 (q-Holonome Differential- und Rekursionsgleichungen)
Beispiele q-holonomer Differential- und Rekursionsgleichungen q-holonomer Funktionen Potenzen
qHolonomicDE(x^k,F(x));
O O
(15) O
O O
O
(20) O
(16)
O
(18)
O
(11)
(12)
(14) O
(21) (13) (10)
O
(19) (17) qK1 x Dqx F x C KqkC1 F x = 0
qHolonomicRE(x^k,F(x));
F x qkKSqx F x = 0 Kleine q-Exponentialfunktion
qHolonomicDE(qexp(x,q),F(x));
qK1 Dqx F x CF x = 0 qHolonomicRE(qexp(x,q),F(x));
Sqx F x C K1Cx F x = 0 Große q-Cosinusfunktion
qHolonomicDE(qCos(x,q),F(x));
q x qK1 qC1 Dqx F x C qK1 2 1Cq2 x2 Dqx,x F x Cq F x = 0 qHolonomicRE(qCos(x,q),F(x));
K1Kq Sqx F x Cq F x C 1Cq2 x2 Sqx,x F x = 0 q-Orthogonales Polynom: q-Laguerre
qHolonomicDE(qLaguerreL(n,a,x,q),F(x));
qK1 qa q2 xKq qn qa xCqa qCq qa xK1 Dqx F x Cx qa q qK1 2 q x C1 Dqx,x F x Kqa q qnK1 F x = 0
qHolonomicRE(qLaguerreL(n,a,x,q),F(x));
Kq qn qa xKqaK1 Sqx F x CF x Cqa q xC1 Sqx,x F x = 0
Sitzung 03 (q-Holonome Rekursionsgleichung für die verallgemeinerte q- hypergeometrische Funktion)
qHolonomicRE(qphihypergeom([a,b,c],[d,e],x,q),F(x));
q q x cCq x bCq x aKqKeKd Sqx F x Kq2 K1Cx F x C x q2 a b c Kd e Sqx,x,x F x C Kx q2 b cCd eKx q2 a cKx q2 a bCq e
Cd q Sqx,x F x = 0
Sitzung 04 (Summen-, Produkt und Kompositionsalgorithmus)
RE1:=qHolonomicRE(qsin(x,q),F(x));
RE1:= K1Kq Sqx F x Cq x2C1 F x CSqx,x F x = 0 RE2:=qHolonomicRE(qpochhammer(x,q,k),F(x));
RE2:= K1Cx Sqx F x C 1Kx qk F x = 0 qSumRE(RE1,RE2,F(x));
q K1CxK2 qCq2 x3Cq xKq2 x qk 2Cx q2 qkCq x qkCqkKq x2K2 q2 x2Kq3 x4 Kq4 x4Cq4 x3Cx3 q5Kq4 x2Kq5 x4Cq5 x5C3 q3 xK qk 2C2 x q2K2 q2
Kq3 x qk 2Cx q qk 3Cq4 x3 qk 2K2 q3 x2Cx q4Kq5 x2Kq3C2 qk q Kq2 qk 2C2 qk q2Cqk q3Kq qk 2C2 q3 x3Kq x qk 2Kx2 qk q3Kx2 qk q4 Cx q2 qk 3Kx2 q2 qk 2 Sqx F x Kq3 x2C1 K1Cx qk q3 x2Kx q2Kq x CqKqk qKqkC qk 2C1 F x C K1Cx q2 q x2KxKq xC1CqC qk 2
(29) (24)
O O O
(31) (26) (27) O
(25)
O O
(28) O
(22)
O
(30) (21)
O
O
(10)
(23) O
Kqk qKqk Sqx,x,x F x C 1KxC2 qK2 q xKqkCq x2C2 q2 x2Kq4 x3
Cq4 x2K3 q3 xC qk 2K2 x q2C2 q2C2 q3 x2Cx qk q4K2 x q4Cq5 x2Cq3 K2 qk qCq2 qk 2K2 qk q2Kqk q3Cq qk 2Kq3 x3Cx2 qk q3Cx2 qk q4Kx3 qk q5 Cx qk q3Kx q2 qk 3 Sqx,x F x = 0
qProductRE(RE1,RE2,F(x));
KqC1 K1Cx K1Cq x qk Sqx F x Cq x2C1 K1Cx qk K1 Cq x qk F x C K1Cq x K1Cx Sqx,x F x = 0
qCompositionRE(RE1,F(x),a*x^2);
1Cq2 q3 a2 x4K1 Sqx F x Cq2 1Ca2 x4 q2 a2 x4C1 F x CSqx,x F x
= 0
qHolonomicRE(qsin(a*x^2,q),F(x));
1Cq2 q3 a2 x4K1 Sqx F x Cq2 1Ca2 x4 q2 a2 x4C1 F x CSqx,x F x
= 0
Sitzung 05 (q-Petkovšek-Algorithmus)
RE1:=(q^(k+2)-1)*A(k+2)+(q^(2*k+2)*(1+q)-q^(k+1))*A(k+1)+q^(3*
k+2)*A(k)=0;
RE1:= qkC2K1 A kC2 C q2C2 k qC1 KqkC1 A kC1 Cq2C3 k A k = 0 qPetkovsek(RE1,A(k));
K1 k qbinomialk, 2 , qk2
qpochhammer q,q,k
RE2:=A(k+2)-(1+q)*A(k+1)+q*(1-q^(2*k+1))*A(k)=0;
RE2:=A kC2 K qC1 A kC1 Cq 1Kq1C2 k A k = 0 qPetkovsek(RE2,A(k));
qPetkovsek(RE2,A(k),quadratic);
qpochhammer K q,q,k ,qpochhammer q,q,k
Sitzung 06 (q-FPS-Algorithmus)
convert(qexp(x,q),qFPS);
k= 0
>
N xk
qpochhammer q,q,k infolevel[qFPS]:=4:
convert(qsin(x,q),qFPS);
convert/qFPS: q-holonomic recurrence equation of order 2 KqK1 Sqx F x Cq x2C1 F x CSqx,x F x = 0
convert/qFPS: q-holonomic recurrence equation for series coefficients of order 2
A k C qk q2K1 qk qK1 A kC2 = 0 convert/qFPS: solution of recurrence equation
0, K1 k x1C2 k qpochhammer q,q, 1C2 k
O O
(31)
O
(34) O
O
(37) (32) O
O
(21)
(36) O
(10)
(35) (33)
k= 0
>
N K1 k x1C2 k
qpochhammer q,q, 1C2 k infolevel[qFPS]:=0:
convert(qCos(x^2,q),qFPS);
k
>
= 0N K1 k qk K1C2 k x2 2 k qpochhammer q,q, 2 k
convert(qphihypergeom([a,b],[c,d],x,q),qFPS,x);
k= 0
>
N K1 k q
1
2 kK1 k
qpochhammer a,q,k qpochhammer b,q,k xk qpochhammer c,q,k qpochhammer d,q,k qpochhammer q,q,k Anwendungsbeispiel: Beweis von q-Identitäten
convert(qsin(x,1/q),qFPS);
k= 0
>
N
K K1 k q kC1 1C2 k x1C2 k qpochhammer q,q, 1C2 k convert(-qSin(q*x,q),qFPS);
k= 0
>
N
K K1 k q kC1 1C2 k x1C2 k qpochhammer q,q, 1C2 k oder
convert(x/(1-q)*qphihypergeom([0,0],[q^3],-x^2,q^2),qFPS);
k= 0
>
N K1 k x1C2 k
qpochhammer q,q, 1C2 k convert(qsin(x,q),qFPS);
k= 0
>
N K1 k x1C2 k qpochhammer q,q, 1C2 k