  This study proposes an approach toward the first principles electronic
structure calculation with the aid of symbolic-numeric solving. The symbolic
computation enables us to express the Hartree-Fock-Roothaan equation and the
molecular integrals in analytic forms and approximate them as a set of
polynomial equations. By use of the Grobner bases technique, the polynomial
equations are transformed into other ones which have identical roots. The
converted equations take more convenient forms which will simplify numerical
procedures, from which we can derive necessary physical properties in order, in
an a la carte way. This method enables us to solve the electronic structure
calculation, the optimization of any kind, or the inverse problem as a forward
problem in a unified way, in which there is no need for iterative
self-consistent procedures with trials and errors.
