Polinomial Lagrange Di Matlab

Berikut ini saya berikan contoh curve-fitting dengan menggunakan polinomial Lagrange dalam script MATLAB. Semoga bisa bermanfaat bagi pembaca semua yang memerlukan

function lagrangePolinomial()
clc; 
X =  [7.083 7.167 7.25 7.5] ;  
Y1 = [0.11 0.003 0.009 0.005];
xx = 7.083:.005:7.5; 
m = length(xx); 
for i=1:m 
    k(i) = lagrange(X,Y1, xx(i)); 
end
plot(X,Y1,'or',  'markersize', 6, 'MarkerEdgeColor','k','MarkerFaceColor','r'); 
hold on; 
plot(xx,k, 'b.'); 
hold on; 


function sum = lagrange(A,B, p)
a = length(A); 
sum = 0; 
for i=1:a
    prod = B(i); 
    for j=1:a
        if i~=j
            prod = prod * ( p - A(j))/ (A(i) - A(j) );  
        end
    end
    sum  = sum + prod; 
end

Hasil eksekusinya bisa dilihat pada gambar berikut: