Тексты файл-функций на языке Math lab
function
otv=R(Y1,Z);
otv=(Y1^2+Z^2)^0.5;
function
otv=R1(L,X,Y1);
otv=(LmX(L,X)^2+Y1^2)^0.5;
function
otv=R2(L,X,Y1);
otv=(LpX(L,X)^2+Y1^2).^0.5;
function
otv=LmX(L,X);
otv=L-X;
function
otv=LpX(L,X);
otv=L+X;
function
otv=R1pZ(L,X,Y1,Z);
otv=(R1(L,X,Y1)^2+Z^2)^0.5;
function
otv=R2pX(L,X,Y1);
otv=(R2(L,X,Y1)^2+X^2)^0.5;
function
otv=R2pZ(L,X,Y1,Z);
otv=(R2(L,X,Y1)^2+(Z^2))^0.5;
function
otv=LpXpR(L,X,Y1,Z);
otv=(LpX(L,X)^2+R(Y1,Z)^2)^0.5;
function
otv=LmXpR(L,X,Y1,Z);
otv=(LmX(L,X)^2+R(Y1,Z)^2)^0.5;
function
otv=Hx(L,X,Y1,Z);
otv=(-1)/(4*3.14)*(Y1/(R1(L,X,Y1)^2)*(1-Z/R1pZ(L,X,Y1,Z))-Y1/(R2(L,X,Y1)^2)*(1-Z/R2pX(L,X,Y1)));
function
otv=Hy(L,X,Y1,Z);
otv=1/(4*3.14)*(Z/(R(Y1,Z)^2)*(LpX(L,X)/LpXpR(L,X,Y1,Z)-LmX(L,X)/LmXpR(L,X,Y1,Z))-(LmX(L,X)/(R1(L,X,Y1)^2)*
(1-Z/R1pZ(L,X,Y1,Z))+LpX(L,X)/(R2(L,X,Y1)^2)*(
1-Z/R2pZ(L,X,Y1,Z))));
function
otv=Hz(L,X,Y1,Z);
otv=Y1/(4*3.14*(R(Y1,Z)^2))*(LpX(L,X)/LpXpR(L,X,Y1,Z)+
LmX(L,X)/LmXpR(L,X,Y1,Z));
20