ϵ = 7;
lasolucion = 
        NDSolve[{ψ''[x] + (ϵ - (x^2)) ψ[x] == 
              0, ψ[0] == 0, ψ[1] == 
              1}, ψ, {x, -20, 20}];
Plot[{Evaluate[ψ[x] /. lasolucion], (x^2)}, {x, -5, 
        5}]