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/ketpicscifiles6/Intersectlines.sci | |
parent | 37bfdbee2c91bab6f0abd54503a0b9c669c9109e (diff) |
add ketpic
Diffstat (limited to 'texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectlines.sci')
-rwxr-xr-x | texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectlines.sci | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectlines.sci b/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectlines.sci new file mode 100755 index 00000000..ed92a5aa --- /dev/null +++ b/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectlines.sci @@ -0,0 +1,22 @@ +// 09.10.11 +// 09.10.12 + +function Out=Intersectlines(L1,L2) + Eps0=10^(-6); + P=Ptstart(L1); + a1=P(1); a2=P(2); + V=Ptend(L1)-P; + b1=V(1); b2=V(2); + P=Ptstart(L2); + c1=P(1); c2=P(2); + V=Ptend(L2)-P; + d1=V(1); d2=V(2); + Den=b1*d2-d1*b2; + if abs(Den)<Eps0 + Out=[]; + return; + end; + x=-(d1*b2*a1-c1*b1*d2-d1*a2*b1+d1*b1*c2)/Den; + y=-(-a2*b1*d2+b2*a1*d2-b2*c1*d2+c2*d1*b2)/Den; + Out=[x,y]; +endfunction; |