Esta investigación trata sobre cálculos de precisión de campos magnéticos estacionarios de cuerpos volumétricos. La analogía electrostática permite el uso de un potencial magnético escalar, lo que reformula la tarea original como un problema de valores límite para la ecuación de Laplace. Abordamos esto con el método de elementos de contorno, específicamente en rangos de distancia cercanos a la superficie magnetizada, donde se sabe que los métodos numéricos estándar existentes tienen dificultades. Este trabajo presenta un enfoque basado en fórmulas de cuadratura mejoradas para el potencial de capa simple y su derivada normal. Pruebas numéricas confirman mejoras significativas en el cálculo del campo a cualquier distancia desde la superficie del imán.