A phenomenological theory of bone remodeling was developed with improved spatial stability. The improved stability was created by changing the nature of the remodeling differential equation to have an exponential character. The new theory was used to predict the distribution of apparent density in the proximal femur. The mesh, applied loads, results and algorithms for the model have been placed in an Internet accessible archive.