def model(x, nu, a):
    """
    Model for fitting Figure 3.

    Args:
    x (numpy array): x coordinates
    nu (float): the exponent
    a (float): the constant of proportionality

    Returns:
    yy (numpy array): the y values for the model y = ax^{nu}
    """
    yy = a * x ** nu
    return yy