Hi all,
I am trying to perform some operations using Bessel function and plot 'v' against 'r' (as defined in the adjacent image).
But when I am defining the function 'F' (from image) and evaluating it at certain points (which are roots of J2), I am getting unexpectedly large numbers and MATLAB is not performing simple mathematical operations on them on it's own.
My program for this is:
bessj2 = inline('besselj(2,x)');
for n = 1:10
z(n) = fzero(bessj2,[(n-1) n]*pi);
end
z(1)=[];
syms x;
F1=diff(besselj(1,x),x);
f1=besselj(0, z) - besselj(1, z)/z;
syms r t;
F=@(x,t)besselj(1,x).*(exp(x.*x*t/100)-1)./(x.*x);
f=eval(subs(F,x,z));
F2=besselj(1,z).*(1./z-1)./z;
F3=@(x,r,t)(besselj(1,x.*r).*exp(-x.*x*t./100).*(f+F2))./(f1.^2);
f3=eval(subs(F3,x,z));
f4=sum(f3);
v=@(r,t)2*f4;
ezplot(v,[0,1,0,10])
Please help me figure out a way to solve this problem.
Thank you.
