This script can be pasted into the Eigenmath app.

-- Superposition demo ket0 = (1,0,0,0) ket1 = (0,1,0,0) ket2 = (0,0,1,0) ket3 = (0,0,0,1) "Energy operator" Ehat = 0.38 "eV" outer(ket1,ket1) + 1.50 "eV" outer(ket2,ket2) + 6.55 "eV" outer(ket3,ket3) Ehat "State vector" xi = 1/2 ket0 + 1/2 ket1 + 1/2 ket2 + 1/2 ket3 xi "Expected energy" dot(xi,Ehat,xi)

Run