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.

1 次查看(过去 30 天)
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;

采纳的回答

Walter Roberson
Walter Roberson 2019-1-25
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)

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by