DSolve[Derivative[1][y][x] - 2*a*x*Derivative[2][y][x] + a^2*Derivative[2][y][x]^2 == 0,y[x],x]