use of polyxpoly to find roots

1 view (last 30 days)
vipul kumar
vipul kumar on 23 Apr 2020
Edited: Pravin Jagtap on 28 Apr 2020
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)

Answers (1)

Pravin Jagtap
Pravin Jagtap on 27 Apr 2020
Edited: Pravin Jagtap on 28 Apr 2020
Hello Vipul,
Refer to the following documentation for understanding how 'ployxpoly' works:
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.

Categories

Find more on Get Started with MATLAB 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!