Derivative of Bessel function at x=0.
14 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Selçuk Sehitoglu
el 15 de Mayo de 2022
Comentada: Star Strider
el 21 de Mayo de 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 comentarios
Respuesta aceptada
Star Strider
el 15 de Mayo de 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 comentarios
Más respuestas (0)
Ver también
Categorías
Más información sobre Bessel functions en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!