Derivative of Bessel function at x=0.
13 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Selçuk Sehitoglu
le 15 Mai 2022
Commenté : Star Strider
le 21 Mai 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 commentaires
Réponse acceptée
Star Strider
le 15 Mai 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 commentaires
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Bessel functions dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!