@staticmethod
def evaluate(x, y, amplitude, x_0, y_0, gamma, alpha):
    """Two dimensional Moffat model function"""
    rr_gg = ((x - x_0) ** 2 + (y - y_0) ** 2) / gamma ** 2
    return amplitude * (1 + rr_gg) ** -alpha
