
mu0=4*Pi*10^-7;
c=3*10^8;
resistance=mu0*c/(2*Pi)*NIntegrate[Cos[Pi/2*Cos[theta]]^2/Sin[theta],{theta,0,Pi}]
            

p=10^4;
intensite=Sqrt[2*p/resistance]
            

ps[theta_]:=Cos[Pi/2*Cos[theta]]^2/Sin[theta]^2;
            
PolarPlot[ps[theta-Pi/2],{theta,0,2*Pi},PlotRange->{{-1,1},{-1,1}}]
