function [waste,radius] =q2wr(q, lambda)
% convert Gaussian beam q parameter to waste and radius
	waste=sqrt (-lambda./pi./imag(1./q));
	radius=(1./real(1./q));
endfunction