function out = npartition(x_p,y_p,n,C)
out = zeros(1,n);
if nargin < 4
    C = @(s) s;
end
for i=0:n
    out(i+1)=ntstar(x_p,y_p,C(i/n));
end