Kindly check if my code for the attached equation is correct as the result obtained is incorrect. If I find that the code is correct, I'll recheck the values I defined.

The equation is attached below.

CODE:

syms x 
fun=(besselj(x,sqrt(a/b)).*besselj(x,sqrt(b/a)).*besselj(x,(gamma/sqrt(a*b))).*exp(-x.*(d./sqrt(a.*b))))
q=int(fun,x,0,Inf)
M1=pi*u*sqrt(a*b)*q;syms x 
fun=(besselj(x,sqrt(a/b)).*besselj(x,sqrt(b/a)).*besselj(x,(gamma/sqrt(a*b))).*exp(-x.*(d./sqrt(a.*b))))
q=int(fun,x,0,Inf)
M1=pi*u*sqrt(a*b)*q;

Answer 1

[Walter Roberson]

J1(value) in MATLAB becomes besselj(1, value)

J0(value) in MATLAB becomes besselj(0, value)

What you coded is the equivalent of J_x(value)

[Swathi S]

Sir, thank you for your speedy response. I am getting the correct result.