N=input('Enter the number of N:');pointx=ones(1,N+1);pointy=ones(1,N+1);%*******************for i=1:N+1    pointx(i)=cos(((2.*i-1)./(2.*N+2)).*pi);    pointy(i)=(1+25.*(pointx(i).^2)).^(-1);end%**********************x=linspace(-1,1,1000);y=(1+25.*(x.^2)).^(-1);xi=linspace(-1,1,1000);yi=lagrange_interpolating(xi,pointx,pointy);plot(x,y,xi,yi,'--');