import numpy as np


def grav(lat, h):
    g0 = 9.7803253359 * ((1 + 0.00193185265241 * np.sin(np.deg2rad(lat)) ** 2) / (
        np.sqrt(1 - 0.00669437999013 * np.sin(np.deg2rad(lat)) ** 2)))

    res = g0 - (3.0877 * 1e-6 - 4.3 * 1e-9 * np.sin(np.deg2rad(lat)) ** 2) * h + 7.2 * 1e-13 * h ** 2

    return res