interpola1 
(* interpolazione polinomiale lineare *)

(* valori della funzione polinomiale noti per *)
(* x1<x2<x3<x4<x5:determinare valori funzione *)
(* per xn prossimi ai precedenti valori di x  *)

(* generazione polinomio interpolante         *)
interpola=InterpolatingPolynomial[{1,2,6,24,120},x]

1 + (1 + (3/2 + (11/6 + (53*(-4 + x))/24)*(-3 + x))*(-2 + x))*(-1 + x)
;[o]
          3    11   53 (-4 + x)
1 + (1 + (- + (-- + -----------) (-3 + x)) (-2 + x)) (-1 + x)
          2    6        24
(* cerca valore per valori sucessivi di x *)
interpola/.x->1       (* 1    *)
interpola/.x->1.5     (* -0.25 *)
interpola/.x->2       (* 2    *)
interpola/.x->2.5     (*  4.17 *)
interpola/.x->3       (* 6    *)
interpola/.x->4       (* 24   *)
interpola/.x->5       (* 120  *)
interpola/.x->6	      (*  411  *)
1
;[o]
1

-0.2578125
;[o]
-0.257813

2
;[o]
2

4.1796875
;[o]
4.17969

6
;[o]
6

24
;[o]
24

120
;[o]
120

411
;[o]
411

ritorna a inizio