1 2 3 4 5
function q=wr2q(waste,R,lambda) % convert Gaussian beam waste and radius to q parameter q=1./(1./R-1i.*lambda./pi./(waste.^2)); endfunction