Multi-stage distillation column for the concentration of nitric acid was completely analyzed. Composition, hold up and enthalpy changes （CHE） are considered in the MESH equations and thus for binary systems, the steady states can be expressed as the solutions of 3(N+1) nonlinear algebraic equations for column of N stages and a reboiler stage. Here, solving the hold up and enthalpy balance equations analytically the liquid and vapor flow rates are obtained as the functions of concentration of each stage and the number of equations reduced to (N+1). For this system, Newton-Raphson method and implicit function theorem are applied to obtain the steady state behavior depending on the change of opeartion parameters such as feed stage, feed concentration, temperature, reflux flow rate and reboiler heating rate, and consequently the optimal operating conditions are suggested.