Interpolacja Lagrange'a

Witam, nie jestem pewien czy umieszczam ten temat w odpowiednim miejscu ale lepszego nie znalazłem. Mam do napisania trzy programiki, eliminacja Gaussa, interpolacja Lagrange'a i metoda Sympleksu. Z gaussem sobie poradziłem, ale problemy zaczeły się przy interpolacji. Wiem o co chodzi w interpolacji ale nie mam koncepcji jak się zabrać ze ten program. Idea jest taka: najpierw wprowadza się z klawiatury n - czyli ilośc węzłów a nastepnie n razy xi oraz fi. Program ma znaleŹć i wypisać wielomian interpolacyjny Lagrange'a w postaci stringu np

n = 4;
x0=-2
y0=-3
x1=-1
y1=3
x2=1
y2=3
x3=2
y3=3

Program powinien wypisać:

w(x) = 1/2x^{3 - x}2 - 1/2x + 4

Proszę o algorytm lub cokolwiek co może być pomocne.

Ja nie rozumiem z czym masz problem? Wiesz na czym polega interpolacja, ale nie znasz algorytmu? To jak to wiesz, ale nie wiesz? ;)

Tutaj masz przykład jak nie rozumiesz na czym to polega:
http://sirjoker.w.interia.pl/mat/numer/_Int_Lagrange.pdf
Zrozumiesz i zaprogramujesz.

Powodzenia!