  We consider conformal N=2 super Yang-Mills theories with gauge group SU(N)
and Nf=2N fundamental hypermultiplets in presence of a circular 1/2-BPS Wilson
loop. It is natural to conjecture that the matrix model which describes the
expectation value of this system also encodes the one-point functions of chiral
scalar operators in presence of the Wilson loop. We obtain evidence of this
conjecture by successfully comparing, at finite N and at the two-loop order,
the one-point functions computed in field theory with the vacuum expectation
values of the corresponding normal-ordered operators in the matrix model. For
the part of these expressions with transcendentality zeta(3), we also obtain
results in the large-N limit that are exact in the 't Hooft coupling lambda.
