Аннотация:We consider an application of the barycentric method for the numerical solution of boundary value problems of mathematical physics. The main assumptions are that the system of partial differential equations of the boundary value problem is solvable in the approximation of the Galerkin method and the boundary of the domain of analysis is piecewise linear. A distinctive feature of the barycentric method is the formation of a global system of basis functions for the area of analysis via barycentric coordinates. Solutions for determining barycentric coordinates are given. A comparison is made of the rate of convergence of the barycentric method and grid methods when solving some typical boundary value problems of mathematical physics.