Вы можете использовать fzero
, чтобы найти нули тогда и только тогда, когда кривая пересекает ось X.Так что
fzero(@(x)(2*x*besselj(1,x)-besselj(0,x)),x0)
даст вам ноль, близкий к x0
.
См. мой ответ для объяснения того, когда и почему fzero
не удаетсяхотя это не относится к функции Бесселя, но это полезно знать) и компромиссы с fsolve
.