summaryrefslogtreecommitdiff
path: root/axial/w_.m
blob: 6d4a2d448b9751efbaac702d8185d2cf901aa9c3 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
%---------------------------------------------------------------
% Returns the beam width given the q factor and wavelength 
% of a Gaussian beam. Note that q and lambda are assumed 
% to have the same units.
%
% SYNTAX: w=w_(q <,lambda>);   
%
% q      = q-factor of the beam at the position where R and w are to
%          be found. q can be a vector
% lambda = wavelength. Can be a vector or scalar.
%
% If both q and lambda are vectors, they must be the same size.
%
%---------------------------------------------------------------
% SYNTAX: w=w_(q <,lambda>);
%---------------------------------------------------------------

function wout=w_(q,varargin)

if nargin>=2, lambda=varargin{1}; else lambda=1064e-9; end

wout=sqrt(lambda/pi .* imag(q).*(1+real(q).^2./imag(q).^2));