In silico modeling of Gene Regulatory Networks has gained a lot of attention recently as it gives a very powerful tool to experimental biologists to gather the knowledge gained from different biological experiments and understand the dynamics of the overall system. One of the key dynamics that is often interesting is the steady states of the networks which biologically corresponds to the cellular states. In our previous paper, we gave an efficient method called GenYsis to compute these steady states in Boolean representation of Gene Regulatory Network. It has been observed that protein may be expressed at more then two level of expression. This may result in different cellular outcomes. To address this issue, we present here a multiple-level modeling methodology that allows us to be more accurate. In this paper we extend our software GenYsis to model gene regulatory networks where each node in the network may take multiple values.