Planar Imaging of a square

generic function w/o oblique angle

\(\(intensity[y_, z_, d_, w_, l_, μ_] := 
	
       If[Abs[y] < \((d*l)\)/\((2\ z + \ w)\), Exp[\(-μ\)\ w], 
		
         If[Abs[y] > \((d*l)\)/\((2\ z - w)\), 1.0, 
			
           Exp[\(-μ\) \((d*l/\((2\ Abs[y])\)\  + \ w/2\  - \ z)\)]]];
     \)\)

\(Plot[intensity[y, 2, 4, 1, 1, 1], {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotStyle -> Thickness[0.01]}]\)

generic function w oblique angle

\(\(intensityobl[y_, z_, d_, w_, l_, μ_] := 
	
       If[Abs[y] < \((d*l)\)/\((2\ z + \ w)\),
         \((d/Sqrt[d^2 + y^2])\)^3
            Exp[\(-Sqrt[1\  + \ \((y/d)\)^2]\)\ μ\ w], 
		
         If[Abs[y] > \((d*l)\)/\((2\ z - w)\), \((d/Sqrt[d^2 + y^2])\)^3, 

           		\((d/Sqrt[d^2 + y^2])\)^3	
             Exp[\(-Sqrt[1\  + \ \((y/d)\)^2]\) μ
                 \((d*l/\((2\ Abs[y])\)\  + \ w/2\  - \ z)\)]]]; 
\)\)

\(Plot[intensityobl[y, 2, 4, 1, 1, 1], {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotStyle -> Thickness[0.01]}]\)

\(Plot[intensityobl[y, 2, 4, 1, 1, 1], {y, \(-10\), 10}, 

       		{PlotRange -> {0, 1}, \ PlotStyle -> Thickness[0.01]}]\)

Variation of L

w/o oblique

\(Plot3D[intensity[y, 2, 4, 1, L, 1], {L, 0, 1}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

w/ oblique

\(Plot3D[intensityobl[y, 2, 4, 1, L, 1], {L, 0, 1}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

Variation of μ

w/o oblique

\(Plot3D[intensity[y, 2, 4, 1, 1, μ], {μ, 0, 1}, {y, \(-2\), 2},
       
		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

w/ oblique

\(Plot3D[
       intensityobl[y, 2, 4, 1, 1, μ], {μ, 0, 1}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

Variation of w

w/o oblique

\(Plot3D[intensity[y, 2, 4, w, 1, 1], {w, 0, 1}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

w/ oblique

\(Plot3D[intensityobl[y, 2, 4, w, 1, 1], {w, 0, 1}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

Variation of d

w/o oblique

\(Plot3D[intensity[y, 2, d, 1, 1, 1], {d, 2.5, 4}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

w/ oblique

\(Plot3D[intensityobl[y, 2, d, 1, 1, 1], {d, 2.5, 4}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

Variation of z

w/o oblique

\(Plot3D[intensity[y, z, 4, 1, 1, 1], {z, 0.5, 3.5}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)

w/ oblique

\(Plot3D[intensityobl[y, z, 4, 1, 1, 1], {z, 0.5, 3.5}, {y, \(-2\), 2}, 

       		{PlotRange -> {0, 1}, \ PlotPoints -> {64, 64}}]\)