dsolve(2*diff(diff(y(x),x),x)*y(x)+diff(y(x),x)^2+1=0,y(x))