A continuous time dynamic system of an unsteady state fugacity model is presented. Properties of this model as stability are studied. In order to evaluate numerical results a discretization preserving the stability and yielding the positivity property of the model is used. Finally, algorithms to determine the values of the fugacities, the concentrations and the dissipation time are given. The above study is illustrated with numerical results in a theree compartmental environmental system