%Refractive Indexnc=1.0;                     %Region I Coverng=2.0;                     %Region II RI Cylindernr=2.0;                     %Region II RI Backgroundns=4.0;                     %Region III Substrate%Setup Sourcelambda0=3;                  %Freespace wavelengthk0=2*pi/lambda0;            %Freespace wavenumberrad=pi/180;                 %Radianstheta0=0;                   %Polar anglephi0=0;                     %Azimuthal anglepsi=0;                      %Polarizatin angle 0->TE, 90->TMtheta0=theta0*rad;phi0=phi0*rad;psi=psi*rad;depth=0.5;