Exact superpotentials for theories with flavors via a matrix integral

We extend and test the method of Dijkgraaf and Vafa for computing the superpotential of N=1 theories to include flavors in the fundamental representation of the gauge group. This amounts to computing the contribution to the superpotential from surfaces with one boundary in the matrix integral. We compute exactly the effective superpotential for the case of the gauge group U(Nc), Nf massive flavor chiral multiplets in the fundamental and one massive chiral multiplet in the adjoint, together with a Yukawa coupling. We compare up to sixth order with the result obtained by standard field theory techniques in the already nontrivial case of Nc=2 and Nf=1. The agreement is perfect