この説明書では 1変数代数的局所コホモロジー類用のパッケージ taji_alc について説明する. 数学的解説や背景については, 解説記事 "1変数代数的局所コホモロジー類用に対する Risa/Asir 用パッケージ taji_alc" (Risa/Asir Journal (2007)) およびその参考文献を参照.
本セクションの関数を呼び出すには,
import("taji_alc.rr")$
を実行してプログラムをロードする.
taji_alc.cpfd
[235] taji_alc.cpfd(x^3-x-1,x^4+2*x^3+2*x^2+2*x+1); [[[1/2*x-1,[x^2+1,1]]],[[-1/2,[x+1,2]],[1/2,[x+1,1]]]] [236] taji_alc.cpfd(x^3-x-1,x^4+2*x^3+2*x^2+2*x+1|switch=1); [[[[x-2,2],[x^2+1,1]]],[[[-1,2],[x+1,2]],[[1,2],[x+1,1]]]] [237] taji_alc.cpfd(x^3-x-1,x^4+2*x^3+2*x^2+2*x+1|switch=10); [[1/2*x-1,[x^2+1,1]],[1/2*x,[x+1,2]]] [238] taji_alc.cpfd(x^3-x-1,x^4+2*x^3+2*x^2+2*x+1|switch=11); [[[x-2,2],[x^2+1,1]],[[x,2],[x+1,2]]]
ChangeLog
taji_alc.snoether
[296] taji_alc.snoether(1,[[x^3-x-1,3]]); [[x^3-x-1,[9/529*x^2-27/1058*x+11/1058,-81/529*x^2-9/529*x+135/529,-49 05/12167*x^2+4563/12167*x+3270/12167]]] [299] taji_alc.snoether(1,[[x^3-x-1,3]]|switch=1); [[x^3-x-1,[[18*x^2-27*x+11,1058],[-81*x^2-9*x+135,529],[-4905*x^2+4563 *x+3270,12167]]]] [297] taji_alc.snoether(1,[[x^3-x-1,3]]|switch=10); [[x^3-x-1,[[414*x^2-621*x+253,-3726*x^2-414*x+6210,-9810*x^2+9126*x+65 40],24334]]] [298] taji_alc.snoether(1,[[x^3-x-1,3]]|switch=20); [[x^3-x-1,[[[18*x^2-27*x+11,529],[-162*x^2-18*x+270,529],[-9810*x^2+91 26*x+6540,12167]],2]]] [241] taji_alc.snoether(x^3+1,x^18-2*x^14+x^10-x^8+2*x^4-1|switch=10); [[x^4+x^3+x^2+x+1,[[-2*x^2-x-2],50]],[x^4-x^3+x^2-x+1,[[-2*x^3+4*x^2-x -2],50]],[x^2+1,[[-x+1,8*x+5],32]],[x+1,[[-6,-39],320]],[x-1,[[2,-24,6 7],320]]]
ChangeLog
taji_alc.laurent_expansion
[354] taji_alc.laurent_expansion(x,(x-1)^3); [[x-1,[1,1,0]]] [356] taji_alc.laurent_expansion(x^5+x^4+x^3+x^2+x+1,(x^4+1)^3); [[x^4+1,[1/64*x^2+1/64*x,1/16*x^3+1/16*x^2-3/128*x-5/128,-5/128*x^3-1/ 8*x^2-3/16*x]]]
taji_alc.snoether
ChangeLog
taji_alc.residue
[219] taji_alc.residue(1,x^4+1); [[x^4+1,-1/4*x]]この例で言うと, 求めた留数多項式-1/4*xに, x^4+1の(4つある)零点をそれぞれ代入したものが個別の留数値である.
[221] taji_alc.residue(x^8,[[x^3-x-1,3]]); [[x^3-x-1,-2243/12167*x^2+2801/12167*x+5551/12167]] [222] taji_alc.residue(x^2+x,[[x+1,3],[x-1,3],[x^2+3*x-1,2]]); [[x^2+3*x-1,-284/4563*x-311/1521],[x-1,89/432],[x+1,7/432]] [223] taji_alc.residue(x^2+x,[[x+1,3],[x-1,3],[x^2+3*x-1,2]]|switch=1) ; [[x^2+3*x-1,[-284*x-933,4563]],[x-1,[89,432]],[x+1,[7,432]]] [234] taji_alc.residue(x^2+x,[[x+1,3],[x-1,3],[x^2+3*x-1,2]]|switch=1, pole=[x+1]); [[x+1,[7,432]]] [225] taji_alc.residue(x^3+1,x^18-2*x^14+x^10-x^8+2*x^4-1); [[x^4+x^3+x^2+x+1,-1/25*x^2-1/50*x-1/25],[x^4-x^3+x^2-x+1,-1/25*x^3+2/ 25*x^2-1/50*x-1/25],[x^2+1,1/4*x+5/32],[x+1,-39/320],[x-1,67/320]]
ChangeLog
taji_alc.invpow
[236] taji_alc.invpow(3*x^2-1,x^3-x-1,1); -6/23*x^2+9/23*x+4/23 [237] taji_alc.invpow(3*x^2-1,x^3-x-1,1|switch=1); [-6*x^2+9*x+4,23] [238] taji_alc.invpow(3*x^2-1,x^3-x-1,30|switch=1); [1857324483*x^2-2100154824*x-477264412,266635235464391245607]
ChangeLog
taji_alc.rem_formula
[228] taji_alc.rem_formula([[x^2+1,1,z]]); [[-1/2*z*x+1/2]]の正しい見方は, x^2+1の零点をa1,a2とおいたときに, zにa1とa2を代入した, r(x)=(-1/2*a1*x+1/2)+(-1/2*a2*x+1/2) である. しかし出力では, 零点の和の部分を便宜上省略して返す.
[583] taji_alc.rem_formula([[x-1,1,z1],[x-2,1,z2]]); [[-x+2],[x-1]] [584] taji_alc.rem_formula([[x-1,1,z1],[x-2,1,z2]]|switch=20); (-p^(0)(z1)+p^(0)(z2))*x+2*p^(0)(z1)-p^(0)(z2) [587] taji_alc.rem_formula([[x-1,2,z1]]); [[x-1,1]] [588] taji_alc.rem_formula([[x-1,2,z1]]|switch=20); p^(1)(z1)*x-p^(1)(z1)+p^(0)(z1) [494] taji_alc.rem_formula([[x-1,3,z1]]|switch=20); 1/2*p^(2)(z1)*x^2+(-p^(2)(z1)+p^(1)(z1))*x+1/2*p^(2)(z1)-p^(1)(z1)+p^( 0)(z1) [229] taji_alc.rem_formula([[x+1,2,z1],[x^3-x-1,1,z2]]); [[-x^4-x^3+x^2+2*x+1,-2*x^4-3*x^3+2*x^2+5*x+3],[(-1/23*z2^2-10/23*z2+1 6/23)*x^4+(-12/23*z2^2-5/23*z2+31/23)*x^3+(-5/23*z2^2+19/23*z2-12/23)* x^2+(22/23*z2^2+13/23*z2-53/23)*x+16/23*z2^2-1/23*z2-26/23]] [230] taji_alc.rem_formula([[x+1,2,z1],[x^3-x-1,1,z2]]|switch=20); (-1/23*p^(0)(z2)*z2^2-10/23*p^(0)(z2)*z2-2*p^(0)(z1)+16/23*p^(0)(z2)-p ^(1)(z1))*x^4+(-12/23*p^(0)(z2)*z2^2-5/23*p^(0)(z2)*z2-3*p^(0)(z1)+31/ 23*p^(0)(z2)-p^(1)(z1))*x^3+(-5/23*p^(0)(z2)*z2^2+19/23*p^(0)(z2)*z2+2 *p^(0)(z1)-12/23*p^(0)(z2)+p^(1)(z1))*x^2+(22/23*p^(0)(z2)*z2^2+13/23* p^(0)(z2)*z2+5*p^(0)(z1)-53/23*p^(0)(z2)+2*p^(1)(z1))*x+16/23*p^(0)(z2 )*z2^2-1/23*p^(0)(z2)*z2+3*p^(0)(z1)-26/23*p^(0)(z2)+p^(1)(z1) [231] taji_alc.rem_formula([[x^3-x-1,2,z]]|switch=10); [[[(3/23*z^2-4/23)*x^2+(-1/23*z+3/23)*x-4/23*z^2+3/23*z+4/23,(162/529* z^2-174/529*z-108/529)*x^2+(-105/529*z^2+54/529*z+70/529)*x-108/529*z^ 2+116/529*z+72/529],[(-6/23*z^2+9/23*z+4/23)*x^2+(9/23*z^2-2/23*z-6/23 )*x+4/23*z^2-6/23*z+5/23]]]
ChangeLog
taji_alc.solve_ode_cp
[287] taji_alc.solve_ode_cp(x*(x-3)^2,z,0); [[[x-3,0],[x,1]],[[x-3,-z+2/3],[x,-2/3]],[[x-3,1/3*z-1/9],[x,1/9]]] [289] taji_alc.solve_ode_cp((x^3-x-1)^2,z,0|switch=1); [[[x^3-x-1,[(92*z+200)*x^2+(-69*z-254)*x-92*z+43,529]]],[[x^3-x-1,[(92 *z+420)*x^2+(-46*z-216)*x-161*z-280,529]]],[[x^3-x-1,[(-69*z-195)*x^2+ (23*z+327)*x+23*z+130,529]]],[[x^3-x-1,[(-161*z-270)*x^2+(69*z+290)*x+ 184*z+180,529]]],[[x^3-x-1,[-105*x^2+(-23*z+54)*x+69*z+70,529]]],[[x^3 -x-1,[(69*z+162)*x^2-174*x-92*z-108,529]]]] [277] taji_alc.solve_ode_cp(x^2-4,z,0); [[[x+2,1/2],[x-2,1/2]],[[x+2,-1/4],[x-2,1/4]]] [278] taji_alc.solve_ode_cp(x^2-4,z,0|data=[1,-1]); [[x+2,3/4],[x-2,1/4]] [279] taji_alc.solve_ode_cp(x^2-4,z,0|data=[c0,c1]); [[x+2,1/2*c0-1/4*c1],[x-2,1/2*c0+1/4*c1]]
taji_alc.solve_ode_cp_ps
[345] taji_alc.solve_ode_cp_ps((x-2)*(x+3),z,[[x-1,1]]); [[x+3,1/20],[x-1,-1/4],[x-2,1/5]] [346] taji_alc.solve_ode_cp_ps((x-2)*(x+3),z,[[x-1,1]]|switch2=1); [[x-1,-1/4]] [347] taji_alc.solve_ode_cp_ps((x-2)*(x+3),z,[[x-2,1]]); [[x+3,1/25],[x-2,1/5*z-1/25]] [348] taji_alc.solve_ode_cp_ps((x-2)*(x+3),z,[[x-2,1]]|switch2=1); [[x-2,1/5*z-1/25]] [349] taji_alc.solve_ode_cp_ps((x-2)*(x+3),z,[[x+1,1],[x-2,1]]|switch2 =1); [[x+1,-1/6],[x-2,1/5*z+2/75]] [350] taji_alc.solve_ode_cp_ps((x^3-x-1)*(x-3)^2,z,[[x-3,2],[x-1,3*z^2 +1]]); [[x-1,[-6*z^2-36*z-119,8]],[x^3-x-1,[42291*x^2+55504*x+32313,12167]],[ x-3,[4232*z^2-4278*z-4295,97336]]]
ChangeLog
taji_alc.fbt
[235] taji_alc.fbt(1,(x^3-x-1)^3,z); [[x^3-x-1,(9/529*z^2-81/529*z-4905/12167)*x^2+(-27/1058*z^2-9/529*z+45 63/12167)*x+11/1058*z^2+135/529*z+3270/12167]]
taji_alc.residue
ChangeLog
taji_alc.invfbt
[8] taji_alc.invfbt([[x^3-x-1,2*x^2*z^2+x*z+1],[x^2+1,z*x+z^2]],z|swit ch=1); [3*x^14+14*x^12+39*x^11+33*x^10+179*x^9+206*x^8+350*x^7+223*x^6+126*x^ 5+176*x^4+107*x^3+101*x^2+15*x-4,[[x^2+1,3],[x^3-x-1,3]]] [9] taji_alc.fbt(3*x^14+14*x^12+39*x^11+33*x^10+179*x^9+206*x^8+350*x^ 7+223*x^6+126*x^5+176*x^4+107*x^3+101*x^2+15*x-4,[[x^2+1,3],[x^3-x-1,3 ]],z); [[x^3-x-1,2*z^2*x^2+z*x+1],[x^2+1,z*x+z^2]]
taji_alc.fbt
ChangeLog
Jump to: t
@vfill @eject
This document was generated on 23 February 2009 using texi2html 1.56k.