>> orbit2b([0 200],[2 0.2 2 -0.2 0 -0.02 0 .02],[0.03 0.3], 10000, 50,'twobodyfig2.gif')