A mathematical eutrophication model is proposed, which was obtained by combining a pre-existing three-dimensional model based on eddy-diffusive transport, embodying tidal agitation, and a ‘trophic dynamic reactor’ interconnecting eight intensive physico-chemical and biological variables. Parameters included in the phytoplankton growth equation, K N, K P, μ max, were obtained through several ‘in situ’ investigations performed under varying lagoon conditions. The usual statistical treatments for analyzing collected data, combined with nonlinear multiple variability, were found unsuitable. Therefore, three statistical methods were elaborated: algebraic resolution, including analyses of normal and log-normal distributions of the results; multiple regression analyses with constrained minima through a secular equation of the sixth degree; and analysis of the hypersurfaces of variance and reconstruction, from computer maps, of three-dimensional profiles with paraboloid-like base. Comparisons of results obtained with these methods permit delimitation of variability ranges of constants. Further, the values adopted in the model to simulate trophic equilibria, reproduced real lagoon conditions.