Making y = z^2 we obtain
z (x z + 2 z'-1) = 0 or z=0 and x z + 2 z'-1=0
We discard z=0 due to the initial conditions so we follow with
2 z'+x z-1=0 which is a linear non-homogeneous differential equation with solution
z = C e^(-(x^2/4))+ e^(-(x^2/4)) int_0^(x/2)e^(-xi^2) d xi and then
y = pm sqrt( e^(-(x^2/4))(C+ int_0^(x/2)e^(-xi^2) d xi))
and
y(0) = pm sqrt(C) = 1 then C = 1 and
y = pm sqrt( e^(-(x^2/4))(1+ int_0^(x/2)e^(-xi^2) d xi))