The structures of crystalline vitreous and liquid $SiO_2$ were Monte carlo simulated employing the potential energy function comprising Lennard-Jones 2-body and Axilrod-Teller 3-body potentials. Although the Si-O-Si angular distribution functions obtained in the simulation appear to be higher than the experimental results the other simulation results including SiO, O-O and Si-Si radial distribution functions and O-Si-O anglular distribution functions agree well with experimental data within acceptable limits. Themost important outcome in this study is that various $SiO_2$forms were successfully reproduced with the same potential energy function.