lab3 / int

x=-2:0.4:2;
y_r=frunge(x);
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(2);
pogr=L-y1_r;
plot(x1, y1_r,x1,L,x1,pogr)
hold on
plot(x,y_r,'x')