The main objective of this paper is to determine the aerodynamic coefficients of a steady supersonic wing, using boundary element method based on the solution of potential flow. In this regard, the governing partial deferential equation of flow field is first derived and linearized. Then, using boundary element method, in conjunction with the first order shape functions, the system of discrete equations are obtained. However, there are some numerical problems in the evaluation of integral coefficients due to singularity caused by intersection of Mach forecone and the edge of some elements in computational domain. Of course, this problem can be eliminated using some numerical techniques. Finally, the steady supersonic flow analysis is performed for two different wings with identified geometries. Comparison between the numerical results and those obtained using the exact solution indicates good agreements.