>> orbit3b([0 100],[2 .2 2 -.2 0 0 0 0 -2 -.2 -2 .2], [.03 .30 .03], 50000,500,'case1.gif')