00001 double precision function GAUSSIAN(axis,pars) 00002 00003 double precision axis, pars(2) 00004 00005 GAUSSIAN = exp(-0.5*(axis-pars(1))**2/pars(2)**2) 00006 00007 return 00008 00009 end