lab3 / lagrangevariant

clc
clear all
%%интреполяция с использованием Лагранжа
n=11;
a=-2;
b=2;
h = (b-a)/(n-1);
x = a:h:b;
x1 = a:0.01:b;
y_r = frunge1(x);
y1_r = frunge1(x1);
N = length(x1);
L=zeros(1,N);

for ii=1:n
p=ones(1,N);
for jj=1:n
if jj~=ii
p=p.*(x1-x(jj))/(x(ii)-x(jj));
end
end
L=L+y_r(ii)*p;
end
figure(1)
pogr=L-y1_r;
plot(x1,y1_r,x1,L,x1,pogr)
hold on
plot(x,y_r,'x')

xx=linspace(-2,2,100);
yy=zeros(size(xx));
%% используя функции матлаб
figure(2)
c = polyfit(x,y_r,n-1);
y_ML = polyval(c,x1);
plot(x1,y_ML, x1,y1_r)