Fix docs for σ^± and extend tests to higher spins

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.