clear all; close all;

r = linspace(0,10,1e3); % Integration grid
dr = r(2)-r(1);
R = 1.25; % m
f = besselj(0,2.4048/R*r).*(r<R);

subplot(211);
plot(r,f);
grid on;

subplot(212);
rho = linspace(0,4,1e3);
F = zeros(size(rho));
for i = 1:length(rho)
    F(i) = 2*pi*sum(r.*f.*besselj(0,2*pi*r*rho(i)))*dr;
end
F = F/max(F);
plot(rho,F);