'ProgramMode:RUN
'diffusion at parabel.
7Exp77->A~Z:7Exp77->a~z:Rad:'#_Mat _0
{500,1}->Dim Mat I'screen
1->f'Focal length
0.555Exp(-)6->l'lambda 555nm
0.2->A'aperture
0.000005->w
AA/(4f)->G
10->S'steps in each direction
For 0->i To 383'every pt at scr..
 i*w/383->X
 0->Z
 For (-)A/2->x To A/2 Step A/(S+1)  'every source x
  For (-)A/2->z To A/2 Step A/(S+1) 'every source z
   Sqrt(xx+zz)->r
   If  r<=A/2
    rr/4f->y
    G-y->g
    f-y->Y
    G+Sqrt(X^<2>+f^<2>)->Q
    g+Sqrt((x-X)^<2>+Y^<2>+(z-Z)^<2>)->q
    Q-q->s
    cos (2pis/l)/qq->a'amplitude x to X
    I[i]+a->I[i]
    '{X,x,z,r,y,g,s,a}_Disps_
    
   IfEnd
  Next
 Next
Next
Gosub G
0->M'max ampl
For 0->i To 383
 I[i]^<2>->a
 a>M=>a->M
Next
For 0->i To 383
 I[i]^<2>->a
 F-Line i,0,i,100*a/M
Next
Gosub S
Stop



Lbl G
ClrGraph:CoordOff:AxesOff:GridOff:LabelOff:S-L-Normal: Plot/Line-Color RGB(0,0,230):Rad:'#_SketchThin _
'_AxesOn_
' x0 x\1 dx y0 y\1 dy
ViewWindow 0,383,0,0,100,0
Return 
'Print.XY
Lbl P
Text 1,9,"                               "
Text 1,9,X
Text 1,179,Y
Return 
'Stop
Lbl S
Plot :Text 1,375,">"Disps
While Getkey<>31:WhileEnd'exit=47
Return 
'stop at test
Lbl s
Getkey<>47=>Return 
Stop