Ej9

f=@(x) 1./(1+25*x.^2);
x=linspace(-1,1,500);
y_true=f(x);
figure;
plot(x,y_true,'r','linewidth',1);
hold on;
N=10;
xint=linspace(-1,1,N+1);
yint=zeros(1,n1);
y=f(xint);
n1=length(x);
for i=1:n1
    yint(i)=Lagrange(xint,y,x(i));
end
plot(x,yint,'k')
hold on

p=polyfit(xdata,ydata,N);
y_fit=polyval(p,x);
poly_10=plot(x,y_fit,'b','linewidth',1);
hold on;
figure
plot(x,y_true,'r','linewidth',1);
hold on
N=20;
xdata=linspace(-1,1,N+1)';
ydata=f(xdata);
p=polyfit(xdata,ydata,N);
y_fit=polyval(p,x);
poly_20=plot(x,y_fit,'g','linewidth',1);
hold on;
fprintf('\nEn la figura 1 se representan:\nEn rojo la función original.\nEn negro la función aproximada usando Lagrange (10 puntos).\nEn azul, la función aproximada usando polyfit (10 puntos).\n\nEn la figura 2 se observa en rojo la función original y en verde \nla función aproximada usando polyfit con 20 puntos.\n\n ')

Undefined function or variable 'n1'.

Error in Ej9 (line 9)
yint=zeros(1,n1);