function out=ntstar(x_p,y_p,s)
out = newton(@(t) arclength(x_p,y_p,t)-s,@(t) sqrt(x_p(t)^2+y_p(t)^2),.5,.5e-3);