05/04/09 15:39:48
>>329 最近 mAXIMAも6。0目指していろいろ主要部分で変化の
最中です。そのなかで gcdのmainのアルゴリズムがspmodを使うように
なるようです。そこで これを使うとalgsysのbugの一部がとれて
>329で指摘されていたH[4,4]のjordan標準型が数式で求まります。
(%i1) load("diag")$
(%i2) H[I,J]:=1/(I+J-1)$
(%i3) gcd:'spmod$ ;これが肝心
(%i4) display2d:false$
(%i5) jordan(GENMATRIX(H,4,4));
(%o5) [[-SQRT(-((357210000*(SQRT(4621318097)*%I/(750141000*SQRT(3))
+3805076179/6751269000000)
^(2/3)
-375521400*(SQRT。。。。。。。以下略
こうやってゆっくりですが着実に改良されているようです。