interpolazione polinomiale di lagrange in C

[fatoldsun]

Ciao a tutti,

avrei la necessità di implementare in C un algoritmo che, data una funzione da interpolare, un intervallo su cui farlo ed il numero di nodi sui quali si intende interpolare, costruisca il polinomio interpolatore.
Credo che il metodo migliore sia quello della rappresentazione del polinomio interpolante mediante la formula di Newton, quella che si basa sulle differenze divise.
Ho provato ad "inventare" qualcosa ma il programmino non funziona.
Sapreste aiutarmi?

vi ringrazio in anticipo,

ciao