Exercício 4.6 (25 dígitos)

> X:=vector(4,[1.10,1.05,1.15,1.00]);
fX:= vector(4,[0.2933,0.2414,0.3492,0.1924]);
poli:=PoliNewton(X,fX,25);
#factor(poli);
PoliNewtonPt(1.09,X,fX,25);

> X:=vector(4,[1.05,1.00,1.10,1.15]);
fX:= vector(4,[0.2414,0.1924,0.2933,0.3492]);
poli:=PoliNewton(X,fX,25);
p:=x->poli:
PoliNewtonPt(1.09,X,fX,25);

> with(plots):
A:=pointplot({[X[1],fX[1]],[X[2],fX[2]],[X[3],fX[3]],[X[4],fX[4]]},axes=BOXED,symbol=circle):
B:=plot(poli(x),x=0.95..1.2):
display({A,B});