The problem of systematically approximate distributions
numerically will be addressed. We shall present a general ``physical
space'' method based on the combined choice of an appropriate smooth
orthonormal basis and an associated quadrature rule. Various applications
will be mentioned, including interpolation, resolution of the Gibbs
phenomenon and PDEs.