A new numerical approach has been developed for vapor solid equilibrium calculations and for predicting vapor solid equilibrium constant and composition of vapor and solid phases in gas hydrate formation. Equation of state methods generally do a good job of determining vapor phase properties,but for solid phase it is much more difficult and inaccurate. This proposed new model calculates vapor solid equilibrium constant and vapor and solid phase composition as a function of temperature and partial pressure. The results of this proposed numerical approach, for vapor solid equilibrium, have a good agreement with the available reported data. This new numerical model also has an advantage to tune coefficients, to cover different sets of experimental data accurately.