Double square

Harmonic

Linear

Square

Even

Scaling

Odd

POTENTIAL

WAVE FUNCTION

Fine = 0 (slider from -0.5 to 0.5, step 0.01)
Coarse = 5 (slider from 0 to 20, step 0.1)
f = false
E: Coarse + Fine
Axis(X) = E
Vscale = 1 (slider from 0.1 to 10, step 0.1)
Vscale_1: Vscale
SHO(x) = Vscale_1 * x^(2)
Linear(x) = Vscale_1 * abs(x)
c = true
d = false
e = false
Width = 1.7 (slider from 0.2 to 5, step 0.1)
Square = If(abs(x) < Width / 2, 0, Vscale)
Separation = 1.1 (slider from 0 to 10, step 0.1)
Double = If(abs(x) < Separation / 2, Vscale, If[abs(x) < Width + Separation / 2, 0, Vscale])
V = If(c, SHO, If[d, Linear, If[e, Square, Double]])
A = Intersect(V, Axis, (-2.2360679775002814, 5.000000000002199))
a = true
Even_1: a
Vindex = 2
CTP = If(e, Width / 2, If[f, Width + Separation / 2, x(A)])
Symm = If(Even_1 ≟ true, 1, -1)
Psiscale_1 = 1
Scale = 1 (slider from 0.1 to 10, step 0.1)
b = false
Odd_1: b
Range = 2.5 (slider from 2 to 5, step 0.5)
dx = 0.02
loc3 = SolveODE(0, E - V(x), E^(2) - E * V(x), 0, E, Scale, Range * CTP, 0.02)
loc1 = SolveODE(0, E - V(x), E^(2) - E * V(x), 0, E, Scale, -Range * CTP, -dx)
loc8 = SolveODE(0, E - V(x), E^(2) - V(x) * E, 0, E + Scale, 0, Range * CTP, dx)
loc6 = SolveODE(0, E - V(x), E^(2) - V(x) * E, 0, E + Scale, 0, -Range * CTP, -dx)
text3: "ENERGY = " + E + ""