Magnetostatic permeability tomography is an imaging technique that attempts to reconstruct the permeability distribution of an object using magnetostatic measurement data. The data for image reconstruction are external magnetic field measurements on the surface of the object due to an applied magnetostatic field. Theoretically, the normal and tangential components of the magnetic field in the surface uniquely define the internal isotropic permeability distributions. However, the inverse permeability problem is an ill-posed nonlinear problem. Regularization is needed for a stable solution. In this paper, we present a numerical method to solve the reconstruction problem in three dimensions using a regularized Gauss-Newton scheme. We have solved the forward problem using an edge finite-element method and we have employed an efficient technique to calculate the Jacobian matrix. The permeability of the object is assumed to be linear and isotropic. We present the reconstruction results for the permeability using synthetically generated data with additive noise.