The electronic states of ZnO doped with Al, Ga and In, which belong to III family elements in periodic table, were calculated using the density functional theory. In this study, the calculation was performed by two Programs; the discrete variational Xa (DV-Xa) method, which is a sort of molecular orbital full potential method; Vienna Ab-initio Simulation Package (VASP), which is a sort of pseudo potential method. The fundamental mixed orbital structure in each energy level near the Fermi level was investigated with simple model using DV-Xa. The optimized crystal structures calculated by VASP were compared to the measured structures. The density of state and the energy levels of dopant elements were shown and discussed in association with properties.