Afsnit 18.3 opgave 3

L0[x_]=(x-x1)(x-x2)/((x0-x1)(x0-x2));
L1[x_]=(x-x0)(x-x2)/((x1-x0)(x1-x2));
L2[x_]=(x-x0)(x-x1)/((x2-x0)(x2-x1));

p2[x_]=f0 L0[x] + f1 L1[x] + f2 L2[x]

x0=1;x1=1.02;x2=1.04;f0=1;f1=0.9888;f2=0.9784; p2[x]

Simplify[p2[x]]

g3=Plot[Gamma[x]-p2[x],{x,1,1.1}]

Mere praecis regning (uden afrundingsfejl)

f1=Gamma[1.02]

f2=Gamma[1.04]

p2[x]

Simplify[p2[x]]

g4=Plot[Gamma[x]-p2[x],{x,1,1.1}, 
PlotStyle->{RGBColor[1,0,0]}]

Show[g3,g4]

Bemaerk at afvigelserne begynder at vokse kraftigt, naar der ekstrapoleres.