The discrete complex image method (DCIM) is applied to compute the complete set of electric, magnetic, and potential Green functions (or integral equation kernels) for planar multilayers with uniaxial anisotropy. It is assumed that the layered medium is of infinite lateral extent and the optic axis is normal to the stratification. Sample results are presented and are shown to be in good agreement with reference data obtained by a rigorous integration method.