A simple theoretical method for calculating electronic band structures of wurtzite materials based on the bond-orbital models is presented. This method can be used to study many problems such as band mixing and effects of external fields (electric field, magnetic field, and unaxial stress, etc.), since it can reproduce fairly accurate lowest conduction-band and top three valence-band structures. This method is very similar to LCAO method; however, it is much simpler and requires less computational effort than LCAO method.