A method of calculating field functions of circular coils by direct numerical integration is explained, and formulae for the integrands required are given. The method is not only more straightforward mathematically than classical methods, but is also faster. A known principle has been extended, enabling a comprehensive treatment to be presented. The method can be used for calculating both magnetic fields of current-carrying coils, and induced EMFS of sense coils: a complete set of the formulae involved is given. Ways of extending the method to the cases of sets of coils and finite coil volumes and finite field volumes are explained.