u'=v v'=-omega*omega*u z'=u*u omega'=0 b u b u' b z b z'-1 init omega=1,v=1 @ total=1.0001,dt=.01,xhi=1 done