diff options
Diffstat (limited to 'texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Partcrv3.sci')
-rw-r--r-- | texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Partcrv3.sci | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Partcrv3.sci b/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Partcrv3.sci new file mode 100644 index 00000000..5b3cd178 --- /dev/null +++ b/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/Partcrv3.sci @@ -0,0 +1,34 @@ +// 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)<Eps + PL=Joincrvs(Out1,Out2); + else + PL=Mix(Out1,Out2); + end; + return; + end; + // to here + Is=ceil(T1); + Ie=floor(T2); + PL=[]; + if T1<Is-Eps + P=(Is-T1)*Fk(Is-1,:)+(1-Is+T1)*Fk(Is,:); + PL=[P] + end; + PL=[PL;Fk(Is:Ie,:)]; + if T2>Ie+Eps + P=(1-T2+Ie)*Fk(Ie,:)+(T2-Ie)*Fk(Ie+1,:); + PL=[PL;P] + end; +endfunction; + |