A progress report on the development of numerical and rigorous solutions of the effective two-band two-dimensional Hubbard model of the high-temperature superconductivity in cuprates [N.M. Plakida et al., PRB 51, 16599 (1995)], within the equation of motion technique of the thermodynamic Green function approach, is given. Four important consequences of the rigorous mean-field solution are worth mentioning: (1) Occurrence, both in the normal and superconducting states, of the spin-charge separation, conjectured by P.W. Anderson on intuitive grounds. (2) Origin of the static d-wave pairing in the hopping singlet conduction. (3) Complete suppression of the s-type superconducting pairing by the crystal symmetry in square two- dimensional lattices and occurrence of a residual s-type term in rectangular lattices, in agreement with the experiment. (4) The energy band hybridization in the superconducting state points to an overall displacement of the energy spectrum toward lower energies, as pointed by very precise optical measurements.