Maxima で綴る数学の旅

東池袋 Adomani 谷山志村対応をいくつかの楕円曲線で確認しました。その際に与えられた楕円曲線の法pでの解の個数を総当たりで求めていました。これを1000倍くらい高速化してみましょう。 (%i1) Nsolve(elc,p):=block([c:0,evelc], for x:0 while x