/* the initial points F^(-k)[1], k = 0, .., 10 . Value at x=1 {1.96727*10^-46, 1.59598*10^-79, 3.3197*10^-77, -3.3217*10^-75, -9.66052*10^-74, 9.06082*10^-72, 4.68389*10^-70, -4.0926*10^-68, 1.52713*10^-67, 1.00167*10^-65, -1.32142*10^-63}$ non-negative singularities: 0, 0.31963094681032395077983870774802840706, 0.67081712157165105809184801420320961377, 1.1864243344575629425373042576712194034, 1.4582281371202052932900597674187667223, 1.7005805566619710685261523441444767466, 17.274073285432277196363444225596156561, 23.14654015969242635782942325005611441 */ Init0=[ 1.96727*10^(-46), 1.59598*10^(-79), 3.3197*10^(-77), -3.3217*10^(-75), -9.66052*10^(-74), 9.06082*10^(-72), 4.68389*10^(-70), -4.0926*10^(-68), 1.52713*10^(-67), 1.00167*10^(-65), -1.32142*10^(-63)]$ // Init=map(number_float_to_rational,Init0)$ // BUG it returns 0. Todo. def fr(A,B) { return( number_float_to_rational(A)*B ); } Init=[ fr(1.96727,10^(-46)), fr(1.59598,10^(-79)), fr(3.3197,10^(-77)), fr(-3.3217,10^(-75)), fr(-9.66052,10^(-74)), fr(9.06082,10^(-72)), fr(4.68389,10^(-70)), fr(-4.0926,10^(-68)), fr(1.52713,10^(-67)), fr(1.00167,10^(-65)), fr(-1.32142,10^(-63))]$ ODE= -1562091597748518867838547156511013572269431069694503640906180320395189264000000000000000*x^5 - 69674715609195743031882560782396374958263663460914884212693609922879579077531200000000000000*x^7 + 775918585957785223397020721300516010930029525700364979740408742014061646629159561000000000000000*x^9 - 22670269590333702664811416061282933611346014242381934943442424385796574964613491669970000000000000000*x^11 + 2609652889608897959554878107778038318866158338248494506110148535081362030048106492200500000000000000000*x^13 - 514110153837998293738330478190806388113190716056326883498219647553278691105293008385225000000000000000000*x^15 + 8442288873899810813737326358691427494463505807238877983746167654958791468625943595475000000000000000000000*x^17 - 55501404010693428072492630370395958813396008779096782798737657863160283569082639356875000000000000000000000*x^19 + 290159083025648634833823749354869137356703113051229869439568526730994675479701328500000000000000000000000000*x^21 - 1190462932642501765842007522280203394864952588548821277975078925376325673117943462500000000000000000000000000* x^23 + 3632790520546052032150757532258931869347721391770915006441034215241703925353331250000000000000000000000000000*x^25 - 9023129919646260165495902319649609385621305660047643033041346868065054435781250000000000000000000000000000000*x^27 + 17524728112880223067977233734995474108768438093218019262737199454421374465625000000000000000000000000000000000* x^29 - 22680865935259028854332457515640233971485627678132855457142957647176556250000000000000000000000000000000000000*x^31 + 17147894004250732620655137011173176003661058710759272217491710958678125000000000000000000000000000000000000000*x^33 - 6733925569390330400094236080214393978981059354706876627069865818750000000000000000000000000000000000000000000* x^35 + 1176331135025279019247262846701678522304985394964106777165200000000000000000000000000000000000000000000000000*x^37 - 228088107171758240237841675911515168849370520627778673550000000000000000000000000000000000000000000000000000*x^39 + 150980451180910435797134789254875293059305639948935100000000000000000000000000000000000000000000000000000000* x^41 - 38463028105340849673134439059728733131004093215125000000000000000000000000000000000000000000000000000000000*x^43 + 2253750558467716287200691513406234683275087250000000000000000000000000000000000000000000000000000000000000*x^45 - 10964475087059478190588339764957995790000000000000000000000000000000000000000000000000000000000000000000*x^47 + 12861785248637785542190702222200000000000000000000000000000000000000000000000000000000000000000000000*x^49 + (12468659250134046454123175769915778729979183359641632852294681950925275200000 + 1433079504675018002419078398156153341950230151804175953670150656002402565920000000*x^2 + 43967691489430841340128991870723500746842456962513577366699342088954897884009600000000*x^4 - 631126985050960985563812220401166981635982463666648177675987863319458131708921534000000000*x^6 + 14748929025234562850593141181352475879141152799612547470322292789056602078128161031800000000000*x^8 - 14399411538212163010661889744126954209215408227792554231934682028007780057533918170327500000000000*x^10 + 9089387665252109477795946585776716323715783146689108532836029748274520187369116264704000000000000000*x^12 - 906752777300392672795985459101246697340245394384662056302404987559680943136322341399606250000000000000*x^14 + 19055113473565495366394834497665161805503047614853461464569515864874489525711160887933750000000000000000*x^16 - 170477037473706911928593301255254361627790955870293889145116932430631901803297756669437500000000000000000*x^18 + 932818360838081240885213393825618452565562095596080960451425530456436486847604262200000000000000000000000*x^20 - 4300421803936860071272390067893301933656627602177793709724946752722692689594899546875000000000000000000000*x^22 + 14945842724584574759012959590534012895304237994226792874683427224022048225312356250000000000000000000000000*x^24 - 39999057335179730268665796626378375655803379534855629429367122242978326391123828125000000000000000000000000* x^26 + 88261256456627602350581480990861036860076848762317155404685433268941116757812500000000000000000000000000000*x^28 - 140177182026412451622676617416968787291376472274200105321635858193489619726562500000000000000000000000000000*x^30 + 134222096585516222427278119937913084280374604169182028146538317600586132812500000000000000000000000000000000*x^32 - 68932123012353452632716111464359543553506706542589438143858773527343750000000000000000000000000000000000000*x^34 + 16580836631404912615510030560401331379855137493658051356239446875000000000000000000000000000000000000000000* x^36 - 2283016649836073337871812157132320551149338189425054705948437500000000000000000000000000000000000000000000*x^38 + 1223892271592464667582972702490331137331341288055643676562500000000000000000000000000000000000000000000000*x^40 - 485586493374940306097073355608418904789628148564364062500000000000000000000000000000000000000000000000000*x^42 + 58409000375848510690659116053750202880309816750000000000000000000000000000000000000000000000000000000000* x^44 - 294674299627927123772208049990049594273250000000000000000000000000000000000000000000000000000000000000*x^46 + 349665170782284524217466318140810000000000000000000000000000000000000000000000000000000000000000000*x^48)* (dx) + (-12468659250134046454123175769915778729979183359641632852294681950925275200000*x - 1433079504675018002419078398156153341950230151804175953670150656002402565920000000*x^3 - 44742409488635156024805773138353042672683846878351509042209301523771076648009600000000*x^5 + 597048383458698221677069933212559627130504618631284906823516278062913895798040334000000000*x^7 - 14344504147170219183425063022549181511592100196834842694079230070863805036398838477050000000000*x^9 + 2931081575711120623615239125858860906587611861932434348088666414246900971934494493670000000000000*x^11 - 1637358072937440422739336591706133901948528846435401881170676542279097683563355620521625000000000000*x^13 + 41131230958345607974969362123341010910861213976786611286021478751677827537645048758037500000000000000*x^15 - 416336645738781199209336136230107634294781850190471859358399282158759244063302799602500000000000000000*x^17 + 2510403796052383094574318799521531934221112157338145585119768366275827410400264478031250000000000000000*x^19 - 11952682479428012977675099852297417662328871521302514191640024525469735141177550950000000000000000000000*x^21 + 45096617067996686108143449088149997351677437018616391226587442729959435369640250000000000000000000000000*x^23 - 129156087101148628237867969040598555588647861908605438681254170093974259038942187500000000000000000000000*x^25 + 301017394750536094623101839272366583325998877065599599864905355905688127134765625000000000000000000000000*x^27 - 540347951310571169588857470635380935921910068840797295566186074293927453906250000000000000000000000000000*x^29 + 622783496417813596330932248239325937538839786957956995584247050420440039062500000000000000000000000000000*x^31 - 400358836647984968435853586276844301397003094994109683197709632663476562500000000000000000000000000000000*x^33 + 124840350240850175123011441675833123652607904227715353276303753906250000000000000000000000000000000000000*x^35 - 16290380122542861703656216876467539993272505169906150339896875000000000000000000000000000000000000000000*x^37 + 5423045948287948218583593171958946267049453386743997235937500000000000000000000000000000000000000000000*x^39 - 3286452192443232952885315428178676920257796131543726562500000000000000000000000000000000000000000000000*x^41 + 565531275133019048800229697880389264317046460757812500000000000000000000000000000000000000000000000000*x^43 - 2882297787444147992722231142886888285559687500000000000000000000000000000000000000000000000000000000*x^45 + 3430897823784419373164419134277818750000000000000000000000000000000000000000000000000000000000000*x^47)*(dx)^2 + (-59907265997775289751143656785092500602231726929760268797039208666336000 - 3595114683558997781856725774149598077713959052302262825139361716520508000000*x^2 - 108232270547471270244240078561371989430002733403240072160912372082624111208000000*x^4 + 4744331640529022734022756669329514688211842344176222502355617836859319609911870000000*x^6 - 140622926998047751086441750197538641826778829077510221894809442522924872648673708500000000*x^8 + 266534956294788248149443302397392910983630993640394604295780869581263818672818719037500000000*x^10 - 94973158065433024110635985897173741494862444734374022497118016723887205709181346622500000000000*x^12 + 23827873159162439010159114478505539625058413990050457537817594157906087849906213199031250000000000*x^14 - 437616643458700615112751981484630022923203883364393036280762603947100342559999519698125000000000000*x^16 + 3251664921305691175532741998400389035020076441764325471884708371486260756512462102062500000000000000*x^18 - 17077012990082196514080106721449133521171568727746092534320574563217103262039894193750000000000000000*x^20 + 73710772351292850947633550969410704531213946430146127999928283300342034743301834375000000000000000000*x^22 - 235513270415408610385010206450100360147081501006745947485494538154547236135267187500000000000000000000*x^24 + 595904847500033642526313593248464738321022411826894431864569147311693909352343750000000000000000000000*x^26 - 1217184481717826581570947847927727425695186796660634734884980581132431575000000000000000000000000000000*x^28 + 1683105507275273389697590667675296355240547029277272078000378097045663281250000000000000000000000000000*x^30 - 1315514987133547623980634147851045779885667288359865624437542086410546875000000000000000000000000000000*x^32 + 496588409235696119371519876536627709553985149870743154527722120312500000000000000000000000000000000000*x^34 - 67123427432323930675731488766429355168497154969624647113296875000000000000000000000000000000000000000*x^36 + 13024369477832940633640812615058580793904005161855092940625000000000000000000000000000000000000000000*x^38 - 12091567047418273556467379350243275828202488855776118750000000000000000000000000000000000000000000000*x^40 + 2575788733166369697244915689800894644681644835546875000000000000000000000000000000000000000000000000*x^42 - 13187258044989182789734142896139514602879625000000000000000000000000000000000000000000000000000000*x^44 + 15718873934606910955728155359834432500000000000000000000000000000000000000000000000000000000000*x^46)*(dx)^3 + (59907265997775289751143656785092500602231726929760268797039208666336000* x + 7751334433603679933231117697454857654373686838849473775904255700162266400000*x^3 + 274705522147969849298846611605544321029752883627962776368041061189334019848000000*x^5 - 3671714511893404401573332096553762750827428726536419605860827113982649889232670000000*x^7 + 97464775492992676193903148605360362629083867536304581279393305709755810361395088000000000*x^9 - 19422144130231500731223975100351592986692279583994172186836069265760482230331340462500000000*x^11 + 10982593685390803589470243979521928785397422596591911988911259990189406637532389665625000000000*x^13 - 298557299619197519371541649455750397557489781062580156975748176161649416709392342343750000000000*x^15 + 3022581539725496651636697674679881933967215211108932614895301787493558119639079796562500000000000*x^17 - 17024766086372191706415849460021936753181428087291805859075740099112939121209781593750000000000000*x^19 + 80364967595277082429575177999423414887423918434663744343686962257688802847557732812500000000000000*x^21 - 289376828561520501300663645873591522495551894432163386621098446678668981451825000000000000000000000*x^23 + 781577824261159970563422278508411916974244684447822417807942011760005410390234375000000000000000000*x^25 - 1748692116826045688351706312398108255256220222130057640640344044255169481250000000000000000000000000*x^27 + 2847106324182793297604047159990567001893664775865710337523757300239771679687500000000000000000000000*x^29 - 2661162178271522124259280325940229267692001213370669632243191566300195312500000000000000000000000000*x^31 + 1181228042264231242076350843133124754511364778342913869997032185156250000000000000000000000000000000*x^33 - 171480260334156627428434464873965477928193873118351338789368750000000000000000000000000000000000000*x^35 + 17011967919885984751865895906609422018639736728155979685937500000000000000000000000000000000000000*x^37 - 26898250022618918082689133642253193789645551494944967187500000000000000000000000000000000000000000*x^39 + 6748882517552408582705488647497483481938229945914062500000000000000000000000000000000000000000000*x^41 - 34651242165357011112447963128709954069336937500000000000000000000000000000000000000000000000000*x^43 + 41339288489957414687237787813985773750000000000000000000000000000000000000000000000000000000*x^45)*(dx)^4 + (74332013768709211305016614744942839578078739214810986176603235040 - 7119905969482943233680002955586033879381196879618071861049859861984000* x^2 - 211653882254235013514780669247448179396827448554615128164237260983671280000*x^4 - 11615339517349724117947530544880536494368961583970410713246807135111525543550000*x^6 + 437867224119260416921528717288407978794720652228596655807696415660573677625882500000*x^8 - 1050228906668621511974276963673141269799404797471019264324809540639862479670818187500000*x^10 + 304607570021202953561097502931381386617859581270334136649208316548662345735534412500000000*x^12 - 99804581422051782689245552652826018816737445847525882586508219712199720979958402343750000000*x^14 + 1721147425226602316361248839739052779756156936948504780831396839768390769305553940625000000000*x^16 - 11568715737200258632689898523150784889084424725443556903557324726973862897088518906250000000000*x^18 + 58180447019918941524837246852479281885152571823271458292222910873303308889157718750000000000000*x^20 - 240873922595583574771097414897987636606036464678740375837572527271668082964460937500000000000000*x^22 + 708162180406355373623684791881424818361140443590158601943857136928826872554687500000000000000000*x^24 - 1692001855747378095034722412661054836788944510819848373867098909773719089843750000000000000000000*x^26 + 3183068436903463105858372066340243687861510062238231018512890062171242187500000000000000000000000*x^28 - 3537774352161607346211798095550557170409297239265036423719467905863281250000000000000000000000000*x^30 + 1819396763802939624987923528671433428102439905247096219282049703125000000000000000000000000000000*x^32 - 287013301068914249175173558974758360964011463596855052063218750000000000000000000000000000000000*x^34 + 9860173570138289684222139146482039918256968283129796000000000000000000000000000000000000000000*x^36 - 38985911094399363081187952975009635374013094175073500000000000000000000000000000000000000000000*x^38 + 11240744954376951611845467192546533967850748478281250000000000000000000000000000000000000000000*x^40 - 57834668996635281252915150903241204625103750000000000000000000000000000000000000000000000000*x^42 + 69040970938051233295458958611183675000000000000000000000000000000000000000000000000000000*x^44)*(dx)^5 + (-74332013768709211305016614744942839578078739214810986176603235040*x - 12849182696442153350034549306111466321362712096968684404629876360128000* x^3 - 561450042036807053563398218286012884351607367631563944009510993792828080000*x^5 + 7407400857832420600488263857472260054599450291399003287978634422870669183550000*x^7 - 224833366418137524112493749070383327106760035168462975408728496105618095521882500000*x^9 + 43729035294302969949850507554678948910638987242015050118666387492834507782070587500000*x^11 - 24963806516544012006872342719694244927261411385406972599581854333406459239776545000000000*x^13 + 648042459088573497082476726594940722773484710753542103098903781784034579711322218750000000*x^15 - 5997983852851620979099154663125613312332389823910966652582726203610950661158234375000000000*x^17 + 30395409505521590114329297928661104461756579846479777586877561388245121907482031250000000000*x^19 - 141301471339239786952540689588817750170447535098358471659156491398281308780468750000000000000*x^21 + 461960572757112428337287616885158330788315131952652104679600631086182534570312500000000000000*x^23 - 1146092219988298295525139888083656164725380511791770042322093445414333203125000000000000000000*x^25 + 2414720913348708068413664634285129768118952512297861059683197312940722656250000000000000000000*x^27 - 3193198194662295556062337205217679770829807793229600484548989872167968750000000000000000000000*x^29 + 1883757413126086328093409753134550228006010459329658008583933984375000000000000000000000000000*x^31 - 322712071938389977399719935881920764658291882002981217834375000000000000000000000000000000000*x^33 - 3230767935256908653237206701634808249078297818345502343750000000000000000000000000000000000*x^35 - 38169328359380776176552792189919853875111887728294531250000000000000000000000000000000000000*x^37 + 12428796466552002074684776877135387143837737621093750000000000000000000000000000000000000000*x^39 - 64049784106246065821694925878412363390406250000000000000000000000000000000000000000000000*x^41 + 76497377905218657200515947652460625000000000000000000000000000000000000000000000000000*x^43)*(dx)^6 + (-23519433397776242227180111846835587384762443900654455601600 + 16449267332228438085307914089397386423131542950081509597315057520*x^2 + 465470723308822039225592937600470820262977374219059225864945532076000* x^4 + 10352595599989366514953565364501935970406942543770515465920560141993930000*x^6 - 512935420417531983047233392805810641026425653408895072099103962496799702400000*x^8 + 1387514068147058093168090738549937571657443393829074199308301118614178703962500000*x^10 - 360396024944620468135240657575323072957271558144565175316143976628305747041775000000*x^12 + 134476987743484730706330453153227873929494476396493565288250833269211921089082500000000*x^14 - 2044515761872732777170900882803477760225716005029275800185768478534997115845312500000000*x^16 + 11461122988963924047941904000981249486804100396629979717204213989331200633484375000000000*x^18 - 56491108892843294295003017226154798413467057321241319411231637503178504820312500000000000*x^20 + 213591573415738539230203668122485443270571341613545139199948736148498664062500000000000000*x^22 - 551619100553301499088815446213159356491905374621543744139837106606558593750000000000000000*x^24 + 1248768110588761821100004111443944194648933380586609463950725790687500000000000000000000000*x^26 - 1967724079740172383678388113808950956561672003015838095482364148437500000000000000000000000*x^28 + 1315686257860497442368116262327660714851777709854367238238937500000000000000000000000000000*x^30 - 242687128001056606873134775430729439621925064373224745187500000000000000000000000000000000*x^32 - 9681864144008925978482887230450352887509280310322531250000000000000000000000000000000000*x^34 - 25362677635525527784323900591372485586565227453562500000000000000000000000000000000000000*x^36 + 9165966489395748489918013623825339761521590937500000000000000000000000000000000000000000*x^38 - 47293216462572638554665404590346797882500000000000000000000000000000000000000000000000*x^40 + 56505193314429853302261102428850000000000000000000000000000000000000000000000000000*x^42)*(dx)^7 + (23519433397776242227180111846835587384762443900654455601600*x + 8328070590674632349697624158916893436228036788188819128219354160*x^3 + 466276415418979629711224336565200497354279737357705097524336419840800* x^5 - 5829875633296841563106719982033797422516904229160863894333580326719610000*x^7 + 205843748484541402889700622259322789313118319572003779747366874674541488850000*x^9 - 38405914001742020959791005333534560601276202361209071340989995423261447034000000*x^11 + 22410523853862492298225766630912203509882425462165104476238325754824928786912500000*x^13 - 464119235517980379534853243263762800119206539049038291868885396240128191903750000000*x^15 + 3431775934349111756450014279003210539513902817058513831033918657683065083593750000000*x^17 - 15462420175582177796083296378288266005718753576123963682567943843909274687500000000000*x^19 + 69819051994673292315440129628171900329098414802550384278449798319846699218750000000000*x^21 - 187223834587086904282987693139471025951778844946232997670486239937304687500000000000000*x^23 + 433516580653624792406210314068737264053449923672923182017513798095703125000000000000000*x^25 - 808873947270326495847340931044435014874907420711665090427003466796875000000000000000000*x^27 + 600410239343004300816665880021321965725389942888263715296835937500000000000000000000000*x^29 - 117296808349650946318481025211737589400434321605276341250000000000000000000000000000000*x^31 - 6920698687567303887606654836799134943715934878540546875000000000000000000000000000000*x^33 - 11118088639795889850421489814325079115451299531953125000000000000000000000000000000000*x^35 + 4357365699395006308086501906144217057616012109375000000000000000000000000000000000000*x^37 - 22502025213696999718213725351819267440625000000000000000000000000000000000000000000*x^39 + 26892109938940920060264390218062500000000000000000000000000000000000000000000000*x^41)*(dx)^8 + (-5879858349444060556795027961708896846190610975163613900400*x^2 - 147509492303777696882288558671519019277228186424497556155939416*x^4 - 2483935618183967450748910839724295972729773227579807605235339871545*x^6 + 194857067264406004057733124298000569575431378369225216852346414078232750* x^8 - 547026623271500074079634960981062111824805923620733863342914866080717356250*x^10 + 132477651430440494809288185540531691953342470442987508986053542060188616500000*x^12 - 53243239523068095558883958882546122368418541206856306047661000643914887815625000*x^14 + 655114522892465660281696092864738726147450715516384149089043045822760311562500000*x^16 - 2595381480964638476313152534957628351919441675289979002626757627262479765625000000*x^18 + 14666483838463096201955803437961419368060944262183107243783544163812265625000000000*x^20 - 40826708439184361791021485677174199684913065708916556021620527691679687500000000000*x^22 + 96417342493144786006418891440456621704515190803086958537894719453125000000000000000*x^24 - 205840745983996648711402806230849387609601903868214248294176601562500000000000000000*x^26 + 163223156708833901816633142389584893537039156712049848068906250000000000000000000000*x^28 - 32925491771810776597904220700137501079249068756975300937500000000000000000000000000*x^30 - 2277138327201511961818790698633798067777919801165000000000000000000000000000000000*x^32 - 2949947049949785386216353134196457039710815566250000000000000000000000000000000000*x^34 + 1212172081367195882331189501138835240398876562500000000000000000000000000000000000*x^36 - 6262795835684900871824871295561085287500000000000000000000000000000000000000000*x^38 + 7485724168428026094624173301750000000000000000000000000000000000000000000000*x^40)*(dx)^9 + (-1959952783148020185598342653902965615396870325054537966800*x^3 - 131697660205884142319893408527546810282955446772548697789459416*x^5 + 1421470589043716328624989696609163894979587096634879235021440653945*x^7 - 61129805456770043083450862039457913617027088553938988755460181048585250* x^9 + 10772926122039332283409944498464588725369785275104858740542250326627341250*x^11 - 6509649718475078603185617098176868699869579256793680990522184710776127750000*x^13 + 80769923227867180715396909596286584392208484711512859682289595101540740625000*x^15 - 318795610517795804321943599012219235385757324686545618727617419368395312500000*x^17 + 1808618526655030898673134053063034114448842351556738675779904160441015625000000*x^19 - 5035609253962181806540360185143653476429159051868881120748572520312500000000000*x^21 + 11891426860465878052733728860975494948413732505696292246872748242187500000000000*x^23 - 25462316967390438488828005102204688693797679058041464814046289062500000000000000*x^25 + 20220658764527978743018342256318382896299274723717902013281250000000000000000000*x^27 - 4082151848593565418007005807914527127652772322983462500000000000000000000000000*x^29 - 283484245312626240222111286965349928685719580314062500000000000000000000000000*x^31 - 365120044308370451775345294640854550874781810937500000000000000000000000000000*x^33 + 150251587733412838208101400910047959283882812500000000000000000000000000000000*x^35 - 776299927665680672948213724712338187500000000000000000000000000000000000000*x^37 + 927891583148259065565120498750000000000000000000000000000000000000000000* x^39)*(dx)^10$ end$