diff options
author | Norbert Preining <norbert@preining.info> | 2017-09-04 09:45:51 +0900 |
---|---|---|
committer | Norbert Preining <norbert@preining.info> | 2017-09-04 09:45:51 +0900 |
commit | 138f77fe71e43adc24ecb030d73df81cc5c4eb7e (patch) | |
tree | cd2d906c634149bbecd2f4bad3d540c16b23befa /texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Makexybdy.sci | |
parent | 37bfdbee2c91bab6f0abd54503a0b9c669c9109e (diff) |
add ketpic
Diffstat (limited to 'texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Makexybdy.sci')
-rw-r--r-- | texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Makexybdy.sci | 111 |
1 files changed, 111 insertions, 0 deletions
diff --git a/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Makexybdy.sci b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Makexybdy.sci new file mode 100644 index 00000000..fb8aa27b --- /dev/null +++ b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Makexybdy.sci @@ -0,0 +1,111 @@ +// 08.09.13 +// 08.09.16 +// 08.10.11 +// 08.10.27 +// 14.07.07 debug ( Kitahara ) + +function EhL=Makexybdy(Fd,Np) + Eps0=10^(-4); + FdL=Fullformfunc(Fd); + Tmp1=Mixop(2,FdL); + Tmp2=Mixop(3,FdL); + Xystr=[Tmp1,Tmp2]; + Tmp=Mixop(5,FdL); + K=mtlb_findstr(Tmp,'='); + Uname=part(Tmp,1:K-1); + Urg=evstr(part(Tmp,K+1:length(Tmp))); + Tmp=Mixop(6,FdL); + K=mtlb_findstr(Tmp,'='); + Vname=part(Tmp,1:K-1); + Vrg=evstr(part(Tmp,K+1:length(Tmp))); + Umin=Urg(1); Umax=Urg(2); + Vmin=Vrg(1); Vmax=Vrg(2); + DrwS=Mixop(7,FdL); + Cflg=0; + EhL=[]; + Tmp=mtlb_findstr(DrwS,'e'); + if length(Tmp)~=0 + Tmp1=strsubst(Xystr,Uname,'('+string(Umax)+')'); + Tmp2=Mixop(6,FdL); + Tmp3='N='+string(Np(2)); + Tmp=Paramplot(Tmp1,Tmp2,Tmp3); + if norm(Ptstart(Tmp)-Ptend(Tmp))<Eps0 + Cflg=1; + Tmp(size(Tmp,1),:)=Tmp(1,:); + end; + EhL=Mixadd(EhL,Tmp); + end; + Tmp=mtlb_findstr(DrwS,'n'); + if length(Tmp)~=0 + Tmp1=strsubst(Xystr,Vname,'('+string(Vmax)+')'); + Tmp2=Mixop(5,FdL); + Tmp3='N='+string(Np(1)); + Tmp=Paramplot(Tmp1,Tmp2,Tmp3); + if Cflg~=0 // 140707 Kitahara + Tmp1=Mixop(Mixlength(EhL),EhL); + Tmp=Joincrvs(Tmp1,Tmp); + if norm(Ptstart(Tmp)-Ptend(Tmp))<Eps0 + Cflg=1; + Tmp(size(Tmp,1),:)=Tmp(1,:); + end; + Tmp1=Mixsub(1:Mixlength(EhL)-1,EhL); + EhL=Mixadd(Tmp1,Tmp); + else + Cflg=0; + if norm(Ptstart(Tmp)-Ptend(Tmp))<Eps0 + Cflg=1; + Tmp(size(Tmp,1),:)=Tmp(1,:); + end; + EhL=Mixadd(EhL,Tmp); + end; + end; + Tmp=mtlb_findstr(DrwS,'w'); + if length(Tmp)~=0 + Tmp1=strsubst(Xystr,Uname,'('+string(Umin)+')'); + Tmp2=Mixop(6,FdL); + Tmp3='N='+string(Np(2)); + Tmp=Paramplot(Tmp1,Tmp2,Tmp3); + if Cflg~=0 //if Cflg==0 + Tmp1=Mixop(Mixlength(EhL),EhL); + Tmp=Joincrvs(Tmp1,Tmp); + if norm(Ptstart(Tmp)-Ptend(Tmp))<Eps0 + Cflg=1; + Tmp(size(Tmp,1),:)=Tmp(1,:); + end; + Tmp1=Mixsub(1:Mixlength(EhL)-1,EhL); + EhL=Mixadd(Tmp1,Tmp); + else + Cflg=0; + if norm(Ptstart(Tmp)-Ptend(Tmp))<Eps0 + Cflg=1; + Tmp(size(Tmp,1),:)=Tmp(1,:); + end; + EhL=Mixadd(EhL,Tmp); + end; + end; + Tmp=mtlb_findstr(DrwS,'s'); + if length(Tmp)~=0 + Tmp1=strsubst(Xystr,Vname,'('+string(Vmin)+')'); + Tmp2=Mixop(5,FdL); + Tmp3='N='+string(Np(1)); + Tmp=Paramplot(Tmp1,Tmp2,Tmp3); + if Cflg~=0 //if Cflg==0 + Tmp1=Mixop(Mixlength(EhL),EhL); + Tmp=Joincrvs(Tmp1,Tmp); + if norm(Ptstart(Tmp)-Ptend(Tmp))<Eps0 + Cflg=1; + Tmp(size(Tmp,1),:)=Tmp(1,:); + end; + Tmp1=Mixsub(1:Mixlength(EhL)-1,EhL); + EhL=Mixadd(Tmp1,Tmp); + else + Cflg=0; + if norm(Ptstart(Tmp)-Ptend(Tmp))<Eps0 + Cflg=1; + Tmp(size(Tmp,1),:)=Tmp(1,:); + end; + EhL=Mixadd(EhL,Tmp); + end; + end; +endfunction; + |