bash-2.05$ gp GP/PARI CALCULATOR Version 2.1.3 (released) i386 running netbsd 32-bit version (readline v1.0 enabled, extended help available) Copyright (C) 2000 The PARI Group PARI/GP is free software, covered by the GNU General Public License, and comes WITHOUT ANY WARRANTY WHATSOEVER. Type ? for help, \q to quit. Type ?12 for how to get moral (and possibly technical) support. realprecision = 28 significant digits seriesprecision = 16 significant terms format = g0.28 parisize = 4000000, primelimit = 500000 ? e=ellinit([0,0,0,877,0]) %1 = [0, 0, 0, 877, 0, 0, 1754, 0, -769129, -42096, 0, -43169672512, 1728, [0.E-28, 0.E-28 - 29.61418578992169421552065349*I, 0.E-28 + 29.61418578992169421552065349*I]~, 0.6814081438510085906928801662, -0.3407040719255042953464400831 + 0.3407040719255042953464400831*I, -4.610441894390844665097696094 + 5.24064267 E-29*I, 2.305220947195422332548848047 - 6.915662841586266997646544141*I, 0.2321585292532384083231903681] ? elltors(e) %2 = [2, [2], [[0, 0]]] ? x0=29604565304828237474403861024284371796799791624792913256602210; ? y0= -256256267988926809388776834045513089648669153204356603464786949; ? z0=490078023219787588959802933995928925096061616470779979261000; ? x=x0/z0 %6 = 375494528127162193105504069942092792346201/6215987776871505425463220780697238044100 ? y=y0/z0 %7 = -256256267988926809388776834045513089648669153204356603464786949/490078023219787588959802933995928925096061616470779979261000 ? ellisoncurve(e,[x,y]) %8 = 1 ? factor(x) %9 = [2 -2] [3 -2] [5 -2] [37 -2] [41 2] [53 -2] [577 -2] [4013 -2] [4337 -2] [9013 2] [133451 -2] [1658244549709897 2] ? factor(y) %10 = [-1 1] [2 -3] [3 -3] [5 -3] [37 -3] [41 1] [53 -3] [73 1] [577 -3] [4013 -3] [4337 -3] [9013 1] [133451 -3] [54629473 1] [481273108222873 1] [1658244549709897 1] [217886952176630497 1] ? p2=ellpow(e,[x,y],2) %11 = [11472589076204841027503850994111129000443596097387106034430620664607790165733763259523861305965114273599875048594705071554210891246465281851434443352179121788135284801/1632747912067990387379510802200492457676725500764669551203697787883829890478413680968072412463676726939085243874796602136393541306114641069922145227273979088677216400, -5322830902999433189837600469352665585715151486153939512041083501445401674611807112151153172975608564285609598366234998016624610977555052048815356867068768609876316734461567620130105953802679615455708438262502174018110970592272694328561860040502527201/65974894679425281678257769430049247939927922262556891824912946365743855271877854039314720909046654276650455589303159739423962721223940611656138453327073352584276122278073688069387361086220666853388582093955892979296028764398060033735270526960712000] ? p3=ellpow(e,[x,y],3) %12 = [56178585520458416830033416237087819970312353993896197856400764697432647423389522874301560531878358591840172908333793709003904282270565151237774766848249596741241450555274408800114286737516703169571371590130085270450736768748507576447865367783169820371866176150201151001566776468334399920489437334734502082992259669346076941299899125079070987765741574192877243272468566515801/47220095524904932780517720012408412715266114775806859040431913120920988616264161305155784394735817570013999549869128649481389366796117177512401699608106475491768020016705731046051525123659113454754207472285871106396435059322309064266894756291295738119601698625533073271577210136975400700356436515631882325039390888861066236313091506967373002089248916389546160097876416636900, 10489679980743678888087859698961841969187004534495848349781201947808154481515731728069068091281742936214293401721446851332523213466660627423603063779858672629119462505967087361633092643180433984786450915940938203443423937163085634084014748325209428465006469283502098553575490814709611434155014009390147612675058094515970638562358392242550958626371961402962815043963976266840023910110795093895058381207523261778688478512349410750353612223396467624522430758028270633296543032755097642448927264081247728248422456299946900158426203900350218520842801702545453634467451/324481762250319636036378354029343581373906393847660669951534425682218861702035729360246416401783716029680085665859474254073716956480126331676081216306048214224885948055668985855545892414865263034650721808168698884893096598295053905106108277948100943214818695073192810487358554312926137229638545099082383854247432285502067787646070216491928208988606387879992606719408869509084413873610160936867507916889305031718119511950851303524241524471898792422592733414328642496150601095843732848866565231327613366937278772600471202895353334072662611123493728462991916847000] ? elladd(e,[x,y],[0,0]) %13 = [5451421280316310258131244624671477764675700/375494528127162193105504069942092792346201, 28915277188538778175376314749786089767316948966312631462931224330/230094066204268982894707440563804275878997592508073394775934301] ? p=[x,y] %14 = [375494528127162193105504069942092792346201/6215987776871505425463220780697238044100, -256256267988926809388776834045513089648669153204356603464786949/490078023219787588959802933995928925096061616470779979261000] ? elladd(e,p2,[0,0]) %15 = [1431919918883627569731830973529831885382488264170615196405642959974118813949568798208999505730644489525577758878196620073617135725462540218321721364319279660769918782800/11472589076204841027503850994111129000443596097387106034430620664607790165733763259523861305965114273599875048594705071554210891246465281851434443352179121788135284801, 1761047670992384685442382450333620397450317334743109390666315846811160849538674684335721908202897725258272811677367544961041851826817251568074293495639181563826661270606307763783576224741536863346025059445730521833800865047119205692512943934306246962660/1228830996217737588102287590776060554318515371622242754020875517655279515292698418912799954150262193876347027814356953913293473891028511946082343726330374216468788285955831007298823470487124924534483184069626885460470725368115288680514419388707567201] ? e.j %16 = 1728 ? e.disc %17 = -43169672512 ? ellglobalred(e) %18 = [49224256, [1, 0, 0, 0], 2] ? quit; Good bye! bash-2.05$