dsolve(diff(y(x),x) = 1/2*I*x*(I-2*(-x^2+4*ln(a)+4*ln(y(x)))^(1/2))*y(x),y(x))