def moffat(x, amplitude=1, center=0., sigma=1, beta=1.):
    """ 1 dimensional moffat function:

    moffat(amplitude, center, sigma, beta) = amplitude / (((x - center)/sigma)**2 + 1)**beta
    """
    return amplitude / (((x - center)/sigma)**2 + 1)**beta