x0 = 0.0
order = 10
td = TaylorScalar{typeof(x0), order}(x0, one(x0));
td*td # TaylorScalar{Float64, 4}((0.0, 0.0, 2.0, 0.0))
td^2 # TaylorScalar{Float64, 4}((0.0, NaN, NaN, NaN))
Divide by 0 when you have a 0 constant term:
|
for j in 1:(i - 1)]...)) / v[1] |