use of polyxpoly to find roots

I have used this code to find the roots but it is giving an error.

clc

clf

b=linspace(0,1);

x=8;

p= x.*(sqrt(1-b)).*(((besselj(1,x.*sqrt(1-b)))./(besselj(0,x.*sqrt(1-b)))));

q = x.*(sqrt(b)).*((besselk(1,x.*sqrt(b)))./(besselk(0,x.*sqrt(b))));

plot(b,p)

hold on

plot(b,q)

[xx,yy] = polyxpoly(b,p,b,q)

Answer 1

[Pravin Jagtap]

Hello Vipul,

Refer to the following documentation for understanding how 'ployxpoly' works:

https://www.mathworks.com/help/map/ref/polyxpoly.html

From your code, I can see that the first value of 'q' vector is 'NaN' (Not a Number) which is not acceptable. Please check how you can handle that corner case.

Hope this will help you.