cond order Taylor expansion coefficients of these profile functions, restricted only by the symmetry properties of the source, and can be trivially expressed by re-scaling the earlier fit parameters. The above is the most direct form of the Buda-Lund model. However, different combinations may also be used to measure the flow, temperature and fugacity profiles [@Csorgo:1995bi; @cs-rev]: $ H_t   \equiv  {b}/{\tau_0} \, = \, \ave{u_t} / R_G
    \, = \, \ave{u_t^\prime} / R_s$ , $H_l   \equiv \gamma_t /\tau_0$, where $ \gamma_t = \sqrt{ 1  + H_t^2 r_t^2}$ is evaluated at the point of maximal emittivity, and & = & = \_r = ,\
& = & = \_s = .

Buda-Lund fits to Au+Au data at $\sqrt{s_{\NN}}= 200$ GeV 
==========================================================

In this section, we present new fit results to BRAHMS data on charged particle pseudorapidity distributions [@Bearden:2001qq], and PHENIX data on identified particle momentum distributions and Bose-Einstein (HBT) radii [@Adler:2003cb; @Adler:2004rq] in Au+Au collisions at $\sqrt{s_{\NN}}= 200$ GeV.

The analysis codes and methods are identical to the ones used to fit the BRAHMS [@Bearden:2001xw], PHENIX [@Adcox:2001mf; @Adcox:2002uc], PHOBOS [@Back:2001bq], and STAR [@Adler:2001zd] data in 0- 5% most central Au+Au collisions at $\sqrt{s_{\NN}} = 130$ GeV, see ref. [@ster-ismd03]. The applied Buda-Lund 1.5 fitting package can be downloaded, together with the detailed fit results, from ref. [@Csorgo-blhome]. This calculation determines the position of the saddle point exactly in the beam direction, but in the transverse direction, the saddle point equations are solved only approximately, as summarized in ref. [@cs-rev].

The new results for $\sqrt{s_{\NN}} = 200$ GeV Au+Au collisions in the 0-30% centrality class are shown in the first column of Table 1. For comparison, we also show the results of an identical fit to $\sqrt{s_{\NN}} = 130$ GeV Au+Au collisions in the 0- 5% centrality class.

  ------------------------------ ------ ----------- ------- -----