def model(self, height, xo, yo, sx, sy, theta, m):
    dx = (self.x - xo)
    dy = (self.y - yo)
    a = (((np.cos(theta) ** 2) / (2 * (sx ** 2))) + ((np.sin(theta) ** 2) / (2 * (sy ** 2))))
    b = (((- np.sin((2 * theta))) / (4 * (sx ** 2))) + (np.sin((2 * theta)) / (4 * (sy ** 2))))
    c = (((np.sin(theta) ** 2) / (2 * (sx ** 2))) + ((np.cos(theta) ** 2) / (2 * (sy ** 2))))
    psf = (height * np.exp((- (((a * (dx ** 2)) + (((2 * b) * dx) * dy)) + (c * (dy ** 2))))))
    return (psf + m)
def model(self, height, xo, yo, sx, sy, theta, m):
    dx = (self.x - xo)
    dy = (self.y - yo)
    a = (((np.cos(theta) ** 2) / (2 * (sx ** 2))) + ((np.sin(theta) ** 2) / (2 * (sy ** 2))))
    b = (((- np.sin((2 * theta))) / (4 * (sx ** 2))) + (np.sin((2 * theta)) / (4 * (sy ** 2))))
    c = (((np.sin(theta) ** 2) / (2 * (sx ** 2))) + ((np.cos(theta) ** 2) / (2 * (sy ** 2))))
    psf = (height * np.exp((- (((a * (dx ** 2)) + (((2 * b) * dx) * dy)) + (c * (dy ** 2))))))
    return (psf + m)