import numpy as np

import matplotlib.pyplot as plt

h0 = 2.27

R = 4

x2 = y2 = np.linspace(-10,10,110)

x2v, y2v = np.meshgrid(x2,y2,indexing ='ij',sparse = False)

h2v = h0/(1 + (x2v**2 + y2v**2)/(R**2))

dhdx, dhdy = np.gradient(h2v)

plt.plot(x2,dhdx)