\\ selfpoint on X_0(11) e = ellinit("11a1"); polx = 25937424601*x^12 - 48202871557476252*x^11 - 627575688471844224310*x^10 - 4586587322380883649178756*x^9 + 189861110415625174383936023*x^8 - 9450307537215069858510760088*x^7 - 13176917774298176001346511796*x^6 + 558725269921007151668021541368*x^5 + 4102246055136279443069069843703*x^4 + 12745046520678761793965586279924*x^3 + 25688385347866823866818727620042*x^2 + 32963462539585829067954484001996*x + 26652074520418260289484123648825 ; xP = x; yP = -198233393415427385924401628684747483074876904559199671167318314586767011062493046/242723107073930774051072510458839219094317233689099365731320478482825253065409535837152287304200772071*x^11 + 9555583792625036082239903192557732249197142256403002590619439428064758991284895427816259018792158/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^10 + 124099832906101711758861444450374027231886837186185598696228279586253137933012329129851750868239851090/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^9 + 905223693647020574535250413840306131809639213064414031893398407577880420991147727633534741394200138650431/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^8 - 66799465785132133735843552836139924160153883851693211746092673681448348697480552683157056379866551776228976/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^7 + 3177845132247246441790637859097852868318834582669305711046948762425650402823668481689289563927006138746216413/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^6 - 61827641699114968240160439478511497751632723479154709970327520176533225296800194124080152850054083268388659346/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^5 + 19366095660038259533857287051809736336912056708912969155539542245956972865136613322001431148742059956843143730/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^4 + 2295653115877889622262896746430222066310696653524268761038482085434057474355944096075921626896885337451761711957/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^3 + 14658886787119486209237448846282661147416376944845689736583659906699777248351554377295198763920107049645754476923/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x^2 + 30853850084165645138413366308778125989983455310969408946747615197702348430656552183837788456758618117505823865102/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671*x + 27738247126182964431456406705411984433495323665504426017895623943801889286631093493428584633333943631122309774299/6295612288650529184643260563209926159240572756385911874253048134815522874842803956962544866507397076157549118671; P = Mod( [xP,yP], polx ); ellisoncurve(e,P)