AbstractSome new formulae are presented for numerical integration over triangular regions. These formulae, which are triangularly symmetric, involve four, six, and seven integration points, respectively. They were derived by determining co‐ordinates and weight factors for which polynomial functions of up to a certain degree are integrated exactly. The effectiveness of the formulae is demonstrated by applying them to the integration of five non‐polynomial functions, two of which are singular within the region of integration. Results obtained with the new formulae are compared with results obtained with other, comparable formulae. It is seen that the new formulae are more accurate than some existing ones with the same numbers of poi