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