  A study of the quantum triangular billiard requires consideration of a
boundary value problem for the Green's function of the Laplacian on a trianglar
domain. Our main result is a reformulation of this problem in terms of coupled
non--singular integral equations. A non--singular formulation, via Fredholm's
theory, guarantees uniqueness and provides a mathematically firm foundation for
both numerical and analytic studies. We compare and contrast our reformulation,
based on the exact solution for the wedge, with the standard singular integral
equations using numerical discretisation techniques. We consider in detail the
(integrable) equilateral triangle and the Pythagorean 3-4-5 triangle. Our
non--singular formulation produces results which are well behaved
mathematically. In contrast, while resolving the eigenvalues very well, the
standard approach displays various behaviours demonstrating the need for some
sort of ``renormalisation''. The non-singular formulation provides a
mathematically firm basis for the generation and analysis of periodic orbit
expansions. We discuss their convergence paying particular emphasis to the
computational effort required in comparision with Einstein--Brillouin--Keller
quantisation and the standard discretisation, which is analogous to the method
of Bogomolny. We also discuss the generalisation of our technique to smooth,
chaotic billiards.
