Is this helpful?

import sys, getopt

opts, args = getopt.getopt(sys.argv[1:], "g:x:y:z:p:")

for opt, arg in opts:

if opt == '-g': g = arg

if opt == '-x': x = arg

if opt == '-y': y = arg

if opt == '-z': z = arg

if opt == '-p': p = arg

def INH (g, x, y, z, p):

from sympy import diff, symbols

import math

g, x, y, z = symbols('g x y z')

a = float((diff(g,x)))

b = float((diff(g,y)))

c = float(abs(diff(g,z)))

p = float(p)

k = math.sqrt((a**2 + b**2)) / (p + c)

k = math.atan(k)

return k

v = INH(g, x, y, z, p)

print(v)

# sample call: py rg_question.py -g 1.0 -x 1.0 -y 1.0 -z 1.0 -p 1.0

It relies on Python module sympy which can do partial differentiation (https://docs.sympy.org/latest/modules/solvers/pde.html). I am not an expert at mathematics so you might need to amend the above, sympy has a pile of different methods which you could use to suit.

Thank You very Much Derek Colly and May GOD Bless you I will try to use the above .