// 08.05.18 Koshikawa // 08.05.19 Changed // 08.00.04 // 09.11.07 for the same curve function KL=IntersectcrvsPp(varargin) // Modified Nargs=length(varargin); G1=varargin(1); G2=varargin(2); Eps=10.0^(-4); if Nargs>2 Eps=varargin(3) end SqEps=10.0^(-10); Eps2=0.1; if Nargs>3 Eps2=varargin(4) end Data1=G1; Data2=G2; if size(Data1,1)==size(Data2,1) Tmp1=Data2(size(Data2,1):-1:1,:); // 09.11.07 Eps0=10^(-6); Tmp2=norm(Data1-Data2); Tmp3=norm(Data1-Tmp1); if Tmp20 Tmp=Op(1,KL1); P=Op(1,Tmp); T=Op(2,Tmp); I=Op(3,Tmp); J=Op(4,Tmp); Tmp=MixS(P,I+T,J); KL=MixL(Tmp); end for N=2:Mixlength(KL1) Tmp=Op(N,KL1); P=Op(1,Tmp); Flg=0; for K=1:Mixlength(KL) Tmp=Op(K,KL); if Vecnagasa2(P-Op(1,Tmp))