      t0 = -6*M*x**2*y/sqrt(x**2+y**2+z**2)**5
                              Cxyz_dyg(1,1)=t0

      t0 = 2*M*x*(x**2-2*y**2+z**2)/sqrt(x**2+y**2+z**2)**5
                              Cxyz_dyg(1,2)=t0

      t0 = -6*M*y*z/sqrt(x**2+y**2+z**2)**5*x
                              Cxyz_dyg(1,3)=t0

      t0 = 4*M*y/(x**2+y**2+z**2)**2*x
                              Cxyz_dyg(1,4)=t0

      t0 = 2*y*M*(2*x**2-y**2+2*z**2)/sqrt(x**2+y**2+z**2)**5
                              Cxyz_dyg(2,2)=t0

      t0 = 2*M*z*(x**2-2*y**2+z**2)/sqrt(x**2+y**2+z**2)**5
                              Cxyz_dyg(2,3)=t0

      t0 = -2*M*(x**2-y**2+z**2)/(x**2+y**2+z**2)**2
                              Cxyz_dyg(2,4)=t0

      t0 = -6*y*M*z**2/sqrt(x**2+y**2+z**2)**5
                              Cxyz_dyg(3,3)=t0

      t0 = 4*M*z/(x**2+y**2+z**2)**2*y
                              Cxyz_dyg(3,4)=t0

      t0 = -2*M*y/sqrt(x**2+y**2+z**2)**3
                              Cxyz_dyg(4,4)=t0

