function g = gaussian(x,pos,wid)
% gaussian(x,pos,wid) = gaussian peak centered on pos, half-width=wid
% x may be scalar, vector, or matrix, pos and wid both scalar
% T. C. O'Haver, 1988
% Example: gaussian([1 2 3],1,2) gives result [0.5000 1.0000 0.5000]
g = exp(-((x-pos)./(0.6006.*wid)) .^2);