def mmf(time, alpha, w0, gamma, rate):
    """
    Computes the Morgan-Mercer-Flodin growth model

    Parameters
    ----------
    time : time
    alpha : upper asymptote
    w0 : the value at time = 0
    gamma : parameter that controls the point of inflection
    rate : growth rate

    See Also
    --------
    mmf_inverse

    References
    ----------
    .. [1] A. Khamiz, Z. Ismail, and A. T. Muhammad, "Nonlinear growth models
           for modeling oil palm yield growth," Journal of Mathematics and
           Statistics, vol. 1, no. 3, p. 225, 2005.
    """
    result = (w0 * gamma + alpha * time ** rate) / (gamma + time ** rate)
    return result
