def fTaylor(f, x0, x1, steps):
    a = 0
    approx = 0
    for n in range(steps):
        df = f.diff(n)
        approx = approx + (df(a).n()/factorial(n))*((x1-x0)**n)
        n = n + 1
    return approx