A high order local approximation scheme is developed in the context of the numerical manifold method.
The stresses derived from the scheme are continuous at all nodes, causing stress polish unnecessary.
The scheme owns high precision but overcomes the linear dependency issue.