Hybrid Numerical/Asymptotic approach for the Analysis of Rectangular Array Greens Function in a Multilayered Dielectric Slab