A general interpolation/fitting method is presented for constructing a smooth molecular potential energy surface from ab initio energy points. The approach is based on generalized discrete variable representation theory, where a damping term has been introduced. The method is illustrated by constructing a one-dimensional potential representing H-2 and a two-dimensional surface representing collinear H-3.