      t0 = 6.D0/35.D0/R**2*M*sqrt(21.D0)*sqrt(pi)
                               dxC3m1(1,1)=t0

      t0 = cmplx(0.D0,2.D0/35.D0)/R**2*M*sqrt(21.D0)*sqrt(pi)
                               dxC3m1(1,2)=t0

      t0 = 0
                               dxC3m1(1,3)=t0

      t0 = 0
                               dxC3m1(1,4)=t0

      t0 = 2.D0/35.D0/R**2*M*sqrt(21.D0)*sqrt(pi)
                               dxC3m1(2,2)=t0

      t0 = 0
                               dxC3m1(2,3)=t0

      t0 = 0
                               dxC3m1(2,4)=t0

      t0 = -8.D0/35.D0/R**2*M*sqrt(21.D0)*sqrt(pi)
                               dxC3m1(3,3)=t0

      t0 = 0
                               dxC3m1(3,4)=t0

      t0 = 0
                               dxC3m1(4,4)=t0

