Reduce the double integral

747622d9-e4b3-4504-9f36-47bc1064de8d 20 3 12 decimal

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

This Smath Fourier Quantum is adapted from Mathcad 11.

I was given writen permission from Frank Rioux to publish

the Mathcad 11 work sheet in my former Mathcad web site.

]]>

Fourier Quantum CoC [Circle of Confusion]

jmGiraud: 2019_05_22

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

Reduce the double integral

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

Observe . . .

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

perfectly symmetric vs the built-in integrator.

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

]]>

inner integrand borrowed from collapsed ...  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 ∫CoC 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 20 : 13 mat

px py Φ f a b n ∫CoC :

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/mAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAAAftSURBVFhHzdnLj05ZFAXw+w9I/BMYSYxEhyCCinjEIyKFIJUqz1IK7VmUVytNPbxpSjVSgsSj01TPGegpAwMDnRhKzCSdSHQ6/fvu6tx80T3mnoQ6Z5+1115n733Prfq+oru7+8CBA729vUeOHDl8+PC+ffs6OztXrFixatWqrq4uW+wnTpz4oRzHjx8/evTooUOH9u/fv2vXLsiNGzdu2LBh9erV8IaJJaMtADBgLhzDgAohWuTwAkEKKjQ7GbZIKnp6eqx5MtneunVrS0sLaqQYjx07hsgw4bl3797t27evX79+zpw5y5Yta2trg9+xYwfw9+UwsWS0BQAGzIUj92Y25MACCQcvNAFkEENS4QdQNG3evHny5MnOioXFESsKIdnnz58vH8LkfAAnT57s6+v7sWlYMtpK7oG5cOSOpDoqgBACsQsqdJTZJakgkP+ePXs6OjpmzZq1bds2aWSMJwT0li1bZs+eLQF4earF6dOn+/v7BwcHh4aGzp49e65pWDLaAgAD5sKROxJUCNHmzAIJJ6jQBJBBDGORTPBZtGiReltySBNAyLOzGrt378YlE4I1SzGJjmrEUu2ycOHIHUnY0CIXQiDhBBWagJzcsgBVYH2wc+dOUDgU0aRL5s2bJ72KrS6VoAQbGBiQDOBQO6KRIzHaAmjGc0eCCiFa5FEWBnMCyCCGpILAJUuWMEk1t1OnTgWkJyZNmsTZKXVM6mUIZin8wYMHPTLt7e2epqlTp35XDhNLRlsAYMDRZyCxRIgWuRBJhKBCE0AGMSQVuCTQtg0+2eYwY8YMqdaAZ86ciSbsTmYX6dKlS9etWweQFuZlyzDJI2ILAAyYi62IQ4UQBkAIgZKOhCaDGJIKjeZYkcw5R5k+fbpMynAKh0tR1MhWURT+N2dRnYsXL165cuXq1as/lcPEktHW/7rkhGiRCyGQrXRtCkUMSY3e4tA4aV+fiTSOGzfuC00coFtbW10wDspy6dIlIoaHh0dGRn7+z2C0BQAGzIUjdySovlAmnKDNGhgLUNtJvqfABaMtWOTciVGAul0mTJjgypZ/mbh27dqNGzcqBeZEVMOy0moOzIUjdySoEKJFLoRAwgkqNAFksJDUuE6lTm9ycyBXiDkHj7czORyiBQsWqLqlGiWqQYGlkOyiag7DxJLRFkCQXCzZkaBCiNZSCIGEE1RoAsyJIamg0UKTbtq0yYtCA0pveemc0wruOhcdo8MpSjRdv3798uXLGNnztLskFcIwyZ1iCwAMOMq4I2FHiBZ5ogjHKDQBZBBDUuFHyjlx4sScg/P58+cTdeHChW4aFrWozq0E8uzcAFUKq6olMbYAwIAD4I4EFUK0AEIIxJKaEJBGJ6khi0bN77IhE0ugksmo5CwOGk2eMhRr165Fqp1v3bp17969Bw8ePH78+JdymFgy2gIAA+bCMcpQIUSLXIikgEVoAhiTsIYsLynyQZOqlM8STrtgRGeYaEm/FGggibl//z4dT548GRsb+61pWDLaAgAD5sKxmQctciFSyiTMkgxiGrLI1BPuPTJRQNAuk54O/2veFMgpUUyZMsXu7du3Hz58KHazjl/L0azSHAyYC0fuSFJotFUIu4IKTQAZxDR6S9vrU9mzUGxtQbiH1m3LQcOGiCcMamEePXr09OlTgf2vakp28+ZNfWOYWDJWAGAuHLkjySHRIhdCIOEEFZoAGGIaTyKNym9bnkGzvXjxYufQHHKORc/qEv2LVAISUm7u3r2rMxwab3lnDZtYMtoCiDIuHLkjQYUQLXIhBEo6hCaADGJIKtwWM2fOVFH1tieZLF7jJjk9RsHmzp174cIFmUjthBwdHfXgkKJqz58//70cJpaMtgCijAtH7khQIUx2hRBIOBOhCSCDGJbCK3P58uUuurxqpNRLilHCsTuZhtCDzopLdUSVAMkQ+M6dO8+ePXv16tWbN2/+KIeJJaMtADBgLhy5I0GVDkMuhEDCpY4EkEEMY+Ft5aGApteeN4MCu1p0ZfncNLLtPS/PYiRVYji6/hD+9evX7969e//+/YdymFgy2gIAA07CuCNBld4whBBIOEGFJoAMYkgqpFHrqSvthrtkzZo1ms4FzTMVnDZtmqfaheTchoqMHz/e5OXLl0RQ8/Hjxz/LYWLJaAsADLh0GuOOBBVCtMiFEEi4XGAGGcSQVKxcuVK21TX3B4Q0asPIcmKF93QgyrkVxVXp3PKhZNJDyqdPn/4qh4kloy0AMOA8ItyRoEKINrIEEi6yCCCDGJIKN5iL3wOcy90L3NMhmbn9NKbCK7ZucEOmHEjZX7x48fbtW7mRpM+fP/9dDhNLRlsAYMApPXckqBCyIxdCIOEEFZoAMoghqa6yalrEmrZ8TS8IxZbG2l2nNX351PRVrek8HQqcbbxSalvrgeoPLEglHMYD/JV+sbGo46+BelBF3WD1+qXZPzJlD45M2nPdu0sYQVmcEp1TYuSv/LpEc7gqlcyFpGp0GCaWjLYAwIC5cOSOBBVCtMhzXQnHIjQBjMT8K0smPTj1+oOMzCSsXn++anvJtHBb1OiPfXn2AteSNCq5p8NDy8LB4ZTAmwHR1/5oJOUENUwUohYfJGl+aKnjwCqBXlLf/mO3mn5IWdOPdGXSa7x2H4BLuKeAwHp9XeCIELX7cgVCA/KElt66fBWl2H4QGGUumFp8ceeRkUYm65zv23/N2d39D0W0XEo03/eBAAAAAElFTkSuQmCC CoC

1 time t0 - min 11

hatch

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 :

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

x f x 2 Ψ abs (2^:

x g 0 :

coordinate-space prob. density

f g 2 - 2 20 hatch 2 1 sys