The implementation of the Natural Element Method (NEM) by means of a Galerkin-based procedure is parallel to that adopted in FEM or Element-Free Galerkin (EFG) method [14], with the key distinction that separates the three is in the construction of the shape functions [画像:$ \phi_I(\mathbf{x})$] and their derivatives.
A computational procedure to evaluate the shape functions (n-n coordinates) is outlined in [3], which is extended by [9] to compute the derivatives of the interpolating function [画像:$ u^h(\mathbf{x})$]. In the application of n-n interpolation to PDEs (fluid-structure interaction), [10] adopted Lasserre's recursive formula [15] to compute the area (volume of a convex polyhedral in $ \mathbf{R}^k$) of the second-order Voronoi cells. The shape functions [画像:$ \phi_I(\mathbf{x})$] as well as their derivatives [画像:$ \phi_{I,j}(\mathbf{x})$] are computed. In [10], it is pointed out that Lasserre's formula [15] is more robust than the one due to [3], which breaks-down if the point [画像:$ X(\mathbf{x})$] lies along the edge of a Delaunay triangle. Expressions for the shape functions and their derivatives are provided [3,16,9,10].
