Surface Plasmons Single Scattering

Initializations

Clear [α,β,e1,λ,C2];
α = 3;
e1=-16;
λ=600*10^-9;
β = ((2π)/λ)(((((e1 )/(e1+1))))^((1/2)));
um=10^-6;

C2[x_,y_]:=Sqrt[((α^2)/(8 π β Sqrt[x^2+y^2]))];
Intensity[x_,y_]:= 1+C2[x,y]^2+2C2[x,y] Cos[β Sqrt[x^2+y^2] - β x + π/4];

Generation of graphics

rad=2.1*um;
X[t_]:=rad*Cos[t];
Y[t_]:=rad*Sin[t];

GridX=N[Table[ArcCos[((π/4+β rad-2π y)/(β rad))],{y,1,7,1}]];
For[i=1,i≤7,i++,AppendTo[GridX,Evaluate[Table[-ArcCos[((π/4+β rad - 2π y)/(β rad))]+2π,{y,1,7,1}]][[i]]]];

P1=Plot[{Intensity[X[t],Y[t]],(β Sqrt[X[t]^2+Y[t]^2] - β X[t] + π/4)/(2π)},{t,0,2π},GridLines->{GridX,{1,2,3,4,5,6,7}},Frame->True,FrameLabel->{"Angle","Units of 2π"},PlotLabel->"Circular Intensity Distribution",FrameTicks->{{{1,2,3,4,5,6,7},None},{{0,π/2,π,3π/2,2π},None}},ImageSize->300,FrameStyle->Thickness[0.004]]
Export["D:/School Documents/AAU/Projects/P7 - Scattering and Focusing of Surface Plasmon Polaritons/Report/Pics/IntDistCircular.pdf",P1];

X[x_]:=x*um;
P2=Plot[Intensity[X[x],0],{x,-4,10},PlotRange->{0,2},Frame->True,AxesOrigin->{0,1},FrameLabel->{"Distance (μm)","Intensity (AU)"},PlotLabel->"Intensity Distribution Along X",ImageSize->300,FrameStyle->Thickness[0.004],FrameTicks->{{Automatic,None},{Automatic,None}}]
Export["D:/School Documents/AAU/Projects/P7 - Scattering and Focusing of Surface Plasmon Polaritons/Report/Pics/IntDistAlongX.pdf",P2];

P3=Plot[Table[((y^2-λ^2(n-1/8)^2)/(2 λ(n-1/8))),{n,1,15,1}],{y,-5*um,5*um},Frame->True, FrameLabel->{"Distance Y (μm)","Distance X (μm)"},PlotLabel->"High Intensity Fringes",ImageSize->330,FrameStyle->Thickness[0.004]]
Export["D:/School Documents/AAU/Projects/P7 - Scattering and Focusing of Surface Plasmon Polaritons/Report/Pics/HighIntFringes.pdf",P3];