精确计算形状函数高阶导数的插值/近似方法

你好,

我想用无网格法解决固体力学中的一些问题。你知道一些插值/近似方法，能够准确地计算形状函数的三阶和四阶导数，并且不难实现吗?众所周知，MLS或PIM在这些衍生品的准确性方面存在问题。
非常感谢你的建议。