Interpolasi Lagrange dengan Matlab

Formula yang digunakan untuk interpolasi Lagrange adalah sebagai berikut:
Sehingga dapat dibuat algoritmanya sebagai berikut:

1. Tetapkan jumlah titik yang diketahui. Untuk menginputkan titik yang diketahui dapat meenggunakan dua array x dan y dengan jumlah data = jumlah titiknya. Dengan dua array akan lebih mudah mengatur perilaku data didalam program. Bisa juga menggunakan banyak array sejumlah titik yang diketahui, sehingga masing-masing pasang data disimpan dalam satu array. Cara ini terlihat lebih sederhana, tetapi lebih sulit dalam mengatur perilaku data. Dalam implementasi ini nantinya akan dipilih cara yang pertama, yakni menggunakan dua array x dan y.
2. Mencari L_i(x) dan P(x). L_i(x) didapat sejumlah titik yang diketahui, sehingga diperlukan perulangan sebanyak titik yang diketahui. Demikian pula P(x) merupakan jumlahan dari perkalian y_i dan L_i(x), sehingga memerlukan perulangan yang jumlahnya sana dengan proses pencarian L_i(x). Untuk mencari L_i(x) diperlukan Q_i(x) dan Q_i(x_i). Karena Q_i(x) merupakan hasil perkalian (x-x_i) sejumlah titik yang diketahui, maka diperlukan perulangan lagi untuk mencarinya. Tetapi yang harus diingat disini adalah bahwa, untuk (x-x_i) tersebut tidak ikut dalam hasil perkalian. Sehingga proses hanya akan dilakukan untuk nilai selain (x-x_i). Untuk Q_i(x_i) dapat dicari setelah Q_i(x) diketahui dengan cara mensubstitusi nilai x_i ke dalam Q_i(x). Setelah Q_i(x) dan Q_i(x_i) diketahui dapat dicari L_i(x). Dan untuk selanjutnya mencari P(x).

Misalnya banyaknya titik yang diketahui adalah b, maka algoritma diatas dapat diperhalus menjadi sebagai berikut:

1. Inputkan b.
2. Dari i = 1 s.d b : Inputkan titik ke i
3. Dari i = 1 s.d b: Cari Q_i(x), Cari Q_i(x_i), Cari L_i(x)
4. Cari P(x)

Implementasinya sebagai berikut:
Output program sebagai berikut: