Написать две программы для интерполяции функции по формулам Лагранжа при равномерном расположении узлов

clear, clc
% координаты узлов
x = input ('введите значение икс в квадратных скобках ')  %[1.2  1.7 2  2.6 3.1 3.3 3.5 3.9];
y = input('введите значение игрек в квадратных скобках '); %[8.278  8.811  9.356 9.621 10.165 10.431 10.701];
x0= input('введите точки');
P_n = zeros(size(x0)); % создаем нулевой массив значений интерполяционного полинома
for t=1:length(x) % в цикле считаем сумму по узлам
    P = ones(length(x),length(x0));
    for i=1:length(x)
        if t~=i
            P(t,:) = P(t,:).*(x0-x(i))./(x(t)-x(i));
        end
    end
    P_n = P_n + y(t)*P(t,:);
   % output L_n;
end
plot(x0,P_n)
hold on
plot(x,y)

clear, clc
% координаты узлов
x = input ('введите значение икс в квадратных скобках ')  %[1.2  1.7 2  2.6 3.1 3.3 3.5 3.9];
y = input('введите значение игрек в квадратных скобках '); %[8.278  8.811  9.356 9.621 10.165 10.431 10.701];
x0= input('введите точки');
y0 = spline(x,y, x0);
% или y0 = interp1(x,y, x0);
plot(x0,y0)
hold on
plot(x,y)

function yy=lagrange2(x,y,xx,n)
  y0=polyfit(x, y, n-1);
  yy=polyval(y0, xx);
end

function yy=lagrange1(x,y,xx)
N=length(x);% узнаем число узлов интерполяции (N=n+1)
 
yy = zeros(size(xx));% создаем нулевой массив значений интерполяционного полинома
% в цикле считаем сумму по узлам
for k=1:N
    % вычисляем произведения
    t=ones(size(xx));
    for j=[1:k-1, k+1:N]
        t=t.*(xx-x(j))/(x(k)-x(j));
    end
    % накапливаем сумму
    yy = yy + y(k)*t;
   
end