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$