We propose a discretization scheme for a numerical solution of elliptic PDE
's, based on local representation of functions, by their Taylor polynomials
(jets). This scheme utilizes jet calculus to provide a very high order of
accuracy for a relatively small number of unknowns involved.