1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
%---------------------------------------------------------------
% PROGRAM: AiryI
% AUTHOR: Andri M. Gretarsson
% DATE: 8/7/03
%
%
% SYNTAX: I=AiryI([a,w0,L,lambda <,N>],r);
% <,...> indicates optional argument
%
% Returns the intensity of an Airy disk as a function of radius
%
% a = radius of aperture
% w0 = radius of Gaussian beam at aperture
% L = propagation distance (axially) from aperture
% lambda = wavelength
% N = order to which to carry the calculation
% Last updated: 8/7/03 by AMG
%
%---------------------------------------------------------------
%% SYNTAX: I=AiryI([a,w0,L,lambda <,N>],r);
%---------------------------------------------------------------
function I=AiryI(params,r);
I=AiryE(params,r).*conj(AiryE(params,r));
|
