A recently proposed parameter free version of a Lagrangian theorem-based density functional approximation (LTDFA) [S. Zhou, Phys. Lett. A 319, 279 (2003)] for hard-sphere fluid is applied to hard-core attractive Yukawa model fluid by dividing bulk second-order direct correlation function (DCF) of fluid under consideration into hard-core part and tail part. The former is treated by the parameter free version of the LTDFA, while the tail part is treated by second-order functional perturbation expansion approximation as done in a recent partitioned DFA [S. Zhou, Phys. Rev. E 68, 061201 (2003)]. Two versions of mean spherical approximation (MSA) for the bulk second-order DCF are employed as input, one is the less accurate plain MSA whose tail part of the second-order DCF is strictly independent of a density argument, the other is the more accurate inverse temperature expansion version of the MSA whose tail part is not strictly independent of the density argument. Calculational results indicate that prediction based on the plain MSA is far more accurate than that based on the inverse temperature expansion version of the MSA. The reason is considered to be that the partitioned DFA requires that the tail part is highly or completely independent of the density argument, the plain MSA, by assuming that the tail part is exactly the potential itself, embodies all of the nonlinearities into the hard-core part which can be treated satisfactorily by the parameter free version of the LTDFA. The present investigation results in a universal method for constructing DFA for nonuniform any nonhard-sphere interaction potential fluids. (C) 2004 American Institute of Physics.