
                f[x_]:=Sin[2*Pi*x]; 
                g[x_]:=1/x;
                s = "bonjour";
            
Plot[f[x],{x,0,5}]
A={{1,2,3},{4,5,6},{7,8,9}}
 
a=0.0001; (*diametre*)    
lambda=0.5*10^-6; 
k=2*Pi/lambda;  
intensite[theta_]:=(2*BesselJ[1,k*a/2*Sin[theta]]/(k*a/2*Sin[theta]))^2; 

Plot[intensite[theta],{theta,-0.02,0.02},PlotRange->{{-0.02,0.02},{0,1}},AxesLabel->{"theta","i"}]

n=101;   
theta = Table[-0.02+k*0.04/n,{k,n}];   
i = Map[intensite,theta]; 
m = {theta,i};   

