added some code, but Moperator and Boperator do not have the correct eigenvalues!