FOSPACK is a package for the discretization, assembly, and solution of systems of PDEs using a First-Order System Least-Squares (FOSLS) formulation. The talk will briefly cover the FOSLS motivation and methodology, the solution of the linear systems using Algebraic Multigrid (AMG), as well as a number of features of the package, including adaptive mesh refinement and automatic linearization of nonlinear equations.