Zero dimensional climate modeling is one method to simplify describe the correlation among the influence elements. This model is consisting of system of two mathematical equations. One of the models was proposed by Boeker-Van Grondelle. It explained the overall description of the parameters that involved on the process occurred between the atmospheric layer and surface of the earth. Solution for the model could be obtained by numerical methods. On this work, Quasi Newton-Raphson was used to solve model proposed by Boeker-Van Grondelle. 0.5 of relative error was used as the stopping criteria. The results showed that this method was able to satisfying the stopping criterion in three steps iteration and the established solution. © 2011 IEEE.