The Boundary integral method (BIM) is applied to a number of planar microstrip multiport junctions. The E fields for coupling ports are expressed as multimode series, enabling direct calculation of S parameters from the resulting matrix equation. Results obtained for cross junctions, T junctions, Y junctions, and 90° bends agree very well with those obtained from full‐wave analysis and other methods. © 1994 John Wiley & Sons, Inc.