Thank you for your reply, mfb.
I appreciate the point that one has to square afterwards.
I think that my problem lies in the evaluation of the expectation value of the charge operator for each meson. What I'm doing for the ρ meson, for instance, is
\langle u \overline{u}-d...