Abstract: The free energy minimization technique is adopted for calculation of theoretical rocket performance under equilibrium flow conditions. According to optimization theory, the direct search method is employed in combination with a computer program. A composite solid propellant ingredient is chosen for maximum specific impulse and maximum characteristic velocity. The numerical results obtained agree fairly well with the results obtained by graphical method for determining the optimum composition of solid propellant systems.