def get_approximated_rate_func(self, n):
    new_v_rate_func = (lambda t: v_rate_func((np.floor((t * n)) / n)))
    (max_integral, err) = scipy.integrate.quad(v_rate_func, 0, 1)

    def result(t):
        (integral, err) = scipy.integrate.quad(new_v_rate_func, 0, t)
        return (integral / max_integral)
    return result
def get_approximated_rate_func(self, n):
    new_v_rate_func = (lambda t: v_rate_func((np.floor((t * n)) / n)))
    (max_integral, err) = scipy.integrate.quad(v_rate_func, 0, 1)

    def result(t):
        (integral, err) = scipy.integrate.quad(new_v_rate_func, 0, t)
        return (integral / max_integral)
    return result