导数计算

今天发现有限元插件不稳定的原因在于用有限差商计算导数的数值误差。

之前的步长都是1e-4，现在步长变为1e-2（0.01）就工作正常了。

下一步要找一种健壮的、同时计算量不太大的导数计算方法。