ええと･･･恥ずかしながらラグランジュ補間なんて初めて聞きました(笑)。
んで、ネットで検索してみたんですが、
原則的に『コンピューターで処理して』曲線を求める手法みたいですね。
つまり散布データから曲線を類推する手法の一種の様です。
しかしながら原理はそんなに難しくも無いですね。
(もっとも散布データが増えれば大変な計算を要するようですが。)
一番分かり易い説明をしてるのがザーっと見たカンジですと、次のサイトです。

http://www.st.chukyo-u.ac.j\(\frac{p}{h}\)atan\(\frac{o}{l}\)eastsqrshor\(\frac{t}{s}\)annten\(\frac{n}{l}\)aglange.html

ここで行われてる例題がちょうどこのちこうさんの質問にフィットしてると思います。

ちなみに解答は

･f(x)＝\(x^{2}\)･･････☆

となります。

見よう見まねの解き方:

｛(1,1),(2,4),(3,9)｝のサンプルを用いたラグランジュ補間法:

この場合はデータは(\(x_{0}\),\(y_{0}\))＝(1,1),(\(x_{1}\),\(y_{1}\))＝(2,4),(\(x_{2}\),\(y_{2}\))＝(3,9)
の三つなので、ラグランジュの補間で求める関数f(x)は次の様に書き表す事が出来る。

･f(x)＝(x－\(x_{1}\))*(x－\(x_{2}\))*\(y_{0}\)/{(\(x_{0}\)－\(x_{1}\))*(\(x_{0}\)－\(x_{2}\))}
＋(x－\(x_{0}\))*(x－\(x_{2}\))*\(y_{1}\)/{(\(x_{1}\)－\(x_{0}\))*(\(x_{1}\)－\(x_{2}\))}
＋(x－\(x_{0}\))*(x－\(x_{1}\))*\(y_{2}\)/{(\(x_{2}\)－\(x_{0}\))*(\(x_{2}\)－\(x_{1}\))}
＝(x－2)*(x－3)*1/{(1－2)*(1－3)}＋(x－1)*(x－3)*4/{(2－1)*(2－3)}
＋(x－1)*(x－2)*9/{(3－1)*(3－2)}

これを整理すれば☆が求められます。
以上です。