From 138f77fe71e43adc24ecb030d73df81cc5c4eb7e Mon Sep 17 00:00:00 2001 From: Norbert Preining Date: Mon, 4 Sep 2017 09:45:51 +0900 Subject: add ketpic --- .../ketlib/ketpicscifiles6/Intersectcrvsf.sci | 47 ++++++++++++++++++++++ 1 file changed, 47 insertions(+) create mode 100755 texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectcrvsf.sci (limited to 'texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectcrvsf.sci') diff --git a/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectcrvsf.sci b/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectcrvsf.sci new file mode 100755 index 00000000..2a2e5b6c --- /dev/null +++ b/texmf-dist/scripts/ketpic/ketlib/ketpicscifiles6/Intersectcrvsf.sci @@ -0,0 +1,47 @@ +// 08.10.26 +// 10.02.16 + +function Out=Intersectcrvsf(varargin) + Nargs=length(varargin); + Eps0=10^(-4); + Crv=varargin(1); + Fd=varargin(2); + N=3;Uveq=1; + if Nargs>=3 + Tmp=varargin(3); + if type(Tmp)==10 + Uveq=Tmp; + N=4; + end; + end; + Np=[50,50]; + if Nargs>=N + Np=varargin(N); + if length(Np)==1 + Np=[Np,Np]; + end; + end; + Eps=[0.05,1]; // + if Nargs>=N+1 + Eps=varargin(N+1); + end; + if length(Eps)==1 // + Eps=[Eps,1]; // + end; // + Out=[]; + for I=1:Numptcrv(Crv)-1 + PtA=Ptcrv(I,Crv); + PtB=Ptcrv(I+1,Crv); + if norm(PtB-PtA)