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/ProjCurve.sci | |
parent | 37bfdbee2c91bab6f0abd54503a0b9c669c9109e (diff) |
add ketpic
Diffstat (limited to 'texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/ProjCurve.sci')
-rwxr-xr-x | texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/ProjCurve.sci | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/ProjCurve.sci b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/ProjCurve.sci new file mode 100755 index 00000000..bbe28c12 --- /dev/null +++ b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/ProjCurve.sci @@ -0,0 +1,32 @@ +// 08.07.11 +// 08.10.07 +// 09.06.02 +// 09.08.23 + +function AnsL=ProjCurve(Curve) + global PHI THETA; + Eps=10^(-6); + AnsL=[]; // Added + SP=sin(PHI); CP=cos(PHI); + ST=sin(THETA); CT=cos(THETA); + for I=1:size(Curve,1) + P=Curve(I,:); + x=P(1); y=P(2); z=P(3); + if x~=%inf + Xz=-x*SP+y*CP; + Yz=-x*CP*CT-y*SP*CT+z*ST; + Tmp=[Xz,Yz]; + if I==1 + AnsL=[Tmp]; + else + Tmp1=AnsL(size(AnsL,1),:); + if Tmp1(1)==%inf | norm(Tmp-Tmp1)>Eps + AnsL=[AnsL;Tmp]; + end; + end; + else + AnsL=[AnsL;%inf,%inf]; + end; + end; +endfunction + |