[![text{Plot}left[text{Evaluate}left[left{sqrt{y(x)},z(x)right}text{/.}, %right],{x,0,60},text{PlotRange}to text{All},text{PlotPoints}to 200right] left(text{sol}=text{NDSolve}left[left{3 y(x)-0.5 z(x)^2=y(x) (9 eta  y(x)+0.5) z'(x)^2text{/.}, eta to 1,-y'(x)=-eta  y(x) y'(x) z'(x)^2+(4 alpha +18) eta  y(x)^2 z'(x)^2+y(x) z'(x)^2text{/.}, {eta to 1,alpha to 1},y(0)=1,z(0)=0.2right},{z,y},{x,0,60}right]right)][1]][1]