'ProgramMode:RUN
'diffusion at slit
7Exp77->A~Z:7Exp77->a~zRad
{500,1}->Dim Mat L'screen
0.1->D'dist of scr 100mm
0.555Exp(-)6->l'lambda 555nm
5Exp(-)6->b'slit width
0.05->w
'every pt at scr...
For 0->i To 383
 i*w/383->Y'mm
 0->y
 Sqrt(D^<2>+Y^<2>)->R
 'every source...
 For b/2->y To (-)b/2 Step (-)b/10
  Sqrt(D^<2>+(Y-y)^<2>)->r
  R-r->s
  yY/Sqrt(D^<2>+Y^<2>)->s'taylor series
  cos (2pis/l)->a'amplitude y to Y
  L[i]+a->L[i]
  Next
Next
Gosub G
0->M'max ampl
3Exp33->mFor 0->i To 383
 L[i]^<2>->A
 A>M=>A->M
Next
For 0->i To 384
 L[i]^<2>->A
 For 0->j To 100*A/M Step 0.5
  Plot i,j
 Next
Next
Gosub S
Stop



Lbl G
ClrGraph:CoordOff:AxesOff:GridOff:LabelOff:S-L-Normal:Rad:'#_SketchThin _
'_AxesOn_
' x0 x\1 dx y0 y\1 dy
ViewWindow 0,383,0,0,100,0
Return 
'Print._X_Y
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