Neumann boundary value problems are solved for stellarator fields generated by external currents. The method can be applied to the construction of external coil configurations of stellarators and to studies of the properties of vacuum fields. A computer code, NESCOIL, has been developed and applied to the Helias and ATF stellarator configurations.