syms x arfa beta gama
arfa=1.12;
gama=0;
beta=0;
a=.20;
b=.04;

for gama=0:.1:.4
for arfa=0:.04:2.25
%I1=vpa(2*a^2*( abs( int(exp(i*beta*x)*exp(-arfa^2*x)*besselj(0,k*x^.5),x,gama^2,1)))^2     );
%I=10*log10(I1)
I1=( exp(-2*arfa^2)+exp(-2*arfa^2*gama^2) )
I=10*log10(I1)
plot(arfa,double(I))
hold on
end
end
%I=int(sin(x)/x,x,0,inf)