def lorentzian(height, center_x, center_y, width_x, width_y):
    """Returns a gaussian function with the given parameters"""
    width_x = float(width_x) ** 2
    width_y = float(width_y) ** 2
    return lambda x, y: height / (1.0 + 4 * (x - center_x) ** 2 / width_x +
        4 * (y - center_y) ** 2 / width_y)
