In this paper, the network equation for the slow neutron capture process （s-process） of heavy element nucleosynthesis is investigated. Dividing the s-process network reaction chains into two standard forms and using the technique of matrix decomposition, a group of analytical solutions for the network equation are obtained. With the analytical solutions, a calculation for heavy element abundance of the solar system is carried out and the results are in good agreement with the astrophysical measurements.
