The equifrequency contours of a two-dimensional square lattice, Photonic Crystal(PhC) composed of an anisotropic
cylinder in anisotropic medium are modeled and theoretically studied. The beam group velocity inside the PhC is
determined as a function of input angle and from this data a quantitative measure of the self-collimation effect is
deduced. The range of input angles for the self-collimation regime are significantly increased above those for a PhC
composed of isotropic constituents.