>> orbit3b([0 200],[2 .20001 2 -.2 0 0 0 0 -2 -.2 -2 .2], [.03 .30 .03], 100000,500,'case2.gif')