Numerical methods for the inverse problem of density functional theory