Hi nelson
Have you tried different parameters for the guassian and lorentzian components?
prm pr1 0.007
prm pr2 0.007
prm pr3 0.007
prm pr4 0.007
Zero_Error(zero, 0.002)
Simple_Axial_Model(axial, 1.07693`_0.04475)
lor_fwhm = pr1 Tan(Th) + pr2/Cos(Th) ;
gauss_fwhm = pr3 Tan(Th) + pr4/Cos(Th) ;
Adrian