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.