Derivative of Bessel function at x=0.
18 views (last 30 days)
Show older comments
Selçuk Sehitoglu
on 15 May 2022
Commented: Star Strider
on 21 May 2022
Hello everyone,
I would like to change the expression of derivative of besselj function, because embedded expression is leading some problems.
syms nu x
y=besselj(nu,x);
d_y=diff(y,x);
subs(d_y,x,0); % At this row i get "Division by zero." error, because derivative is defined as -nJ(n,x)/x - J(n-1,x) in MATLAB and it is undefined at x=0.
However, the answer is available with the following expression (mathematically):
-(besselj(n-1,x) - besselj(n+1,x))/2
Is there a way to use this expression?
Thanks in advance,
Ozi
0 Comments
Accepted Answer
Star Strider
on 15 May 2022
syms nu x
y(x) = besselj(nu,x);
d_y = diff(y,x)
d_y0 = limit(d_y,x,0)
d_y0 = d_y(eps) % Cheat With 'eps'
d_y0 = vpa(d_y0)
.
2 Comments
More Answers (0)
See Also
Categories
Find more on Bessel functions in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!