VectorPlot[{y,Sqrt[x]}, {x, 0,50}, {y,0,50}]

DSolve[y''[x]==Sqrt[y[x]],y[x],x]

?Hypergeometric2F1

Integrate[Sqrt[y/(c1-y)], y]