Matlab code for tri-resonator structure

I'm working on tri-resonator structure and want to generate coupling matrix of 3x3 including externa quality factors qe1, qe2=qe3.

Can any one help? many thanks