Poisson equations were modified for the electric field and electric field gradient. The solvent density and dielectric polarization were solved. The electric potentials at the solute atoms were solved. The solvation free energies of solutes were solved. These results were similar to those obtained from MD simulations.