信谊x;
π=符号(π);
情况= 4 *π* 10 ^ 7;
L1 =信谊(88)*符号(10)^ (8);
L2 =信谊(88)*符号(10)^ (8);
=符号(50)*符号(10)^ 3;
b =符号(50)*符号(10)^ 3;
d =符号(20)*符号(10)^ 3;
% p = 0:0.1 * 10 ^ 2:0.5;
p = 0;
eq1 = besselj (1, x * sqrt (a / b));
eq2 = besselj (1, x * sqrt (b / a));
eq3 = besselj (0, x * p /√(a * b));
eq4 = (exp (- x * d /√(a * b)));
eq5 =π*情况* sqrt (a * b);
eqf = eq1 * eq2 * eq3 * eq4
M = eq5 * int (eqf x 0,正)
vpa (M, 16)