Hello there,
I would like to plot the real part of this symbolic function I have just solved.
Parameters 'a' and 'e' are varying between 5 and 9, and 0.05 and 0.15 respectively.
Unfortunately, it does not work.
Here you go with my code ( solution is a ratio of sum of Bessel functions of order zero).
clcclearset(groot,'DefaultAxesFontSize',14) set(groot,'DefaultLineLineWidth',1.5)syms a eA = [besselj(0,-1i*a*e), bessely(0,-1i*a*e); ... besselj(0,-1i*a), bessely(0,-1i*a)];b = [-1; -1];X = linsolve(A, b);disp(X)a = linspace(5,9);e = linspace(0.05, 0.15);plot(a, e, X)
Best Answer