// 08.05.18 // 08.05.20 function Ptn=Naigai(A,Bdy) V=[1,1]; Call=Mixlength(Bdy); KL=KoutenList(A,V,Bdy); Ptn=zeros(1,Call); for K=1:Mixlength(KL) Ten=Op(K,KL); T=Op(1,Ten); NC=Op(4,Ten); if T<0 Tmp=modulo(Ptn(NC)+1,2); Ptn(1,NC)=Tmp; end end endfunction