def METHOD_NAME(data, grid, **kwargs):
    """
    Kernel Density Estimation with Scipy
    Parameters
    ----------
    data : numpy.array
        Data points used to compute a density estimator. It
        has `n x p` dimensions, representing n points and p
        variables.
    grid : numpy.array
        Data points at which the desity will be estimated. It
        has `m x p` dimensions, representing m points and p
        variables.
    Returns
    -------
    out : numpy.array
        Density estimate. Has `m x 1` dimensions
    """
    from scipy.stats import gaussian_kde
    kde = gaussian_kde(data.T, **kwargs)
    return kde.evaluate(grid.T)