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).
syms a e
A = [besselj(0,-1i*a*e), bessely(0,-1i*a*e); ...
b = [-1; -1];
X = linsolve(A, b);
a = linspace(5,9);
e = linspace(0.05, 0.15);
plot(a, e, X)