Gaussian basis

Create a 'change of basis' matrix Input is: Lmax number of basisfunctions r radial axis lambda wavelength omega0 angular frequency n refractive index deltar stepsize in the radial direction

Output is: basis 'change of basis' matrix

function basis = GaussianBasis(Lmax,r,z,lambda,omega0,n,deltar)

    % Allocate basis
    basis   = zeros(Lmax,size(r,2));

    % Calculate the Gaussian modes
    for l=1:Lmax
        basis(l,:) = conj( r.* 2.*pi .* deltar .* ...
            GaussianMode(l,r,z,lambda,omega0,n) );
    end

end