'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