When a high-frequency acoustic or electromagnetic wave is scattered by a surface with a conical point, the component of the asymptotics of the scattered wave corresponding to diffractio by the conical point can be represented as an asymptotic expansion, valid as the wave number k.8. The diffraction coefficien is the coefficien of the principal term in this expansion and is of fundamental interest in high-frequency scattering. It can be computed by solving a family of homogeneous boundary value problems for the Laplace-Beltrami-Helmholtz equation ( parametrized by a complex wave number-like parameter.) on a portion of the unit sphere bounded by a simple closed contour , and then integrating the resulting solutions with respect to.. In this paper we give the numerical analysis of a method for carrying out this computation ( in the case of acoustic waves) via the boundary integral method applied on , emphasizing the practically important case when the conical scatterer has lateral edges. The theory depends on an analysis of the integral equation on , which shows its relation to the corresponding integral equation for the planar Helmholtz equation. This allows us to prove optimal convergence for piecewise polynomial collocation methods of arbitrary order. We also discuss efficien quadrature techniques for assembling the boundary element matrices. We illustrate the theory with computations on the classical canonical open problem of a trihedral cone.