---

import math

def gaussian_integral_taylor(a, b, eps=1e-10, max_terms=1000):
    if a == b:
        return 0.0

    total = 0.0
    n = 0
    sign = 1.0
    b_pow = b
    a_pow = a
    fact = 1.0
    denom = 1

    term = sign * (b_pow - a_pow) / (fact * denom)

    while abs(term) > eps and n < max_terms:
        total += term
        n += 1
        sign = -sign
        b_pow *= b * b
        a_pow *= a * a
        fact *= n
        denom = 2 * n + 1
        term = sign * (b_pow - a_pow) / (fact * denom)

    return total

if __name__ == "__main__":
    a = float(input("a = "))
    b = float(input("b = "))
    eps = float(input("eps = "))
    max_terms = int(input("max_terms = "))

    result = gaussian_integral_taylor(a, b, eps, max_terms)
    print(f"Integral = {result:.12g}")