In order to achieve uniform distribution of heat flux over an isothermal heated target surface, a numerical algorithm is developed to obtain an optimized array of four laminar impinging slot jets. Root mean square deviation of the local Nusselt distribution from the desired Nusselt number is considered as the objective function. Jets' widths, jet-to-jet and jet-to-surface spacings, and the overall flow rate are chosen as design variables. Conjugate gradients method along with backtracking line search is applied to optimize the objective function calculated by numerical simulation for three different cases of Nu = 7, 9, and 11. For each of these desired Nusselt numbers, an almost uniform distribution of local Nusselt number with percent of root mean square error less than 2.5% is achieved in fewer than 12 iterations. An experimental study using a Mach–Zehnder interferometer (MZI) has been performed. The measured distribution of local Nusselt number is in good agreement with numerical results in all three optimal configurations.