A method based on properties of simplex co-ordinates is proposed to be used in dual finite-element calculations in magnetostatics. The method considers the source's distribution within the field region and the interelement boundary conditions. The provision of upper and lower bounds, in connection with the obtained high rate of convergence, makes the method extremeley powerful and flexible. The simplicity and applicability of the method is shown by an example, where some previous results are also discussed.