The NEMO package (a collection of Unix programs) does not have a python interface (yet), so just a few comments on an integration of the type of orbit we discussed before, but now using NEMO. If you had preloaded NEMO into your shell before you started the notebook, you can also create cells with the magic
%% script bash
and executing the cell would now execute NEMO (unix) commands. Otherwise, simply cut and paste these commands in a terminal where NEMO was preloaded. Typically this would be something like
source /astromake/opt/nemo/4.0.0/nemo_start.sh
after which your bash shell has a number of new commands.
So here are those commands:
mkorbit out=orb1 x=1 vy=0.1
orblist orb1
# the interface is described here
more $NEMO/src/orbit/potential/data/plummer.c
and uses the same expression as we saw before $$ \Phi = - { M \over { {(r_c^2 + r^2)}^{1/2} } } $$
orbint orb1 orb1a 200 0.1 1 potname=plummer mode=euler
...
### nemo Debug Info: Energy conservation: 0.272269
orbint orb1 orb1b 200 0.1 potname=plummer mode=leapfrog
...
orbint orb1 orb1c 200 0.1 1 potname=plummer mode=rk2
...
### nemo Debug Info: Energy conservation: 0.000465693
orbplot orb1a
orbplot orb1b
orbplot orb1c
To create the graphics in the previous commands, we had to cheat, and convert the postscript file to PNG file using one of the ImageMagick commands
orbplot orb1a yapp=orb1a.ps/vps
convert orb1a.ps orb1a.png
ImageMagick is one of those packages you should install!