summaryrefslogtreecommitdiff
path: root/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Evlptablepara.sci
diff options
context:
space:
mode:
Diffstat (limited to 'texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Evlptablepara.sci')
-rwxr-xr-xtexmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Evlptablepara.sci79
1 files changed, 79 insertions, 0 deletions
diff --git a/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Evlptablepara.sci b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Evlptablepara.sci
new file mode 100755
index 00000000..b9c86eef
--- /dev/null
+++ b/texmf-dist/scripts/ketpic/ketlib/ketpic2escifiles6/Evlptablepara.sci
@@ -0,0 +1,79 @@
+// 08.09.10
+// 13.10.21 ( __ added to varibles )
+// 17.05.21 Mdv, Ndv changed to division number
+
+function [Zval__,Xval__,Yval__]=Evlptablepara(MS__)
+ Nargs__=Mixlength(MS__);
+ Eps__=10^(-3);
+ Tmp__=Mixop(1,MS__);
+ FdL__=Fullformfunc(Tmp__);
+ Mdv=50;
+ Ndv__=50;
+ if Nargs__>=2
+ Tmp__=Mixop(2,MS__);
+ if Mixtype(Tmp__)~=1 Tmp__=Mixop(1,Tmp__); end;
+ if length(Tmp__)>1
+ Mdv__=Tmp__(1,1);
+ Ndv__=Tmp__(1,2);
+ else
+ Mdv__=Tmp__
+ if Nargs__==2
+ Ndv__=Mdv__
+ else
+ Tmp1__=Mixop(3,MS__);
+ if type(Tmp1__)==1 & length(Tmp1__)==1
+ Ndv__=Tmp1__;
+ else
+ Ndv__=Mdv__;
+ end
+ end
+ end
+ end;
+ Tmp__=Mixop(5,FdL__);
+ K__=mtlb_findstr(Tmp__,'=');
+ Uname__=part(Tmp__,1:K__-1);
+ Urange__=evstr(part(Tmp__,K__+1:length(Tmp__)));
+ Tmp__=Mixop(6,FdL__);
+ K__=mtlb_findstr(Tmp__,'=');
+ Vname__=part(Tmp__,1:K__-1);
+ Vrange__=evstr(part(Tmp__,K__+1:length(Tmp__)));
+ U1__=Urange__(1); U2__=Urange__(2);
+ V1__=Vrange__(1); V2__=Vrange__(2);
+ Du__=(U2__-U1__)/(Mdv__); //17.05.21
+ Dv__=(V2__-V1__)/(Ndv__); //17.05.21
+ Xyzstr__=[Mixop(2,FdL__),Mixop(3,FdL__),Mixop(4,FdL__)];
+ I__=1;
+ Zval__=[];
+ for v__=V1__:Dv__:V2__
+ v1__=v__-Eps__/2; v2__=v__+Eps__/2;
+ ZuL__=[];
+ for u__=U1__:Du__:U2__;
+ u1__=u__-Eps__/2; u2__=u__+Eps__/2;
+ Tmpv__=strsubst(Xyzstr__,Vname__,'v__');
+ Tmp__=strsubst(Tmpv__,Uname__,'u1__');
+ P1__=evstr(Tmp__);
+ Tmp__=strsubst(Tmpv__,Uname__,'u2__');
+ P2__=evstr(Tmp__);
+ Tmp1__=Parapt(P1__);
+ Tmp2__=Parapt(P2__);
+ Dxu__=(Tmp2__(1)-Tmp1__(1))/Eps__;
+ Dyu__=(Tmp2__(2)-Tmp1__(2))/Eps__;
+ u1__=u__-Eps__/2; u2__=u__+Eps__/2;
+ Tmpu__=strsubst(Xyzstr__,Uname__,'u__');
+ Tmp__=strsubst(Tmpu__,Vname__,'v1__');
+ P1__=evstr(Tmp__);
+ Tmp__=strsubst(Tmpu__,Vname__,'v2__');
+ P2__=evstr(Tmp__);
+ Tmp1__=Parapt(P1__);
+ Tmp2__=Parapt(P2__);
+ Dxv__=(Tmp2__(1)-Tmp1__(1))/Eps__;
+ Dyv__=(Tmp2__(2)-Tmp1__(2))/Eps__;
+ Tmp__=Dxu__*Dyv__-Dxv__*Dyu__;
+ ZuL__=[ZuL__,Tmp__];
+ end;
+ Zval__=[Zval__;ZuL__];
+ end;
+ Yval__=V1__:Dv__:V2__;
+ Xval__=U1__:Du__:U2__;
+endfunction;
+