NAME: rfgfit PURPOSE: Fit a radial gaussian to the input data (no linear term, width fixed). DESCRIPTION: Fit the equation y=f(r) where: F(r) = A0*EXP(-z^2/2) + A2 and z=r/A1, A1 is constant A0 = height of exp, A1 = sigma (the width). A2 = constant term, The parameters A0, A1, A2 are estimated and then CURVEFIT is called. The gaussian is assumed to be centered at r=0. CATEGORY: Function fitting CALLING SEQUENCE: Result = RFGAUSSFIT(R, Y [, A]) INPUTS: R: The independent variable. R must be a vector. Y: The dependent variable. Y must have the same number of points as R. OUTPUTS: The fitted function is returned. OPTIONAL OUTPUT PARAMETERS: A: The coefficients of the fit. A is a three-element vector as described under PURPOSE. COMMON BLOCKS: rf_gauss - Contains the constant gaussian width to use. SIDE EFFECTS: None. RESTRICTIONS: PROCEDURE: The initial guess of the intensity is the value of the point with the smallest r value. The first point that drops below half of the first is the guess at the width. MODIFICATION HISTORY: DMS, RSI, Dec, 1983. Rewritten to change initial guess, MWB, Lowell Obs., Dec. 1992. Rewritten to remove background function, MWB, Lowell Obs., Jan. 1993.