Reduce the double integral

68828ce6-c809-4e4f-a11a-3bce97235d3e 8 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

Fourier Quantum CoC [Circle of Confusion]

Optimiz => SYMBOLIC M lo hi Scale hi lo - (M M Min - M Max M Min - / (* lo + :

Fourier transform of position wave function into 
momentum representation to get diffraction pattern 
@ the detection screen ... recast Φ(px,py) . px py Φ i px x * y py * + (* - exp y r 2^x 2^- sqrt - r 2^x 2^- sqrt int x r - r int ≡

1. Reduce a double integral

pass the double integrand to maple.
evaluate maple SYMBOLIC i px x * y py * + (* - exp y r 2^x 2^- sqrt - r 2^x 2^- sqrt int x r - r int

Symbolic

i px x * y py * + (* - exp y r 2^x 2^- sqrt - r 2^x 2^- sqrt int x r - r int maple i i px x * r x - (r x + (* sqrt py * + (* - exp - i px x * r x - (r x + (* sqrt py * - (* - exp + (* py / - x r..r - = int

i i px x * r 2^x 2^- sqrt py * + (* - exp - i px x * r 2^x 2^- sqrt py * - (* - exp + (* py / - x r..r - = int

isolate/recast the integrand i py / i r 2^x 2^- sqrt py * x px * + (* - exp * i py / i r 2^x 2^- sqrt - (py * x px * + (* - exp * -

recast/insert the limits i i px x * r 2^x 2^- sqrt py * + (* - exp - i px x * r 2^x 2^- sqrt py * - (* - exp + (* py / - x r - r int

assign the integral
wrt the dependents
Φ(px,py) ... acc 50
.... SYMBOLIC .... px py Φ i i px x * r 2^x 2^- sqrt py * + (* - exp - i px x * r 2^x 2^- sqrt py * - (* - exp + (* py / - vectorize x r - r int :

Fourier transform of position wave function into 
momentum representation to get diffraction pattern 
@ the detection screen ... recast Φ(px,py) . px py Φ i px x * y py * + (* - exp y r 2^x 2^- sqrt - r 2^x 2^- sqrt int x r - r int ≡

β 0.125 :

x q 11 x abs β^+ / :

The mapping function

x q 0.275 0.625 0.275 0.25 2 2 mat 2 1 sys

Observe . . .

For the double integration, this Simpson is 'n' dependent,

perfectly symmetric vs the built-in integrator.

Set 'n=24' representative demo of the CoC project.

]]>

inner integrand  ...  extracted from Maple Φ(px,py) x px py f i i px x * r 2^x 2^- sqrt py * + (* - exp - i px x * r 2^x 2^- sqrt py * - (* - exp + (* py / - :

Numeric

f a b n ∫∫ h b a - n / : 13 / h 2 / a px py f b px py f + (* h a k h * + px py f k 1 n 1 - sum * + (* 23 / h * a h 2 / - k h * + px py f k 1 n sum * + eval 21 line :

a r - : b r : n 24 : 13 mat

recast Simpson double integral px py Φ f a b n ∫∫ :

r 0.275 : (pin hole radius ≡ Δ 51 : (size momentum ≡ N 51 : (canvas size ≡ N Δ ≥ (N<Δ not diffraction ≡ px py 12 mat position momentum ≡ N odd escape divide zero j 0 N range : px Δ - 2 Δ * j * N / + : py Δ - 2 Δ * j * N / + : 12 mat 81 line

..... NUMERIC ..... CoC m 1 N range n 1 N range Fo m 1 - n 1 - el px m el py n el Φ abs (2^: for for fmap Fo q vectorize : fmap - 0255 Scale 31 line :

t0 1 time :

iVBORw0KGgoAAAANSUhEUgAAADIAAAAyCAIAAACRXR/mAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAAAhCSURBVFhHzdnLT5ZHFAbwd6N/hImaVDcYWIAXjIlRQyRWjSbeEBQREfEu4rVqvdIqXgqIFYhETUBoitS9bHTtho2Lblya2J1JFyZdtL/vPcnki24bfScRZ86c8zzPnHNmvk/MWltbjx49evLkyVOnTvl5+PDhPXv2rF69ev369baOHTt2+vTp8+fP/5CPc+fOnTlzhpuQAwcOtLS0NDY2btmyZePGjd/nw8SS0RaHQBYiMBBAAQQLHAUidEiTACG2Mh7WXENTU1NTVVXVtm3bOjo6WM6ePQvIMLE8cuSIGMRLly5dtWqVye7du/ft29fe3s7fMLFktMWBm4kQgV+i8UeEDmkoI4MYkrKk6dChQzt37pw/fz4nKOyOmCDssq9YsUJKZCJQbEnAxYsXfywblowRxY2zEIHCgUQUWOAoELEjRR27oawky+LgwYPbt293OCjHjx+X9vxUpYSLdPo4t/N1dXVhvXz58rVr127cuNHd3f3TF4PRFgdunIUIjByDAgg28BGhQ2qLADKIKcnq7OwUs3Xr1pUrV8q/ZfSBGHD6I/rGhLdMXL16FSXun/ORdKQRlrTLWYhA4UCi/0yAo0CEDilqAsggxjILTfqgra2NUq4XLlzw01xPLFmyZNeuXSdOnGCEnlhTMvBxVhfnNqJAjCmd4S9QOBBQAMEC/4yOADJCWbZjx47ly5cTG9WJzuAksQsWLGCXcB1z/fr1PBfdyCwdUfL379/f0NCwdu3aioqK7/JhYsloiwM3zkIiFoglQLDAUSBKpASwE0NSNmvWLFdU3mxcunSJ9tDkQFIN98qVKw4dgjgIdqy6urpNmzbhdvrU+IFuyWiLAzfOQgSGOFAAwQJHEcqQchBLBjEkZRrNhZTG0KQD9u7du3jxYplkTFiKokY45s2bJ/9OzKI6PT09d+/e/aVsWDLa4sCNsxCBwlnSCYGjQIQOaShjJIakjGrBcVYT8a6rANqTJgEqolXdI+hqcfv27d7e3vv37z948ODXLwajLQ7cOAsRKBwIqKQMBSJ0SMs1kJRJOI3ESr776aLKPws4JwbB1fOTZZmLLf+Y+vr6BgYGkgJzItKwTFrNOQsRKBwIKIBggaNAhA4pagLIYCEpo45qxfbqO1Bzc7NtR3G9ncnhvHIeQx1gee/evRAU+VAy9WJ3y/SyEbeV0RaH0CdEIDsQUADBWqJAhA4pagLIIIakjJVAAe6tJjUBHU+OVuCq/GpvqShojP7+/jt37jgru6jq6movtUMbJpaMtjhw4xxRwoGwAwQbrYkIHVLUokyIISnzl9rL4aJFi/x0jvB2Djdo2bJlbgeLWoB2bmlQAl3iHnEoT2F5Ymxx4MZZCLtwIKAAguUQNWFBmgQQQ1JJFo0+831qyp4jcoXLHp8GcCNPoN0yldq8eTNQHTM4OPjkyZOxsbHx8fHf8mFiyWiLAzfOQgSGMlAAwQJHEaVEipoAMiJhmXJqRg+//o9UGTLsjngVTZwenGGiJX3g61+Jefr06cTExOTk5NTU1B9lw5LRFgdunIUILMcBCxxFlNJATQAZxJCUSZqeWLNmjdwqM4/Q7nb4qXmjQE4JYu7cuc46PDz87Nkz3OU6fs9HuUpzbpyFCBQOJAoNNlGgQ4qaADKIIan0DUKf+l6h0XjIqpzb88pZ6gZA7pRInwmgh4aGVOr58+eI/ZQVJXv48CFPw8SSMTlwFiJQOBBQAHkCR4EIHVJLAsgghqTSu+Vbou14TmKbaufQHHIORdp1if611DpBKTePHz++efNmTU2NQDfOMLFktMUhlAkRKBwIKICWwFEginQIJIAMYkrvllrOmTMnbizVkun516d6Qs6dXs7ZKysrb926pV2idigfPXqkOUhRtenp6df5MLFktMUhlAkRKBwIKIBggaNAhA4pOwFkEFN6eAmsr6+P62pPSn0UMEbCI96BnNVcdSIBKoV4ZGSEjjdv3rx9+/bPfJhYMtriwC1SK1A4EFBx2mgMROiijvEkEcOY+bRyKTS/NNqTSUvxulKqI9s+5+VZXSJVOBxdStDPzMy8e/fu/fv3f+XDxJLRFgducRKBwoGAit4wUCBChxQ1AWRYkpT52qX1XGAbhrd/w4YNqu6BFhkVXLhwoZ4YHR1VIEMOZs+ebSIxRFDz8ePHv/NhYsloiwM3znnQlHAgoACCBY4CETqkOXk3GcSQlK1bt86DIY1SanjKfEWk2gMoUqp1gNth4k5FBT2VEv7y5Uslkx5SPn369E8+TCwZbXHgxjnqKBwIKIAmwFEgQoc02MkghqTMXfDmarf4HOBXW1vLz8lEKoSHTrFdH4/4ixcvlIPRDX/16pVm+vDhgyQR9G8+TCwZbXHgxlmIQOFAQAFkBI4CETqkqAkggxiSiiqroEUsaMsX9IFQ7CI+pwX98CnoR3VBv9gU9Gug5lfLwn1p9pQV8Z8Y/kgav2L9g0yj+YvGYv3zVfa8/crJtUD/2PdIEKglbcuh2+HQLM7qNErgOfkGvxrR/9RxNUx8FBTiF0luhOefRgHEqr1XTvkFMCZlWgG6+K/0a7eC/pLSZ0IRf6Urk752+YroWJykMfrOXP4dyL11esa4ZapjODQyRYHIWYbUyDCxZLTFgVv4CxQOBBRAsMA/oyOADGJIKup/F/Cm1P3UaC6qxEq+c4sxJNzFdo9sfdX/XPGHk4VXzvcK19VLLZLdOUSqS+yyewwdF4rnhzE1u45Jw5IxorhxFiJQOJCIAgscBSJ2pKhjl5iSLBcyKYPCyQe+T03nSxBJHBQ9EefWByZOL////3/cifHqM1n7adsV/cb/zdna+h8p1B6nER7C0QAAAABJRU5ErkJggg== CoC

1 time t0 - min 10.2

hatch, scale

You need that one for scaling grayscale image M lo hi Scale hi lo - (M M Min - M Max M Min - 106 -^+ / (* lo + :

Erzeugt eine Matrix für die Darstellung vertikaler Linien zwischen den Funktionen f(x) und 
g(x) im Bereich von x1 bis x2 mit 2*n Linien. Generates a matrix for 2*N vertical lines (hatching) 
between functions f(x) and g(x) in x-range from x1 to x2 f g x.1 x.2 N hatch dx x.2 x.1 - N / (: h x.1 x.1 f x.1 x.1 g 22 mat : x x.1 x.2 dx - x.1 dx + range h h x dx 2 / + x dx 2 / + g x dx 2 / + x dx 2 / + f x dx + x dx + f x dx + x dx + g 42 mat stack : for h 41 line :

Airy CoC [Circle of Confusion]

x r Ψ 12 r * sqrt / x r - ≥ (x r ≤ (* 0 cases :

r 2 :

u r Φ 12 π * sqrt / (i u * x * - exp x r Ψ * x r - r int * :

Airy x 5 - 5 4.95 - range : j 1 x rows range Airy j el x j el r Φ abs (2^eval : for x Airy augment 31 line :

0 r Φ 1.57 r Φ abs 21 mat 0.798 021 mat

x f x r Ψ abs (2^:

x g 0 :

CoC 1.57 0.25 1.57 0.125 - 22 mat :

coordinate-space prob. density

f g r - r 20 hatch 0 1 0 0.125 - 2 2 mat 2 1 sys

Airy CoC [Circle of Confusion]

Airy CoC 2 1 sys