      t0 = (sqrt(x**2+y**2+z**2)**3+2*M*x**2)/sqrt(x**2+y**2+z**2)**3
                               Cxyz_g(1,1)=t0

      t0 = 2*M*x*y/sqrt(x**2+y**2+z**2)**3
                               Cxyz_g(1,2)=t0

      t0 = 2*M*x*z/sqrt(x**2+y**2+z**2)**3
                               Cxyz_g(1,3)=t0

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

      t0 = (sqrt(x**2+y**2+z**2)**3+2*M*y**2)/sqrt(x**2+y**2+z**2)**3
                               Cxyz_g(2,2)=t0

      t0 = 2*M*y*z/sqrt(x**2+y**2+z**2)**3
                               Cxyz_g(2,3)=t0

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

      t0 = (sqrt(x**2+y**2+z**2)**3+2*M*z**2)/sqrt(x**2+y**2+z**2)**3
                               Cxyz_g(3,3)=t0

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

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

