Differential equations with Maxima in SageMath in Jupyter in CoCalc

Author: Samuel Lelièvre
Date: 2020-03-25
License: CC-BY

Prompted by Ingo Dahn's question on the CoCalc mainling list:

https://groups.google.com/d/topic/cocalc/ThsLbZ59UvQ/discussion

In [1]:

a = maxima('loesung: rk(t - x^2, x, 1, [t, 0, 8, 0.1])')

In [2]:

[[0.0,1.0],[0.1,0.9137945035852622],[0.2,0.8511915444670534],[0.3,0.8076220416607769],[0.4,0.7798078343548364],[0.5,0.7652812115993086],[0.6000000000000001,0.7620928666387918],[0.7000000000000001,0.7686296203599586],[0.8,0.7834997571928738],[0.9,0.8054620963406812],[1.0,0.833384386129039],[1.1,0.86622162622296],[1.2,0.9030076621556891],[1.3,0.942855002832351],[1.4,0.9849589035341213],[1.5,1.028602655744726],[1.6,1.073161868443532],[1.7,1.118106337679893],[1.8,1.162998846308774],[1.9,1.207490859391788],[2.0,1.251315539722856],[2.1,1.294278785459053],[2.2,1.336249099869959],[2.3,1.377147075451045],[2.4,1.416935155456955],[2.5,1.455608169722248],[2.6,1.493184965398314],[2.7,1.529701292830062],[2.8,1.565203976807658],[2.9,1.599746309038329],[3.0,1.633384537470931],[3.1,1.666175296833841],[3.2,1.698173815612652],[3.3,1.729432740882358],[3.4,1.760001437982126],[3.5,1.789925642304432],[3.6,1.819247362095181],[3.7,1.848004951916355],[3.8,1.876233295024788],[3.9,1.903964048771379],[4.0,1.931225920092155],[4.100000000000001,1.958044948404142],[4.2,1.984444781050193],[4.3,2.010446932238823],[4.4,2.036071020584224],[4.5,2.061334983225331],[4.600000000000001,2.086255266401461],[4.7,2.110846993543088],[4.800000000000001,2.135124112605639],[4.9,2.159099524690948],[5.0,2.182785196084777],[5.100000000000001,2.206192255776835],[5.2,2.229331080383641],[5.300000000000001,2.252211368206678],[5.4,2.274842203956185],[5.5,2.297232115472118],[5.600000000000001,2.319389123588096],[5.7,2.341320786116751],[5.800000000000001,2.363034236787538],[5.9,2.384536219840627],[6.0,2.405833120871673],[6.100000000000001,2.426930994430188],[6.2,2.447835588796759],[6.300000000000001,2.46855236829958],[6.4,2.489086533476515],[6.5,2.509443039343724],[6.600000000000001,2.529626611994107],[6.7,2.549641763717191],[6.800000000000001,2.56949280680564],[6.9,2.589183866191267],[7.0,2.608718891034662],[7.100000000000001,2.628101665376626],[7.2,2.647335817946111],[7.300000000000001,2.666424831207797],[7.4,2.685372049722583],[7.5,2.704180687885738],[7.600000000000001,2.722853837100146],[7.7,2.741394472435695],[7.800000000000001,2.759805458820331],[7.9,2.778089556803466],[8.0,2.796249427928215]]

In [3]:

point2d(a)