PHOTONIC BANDS OF DISPERSIVE PERIODIC STRUCTURES: NUMERICAL IMPLEMENTATION OF AUXILIARY FIELD FORMULATION