Integral equation formulation for natural modes of a circular patch antenna in a layered environment