diff options
author | Norbert Preining <norbert@preining.info> | 2018-02-12 14:52:54 +0900 |
---|---|---|
committer | Norbert Preining <norbert@preining.info> | 2018-02-12 14:52:54 +0900 |
commit | 8c3c49fd56de185396400a0e8a636ac03d8239ab (patch) | |
tree | adf96aee3f0bfbdc9df3a9e2dd86332a9dfee9cd /texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/MakeBowdata.sci | |
parent | 5a993b26c76714f0810fdc7e3ae7bdc5eb1afd10 (diff) |
update ketcindy
Diffstat (limited to 'texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/MakeBowdata.sci')
-rw-r--r-- | texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/MakeBowdata.sci | 49 |
1 files changed, 0 insertions, 49 deletions
diff --git a/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/MakeBowdata.sci b/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/MakeBowdata.sci deleted file mode 100644 index 27e70d48..00000000 --- a/texmf-dist/scripts/ketcindy/ketlib/ketpic2escifiles6/MakeBowdata.sci +++ /dev/null @@ -1,49 +0,0 @@ -// -// 08.08.24 -// 10.01.02 - -function Out=MakeBowdata(PA,PB,H) - global BOWSTART BOWEND - BOWSTART=PA; - BOWEND=PB; - Eps=10^(-5); - D=1/2*Vecnagasa(PB-PA); - R=(H^2+D^2)/(2*H); - A1=PA(1); A2=PA(2); - B1=PB(1); B2=PB(2); - if abs(A2-B2)>Eps - x=poly(0,'X'); - y=-(A1-B1)*x/(A2-B2)+(1/2)*(A1^2+A2^2-B1^2-B2^2)/(A2-B2); - Eq1=(A1-x)^2+(A2-y)^2-R^2; - Tmp=coeff(Eq1); - C0=Tmp(1); C1=Tmp(2); C2=Tmp(3); - Det=sqrt(C1^2-4*C0*C2); - Ansx1=(-C1+Det)/(2*C2); - Ansx2=(-C1-Det)/(2*C2); - Tmp=horner(y,Ansx1); - PC=[Ansx1,Tmp]; - Tmp=horner(y,Ansx2); - PC2=[Ansx2,Tmp]; - else - Tmp=0.5*(PA+PB); - PC=Tmp+[0,R-H]; - PC2=Tmp-[0,R-H]; - end; - VA=PA-PC; - VB=PB-PC; - if VA(1)*VB(2)-VA(2)*VB(1)<0 - PC=PC2; - end - AngA=acos((PA(1)-PC(1))/R); - if PA(2)-PC(2)<0 - AngA=-AngA; - end - AngB=acos((PB(1)-PC(1))/R); - if PB(2)-PC(2)<0 - AngB=-AngB; - end - if AngA>AngB - AngB=AngB+2*%pi; - end - Out=Mixmake(PC,R,AngA,AngB); -endfunction |