# Singlet.m

         rho=singlet(spin_system,spin_a,spin_b,state_type)

      spin_a  - the number of the first spin in the singlet state

      spin_b  - the number of the second spin in the singlet state