e=EC('37a')
L = e.Lseries

z0 = float(1)

a0 = L(z0)

def L1(eps):
    eps=float(eps)
    return (L(z0+eps) - a0)/eps

a1 = L1(eval('1e-10'))

def L2(eps):
    eps=float(eps)
    return (L(z0+eps) - a0 - a1*eps)/((eps)^float(2))

a2 = L2(eval('1e-7'))


def L3(eps):
    eps=float(eps)
    return (L(z0+eps) - a0 - a1*eps - a2*eps*eps)/((eps)^float(3))

a3 = L3(eval('1e-2'))