i=[0:1:11];

n=12;

x=cos(((2*i+1)*pi)/(2*n));

y=sin(x);

xr=linspace(x(1),x(12),12);

yr=sin(xr);

xpom=linspace(x(1),x(12));

ypom=sin(xpom);

xpom1=linspace(xr(1),xr(12));

ypom1=sin(xpom1);

[C,L]=newpoly(x,y);

[D,L]=newpoly(xr,yr);

xx=linspace(x(1),x(12),1000);

yy=polyval(C,xx);

yyr=polyval(D,xx);

subplot(4,1,1)

hh1=plot(x,y,'rx',xpom,ypom,'b-',xx,yy,'k-');

legend('punkty','funkcja','wielomian')

subplot(4,1,3)

hh3=plot(xr,yr,'ro',xpom1,ypom1,'g-',xx,yyr,'k-');

Emaxcz=abs((1/factorial(12)).*max(sin(xx)).*(xx-x(1)).*(xx-x(2)).*(xx-x(3)).*(xx-x(4)).*(xx-x(5))...

.*(xx-x(6)).*(xx-x(7)).*(xx-x(8)).*(xx-x(9)).*(xx-x(10)).*(xx-x(11)).*(xx-x(12)));

legend('węzły','funkcja','wielomian');

subplot(4,1,2)

hh2=plot(xx,Emaxcz,'k-');

legend('błąd interpolacji czebyszewa')

xxr=linspace(xr(1),xr(12),1000);

Emaxw=abs((1/factorial(12)).*max(sin(xxr)).*(xxr-xr(1)).*(xxr-xr(2)).*(xxr-xr(3))...

.*(xxr-xr(4)).*(xxr-xr(5)).*(xxr-xr(6)).*(xxr-xr(7)).*(xxr-xr(8))...

.*(xxr-xr(9)).*(xxr-xr(10)).*(xxr-xr(11)).*(xxr-xr(12)));

subplot(4,1,4)

hh4=plot(xxr,Emaxw,'k-');

legend('błąd interpolacji dla węzłów roz. równomiernie')

---