def polynom(v_atm_n, f, t, alpha, beta, rho, volvol):
    """Debug function - to remove"""
    f_ = f ** (1 - beta)
    p = [-beta * (2 - beta) / (24 * f_ ** 2) * t * f ** beta, t * f ** beta *
        rho * beta * volvol / (4 * f_), (1 + t * volvol ** 2 * (2 - 3 * rho **
        2) / 24) * f ** beta, -v_atm_n]
    return p[0] * alpha ** 3 + p[1] * alpha ** 2 + p[2] * alpha + p[3]
