15 lines
347 B
Matlab
15 lines
347 B
Matlab
function Lnf = lagrange_barycentric(nodes,f_nodes,x)
|
|
n = length(nodes);
|
|
m = length(x);
|
|
Lnf = zeros(size(x));
|
|
l = zeros(size(nodes));
|
|
for i = 1:m
|
|
ux = prod( x(i) - nodes([1:n]));
|
|
for j = 1:n
|
|
wi(j) = 1/prod( nodes(j) - nodes([1:j-1,j+1:n]));
|
|
endfor
|
|
Lnf(i) = ux* ((wi./(x(i) - nodes([1:n])))*f_nodes');
|
|
|
|
endfor
|
|
end;
|