Equilibrium Solutions
>
equil:= solve( {rhs(rate_eqn1), rhs(rate_eqn2)}, vars );
>