Re: Fit

*To*: mathgroup at smc.vnet.net*Subject*: [mg76173] Re: Fit*From*: dh <dh at metrohm.ch>*Date*: Thu, 17 May 2007 06:00:14 -0400 (EDT)*References*: <f2ejm8$d2$1@smc.vnet.net>

Hi Mathieu, it depend what you want to do. If it is important that the approximation is good near a point, e.g. 0, one can expand both function in a series and equate the coefficients: Series[f[x,y],{x,0,2},{y,0,2}]//Normal//N Series[Gaussian2D[x,y,rad,amp],{x,0,2},{y,0,2}]//Normal//N in your case we would euqte the zero and quadratic (x^2,y^2 have the same factor) terms. On the otherhand, if the approximation is to be used over a region, one can do it numerical. Choose enough points in the region, calculate the function values at these points and use e.g. FindFit to fit the Gaussion. hope this helps, Daniel Mathieu G wrote: > Hello, > I have an expression of the form: > > 8000000 Sqrt[2/\[Pi]] (-Erf[(-300 + x)/(4000 Sqrt[2])] + > Erf[(300 + x)/(4000 Sqrt[2])]) (-1/2 Sqrt[\[Pi]/2] > Erf[(-300 + y)/(4000 Sqrt[2])] + > 1/2 Sqrt[\[Pi]/2] Erf[(300 + y)/(4000 Sqrt[2])]) > > how can I make a fit using a 2D Gaussian function of the form > > Gaussian2D[x_, y_, Radius_, Amplitude_] := > Amplitude Exp[-1/2 ((x/Radius)^2 + (y/Radius)^2)]; > > Thank you for your help! > MG >