By a proper treatment of the known boundary conditions of a boundary value problem, a complex variable boundary element method (CVBEM) can be used to exactly satisfy the known nodal point boundary values. In this fashion, a numerical model can be developed which generates relative error information along the problem boundary that can be used to reduce the modelling error by either an integrated measure or a maximum relative error measure.