A Gibbs free energy method improved via initiating initial values by computer with random direction method was presented to calculate the equilibrium composition of combustion products. The factors of solving initial number and interval selection affecting the random direction method were discussed. Combustion products of hydrazine and oxygen obtained by improved and original Gibbs free energy methods are same. The advances of improved arithmetic are simple, convergent and initial value matched a constraint condition completely which is generated by computer.