function gaussian,x,width,position

;	Returns Gaussian curve having given width and position.

if n_elements(position) le 0 then position=0.

	if n_params() lt 2 then begin
print,'You must define WIDTH'
return,0
	endif

z=((x-position)/width) > (-5) < 5

return,exp(-z^2*4.*alog(2))

end