N=6;

a=-pi; b=pi;       % INTERVALLO.
h=(b-a)/N;         % PASSO PARAMETRI.

% PUNTI INTP.
t=a:h:b;           % PARAMETRI CURVA INTP.
xt=cos(t);         % ASCISSA PTO IN CURVA.
yt=sin(t);         % ORDINATA PTO IN CURVA.

% TEST PTS.
htest=h/100;       % PASSO PARAMETRO t.
ttest=a:htest:b;   % PARAMETRI CURVA TEST.

% SPLINES INTP. DELLA CURVA PARAMETRICA.
xttest=interp1(t,xt,ttest,'spline');
yttest=interp1(t,yt,ttest,'spline');

% PLOT.
theta=0:htest:2*pi;
plot(xttest,yttest,'k-',cos(theta),sin(theta),'r-');
