def METHOD_NAME(x, mean=0.0, sigma=1.0):
    """
    Utility function to integrate Gaussian density within pixels
    Args:
        x (1D array): pixel centers
        mean (float): mean of Gaussian
        sigma (float): sigma of Gaussian
    Returns:
        array of integals of the Gaussian density in the pixels.
    Note:
        All pixels must be the same size
    """
    x = (np.asarray(x, dtype=float) - mean) / (sigma*np.sqrt(2))
    dx = x[1]-x[0]
    if not np.allclose(np.diff(x), dx):
        raise ValueError('all pixels must have the same size')
    edges = np.concatenate([x-dx/2, x[-1:]+dx/2])
    assert len(edges) == len(x)+1
    y = scipy.special.erf(edges)
    return (y[1:] - y[:-1])/2