      t0 = 1+8*M*q**2/pp**2/R
                                 g_N(1,1)=t0

      t0 = 8*M*q*p/R/pp**2
                                 g_N(1,2)=t0

      t0 = -4*M*(-1+q**2+p**2)*q/R/pp**2
                                 g_N(1,3)=t0

      t0 = -4*M*q/R/pp
                                 g_N(1,4)=t0

      t0 = 1+8*M*p**2/pp**2/R
                                 g_N(2,2)=t0

      t0 = -4*M*(-1+q**2+p**2)*p/R/pp**2
                                 g_N(2,3)=t0

      t0 = -4*M*p/R/pp
                                 g_N(2,4)=t0

      t0 = 1+2*M*(-1+q**2+p**2)**2/pp**2/R
                                 g_N(3,3)=t0

      t0 = 2*M/R*(-1+q**2+p**2)/pp
                                 g_N(3,4)=t0

      t0 = -1+2*M/R
                                 g_N(4,4)=t0

