Parametrizzazione con MatLab

[Vegastar]

Ciao a tutti!!!
Ho veramente tanto bisogno di aiuto. Sto lavorando sull'argomento parametrizzazione in MatLab. Ho alcuni esempi e devo provare che sono giusti. Ma proprio non ci riesco.
Allora, ho dei punti e delle derivate che vengono dal cerchio unitario $ f(s)=(cos(s),sin(s)) $ nei valori $ s_i $ con $ i=0,...,n $ . Usando la parametrizzazione rispetto alla lunghezza della corda (e anche altre, ma mi basta capire questo caso) devo trovare il polinomio $ p_n $ e poi calcolarne la distanza da $ f $ prendendo il massimo della distanza da $ p_n $ a 300 punti equispaziati del cerchio...
Devo usare la cubica di Lagrange sui valori $ (s_0,...,s_3)=(0.0,0.2,1,1.5) $ e la cubica di Hermite sui valori $ (s_0,...,s_3)=(0,0,2,2) $ .
Per Lagrange ho scaricato questa funzione: http://dmpeli.math.mcmaster.ca/matlab/M ... ngeInter.m.
Per Hermite uso pchip. Inoltre, per computare la parametrizzazione rispetto alla lunghezza della corda ho costruito questa funzione

function u=lunghcorda(p)
%
% Funzione che genera il vettore dei parametri dai punti della curva
% rispetto alla lunghezza della corda.
%
% INPUT:
% p: matrice di dimensioni 2xn contenente le coordinate dei punti della 
% curva.
% 
% OUTPUT:
% u: vettore di lunghezza n contente i valori dei parametri per la curva.
%
n=size(p,2);
u=[];
u(1)=0;
for i=2:n
    k=sqrt((p(1,i)-p(1,i-1))^2+(p(2,i)-p(2,i-1))^2);
    u(i)=u(i-1)+k;
end

In Matlab, inserisco i valori $ s=(s_0,...,s_3) $ e costruisco la matrice $ fs=( ( cos(s) ),( sin(s) ) ) $ . Trovo poi il vettore dei valori del parametro $u=lunghcorda(fs)$. A questo punto, però, mi perdo. Forse è stupido, ma ci sono stata sopra talmente tanto che mi sono confusa e non so più come fare...
Cosa devo fare ora per trovare il polinomio, avendo i valori $ u $ del parametro?
Per favore, datemi una mano, non so più dove sbattere la testa!