function out = npartition(x_p,y_p,n)
out = zeros(1,n);
s = arclength(x_p,y_p,1);
for i=1:n
    out(i)=ntstar(x_p,y_p,s*i/n);
end