Given F(x)=\frac{1+4x^2}{1+x^4},

Plot the error of the approximation on the interval [-1, 1] using the data points x_0=-1,~ x_1=0, ~x_2=1.

Now, the formula for the error is given as:

f(x)-p_{2m+1}=(x-x_0)^2(x-x_1)^2\cdots(x-x_m)^2\frac{f^{(2m+2)}(\xi)}{(2m+2)!}

where \xi is in an interval containing x and x_i (i = 0, 1, 2, ..., m).

Ok, so, even though I've found the actual approximating polynomial, p_4 (not 5 since the highest power has zero coefficient), I don't actually need it in this formula, from what I see.

So I can compute the sixth derivative of F(x) in Maple, which is very long to post. And (2m+2)! = 6! = 720. (m=2).

But what do I do about this \xi thing?! I.e. How do I find the correct expression to plot?