In the present paper we propose a new method for constructing a second orde
r Moving Least Squares (MLS) approximation. The method leads to shape funct
ions which are then used for solving Partial Differential Equations (PDE) b
y a collocation method. This work is an extension of the Generalized Finite
Difference Method originally proposed by Liszka and Orkisz (GFDM). However
it differs from GFDM by using a sequence of two first order numerical deri
vations based on linear polynomial basis instead of a second order derivati
on based on a quadratic polynomial basis. This two-stage approach leads to
continuous approximation coefficients using a limited number of surrounding
points and results into quite a simple program structure, very similar to
that of the finite elements. The method is in an early stage of development
so no definitive conclusions may be drawn, however example problems exhibi
t good convergence properties.