  The (classical) Lagrange spectrum is a closed subset of the positive real
numbers defined in terms of diophantine approximation. Its structure is quite
involved. This article describes a polynomial time algorithm to approximate it
in Hausdorff distance. It also extends to approximate the Markov spectrum
related to infimum of binary quadratic forms.
