Created by: femtobit
This resolves #845 (closed) and extends the tests to make sure that σ^± = 1/2 (σ^x ± i σ^y) is also true for higher spins.
σ^± = 1/2 (σ^x ± i σ^y)