Fraunhoferovy ohybove jevy na ctvercovem otvoru

Osa otvoru je totozna se souradnou osou x, strany otvoru jsou rovnobezne se souradnymi osami y, z.
Intenzita ve Fraunhoferove ohybovem obrazci pri ohybu na ctvercovem otvoru je dana vztahem
I=I\(TraditionalForm\`\(\_0*\)\)(((((sin(kβa))/k)βa))^2)*(((((sin(kγa))/k)γa))^2),
kde k=((2  π)/λ) je vlnove cislo,  
      a je polovina sirky ctvercoveho otvoru a β, γ jsou uhly, ktere svira osa ctvercoveho otvoru se spojnici stredu otvoru s bodem, ve kterem je intenzita v ohybovem obrazci pocitana, ve smerech souradnych os y, z
Vykresleme intenzitu I v ohybovem obrazci za nasledujicich predpokladu: I\(TraditionalForm\`\(\(\_0 = 1, \)\ \)\)kβa a kγa se meni v mezich -8,+8
Intenzitu vykreslime postupne ve forme 3D barevneho grafu ve trech ruznych stupnich rozliseni - 70, 110 a 200 bodu/sirka grafu a ve forme "hustotniho" grafu, ktery se blizi fotografiim ohybovych jevu

Plot3D[
      (Sin[x]/x)^2*(Sin[y]/y)^2, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 70, Axes -> True, PlotRange -> {0, 1}]

Plot3D[
      (Sin[x]/x)^2*(Sin[y]/y)^2, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 110, Axes -> True, PlotRange -> {0, 1}]

Plot3D[
      (Sin[x]/x)^2*(Sin[y]/y)^2, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 200, Mesh -> False, Axes -> True, 
      PlotRange -> {0, 1}]

DensityPlot[
      (Sin[x]/x)^2*(Sin[y]/y)^2, {x, -8, 8}, {y, -8, 8}, 
      PlotPoints -> 200, Mesh -> False, Axes -> False]