'Two dimensional Moffat model function'

@staticmethod
def evaluate(x, y, amplitude, x_0, y_0, gamma, alpha):

    rr_gg = ((((x - x_0) ** 2) + ((y - y_0) ** 2)) / (gamma ** 2))
    return (amplitude * ((1 + rr_gg) ** (- alpha)))
