> 
if theta(t) <> 0 then solution:=dsolve([eq1,eq2,eq3,eq4,CI],numeric,maxfun=0): odeplot(solution, [[t, x(t)]], t = 0 .. 100, thickness = 2); odeplot(solution, [[t, z(t)]], t = 0 .. 100, thickness = 2); odeplot(solution, [[t, theta(t)]], t = 0 .. 100, thickness = 2); odeplot(solution, [[t, alpha(t)]], t = 0 .. 100, thickness = 2); #odeplot(solution,[[t,x(t)],[t,alpha(t)],[t,z(t)],[t,theta(t)]], t=0..100, thickness=2); end ;
