An integral algebraic linear programming is defined and an algorithm for solving the programming is given.

The coupled boundary integral equations and their reduced forms are deduced for hybrid mode problems, and these equations are discreted into a system of linear algebraic equations by moment method.

Gauss Integral method is populated due to its high algebraic accuracy in index calculation.