OppC0 = NewOperator("Cr",NF,0) OppA0 = NewOperator("An",NF,0) -- Create the number operator on spin-orbital -- 0 by multiplying a creation with an -- annihilation operator OppN0 = OppC0 * OppA0 -- Identity operator Opp1 = OppC0 * OppA0 + OppA0 * OppC0