// 08.09.13 // 09.03.06 function PL=Partcrv3(T1,T2,Fk) Eps=10^(-4); // Tmp=Mixop(1,Fk); // new part from if T1>T2+Eps Npt=size(Fk,1); Out1=Partcrv3(T1,Npt,Fk); Out2=Partcrv3(1,T2,Fk); Tmp=Fk(1,:)-Fk(Npt,:); if norm(Tmp)Ie+Eps P=(1-T2+Ie)*Fk(Ie,:)+(T2-Ie)*Fk(Ie+1,:); PL=[PL;P] end; endfunction;