      t0 = 0
                               C2m1_3(1,1)=t0

      t0 = 0
                               C2m1_3(1,2)=t0

      t0 = 2.D0/15.D0*sqrt(30.D0)/R*sqrt(pi)*M
                               C2m1_3(1,3)=t0

      t0 = 0
                               C2m1_3(2,2)=t0

      t0 = cmplx(0.D0,2.D0/15.D0)*sqrt(pi)*sqrt(30.D0)/R*M
                               C2m1_3(2,3)=t0

      t0 = 0
                               C2m1_3(3,3)=t0

