円分多項式の例です． (%i9) ROU(x,N):=rat(apply("*",makelist((x^(N/s)-1)^moebius(s),s,listify(divisors(N)))))$ Evaluation took 0.0000 seconds (0.0000 elapsed) using 0 bytes. (%i10) for i:3 thru 22 do print([p:ROU(x,i),mp(p),RR6(cs(nGG(DA)))…

現在のオリジナルプログラムでは，すべての Lagrange Resolvent の冪乗，開冪，和を取って，表題の２つを算出しています．一方，http://ehito.hatenablog.com/entry/20190123/1548255005 で提案した方法は，分解体の（基礎体 K 上の）定義式 P，新たに添加す…