def apply_victorthesis_eps(epsilon, q, s_hat):
    """ Return f(eps, q, s_hat) as defined by victor rule"""
    c = [0, 0.13, 0.09, 0.41, -1.65]
    n = [0, 1, -1, 1]
    return (c[1] * q ** n[1] + c[2] * epsilon ** n[2] + c[3] * s_hat ** n[3] + c[4])