We present an efficient adaptive method which can be used for the numerical solution of various problems of physics and engineering described by partial differential equations.
