DSolve[a^4*y[x] - 2*a^2*Derivative[2][y][x] - lambda*(-b + a*x)*(-(a^2*y[x]) + Derivative[2][y][x]) + Derivative[4][y][x] == 0,y[x],x]