I want to obtain y(x) of following answer with Mathematica:
and my code is
Clear[y, λ, x, z, u, v]
eqn = y''[x] == λ *x*y[x]/Sqrt[-1 + x]
u[x_] := Sqrt[-1 + x]
y[x_] := z[u[x]]
eqn2 = Simplify[eqn /. x -> InverseFunction[u][v], v > 0]
zSolution[v_] = z[v] /. First[DSolve[eqn2, z[v], v]]
ySolution[x_] = zSolution[u[x]]

DSolvefrom the start. But as it stands, there's probably nothing better available within Mathematica. Provide all known assumptions. – Jens May 01 '15 at 16:44NDsolveapproach to the other question. If that's all you need then all this variable substitution is a lot of work for no benefit. – george2079 May 01 '15 at 20:05