10楼: Originally posted by 落日照秋草 at 2014-08-13 16:46:01
就是把二阶、一阶微分转换为求差值。
例如：a''(x)=a(x+1)+a(x)-2a(x)/步长的平方
a'(x)=a(x+1)-a(x-1)/2
具体的翻翻数学物理方法有关差值那部分吧。代码看看徐世良的fortran常用算法集。

10楼: Originally posted by 落日照秋草 at 2014-08-13 16:46:01
就是把二阶、一阶微分转换为求差值。
例如：a''(x)=a(x+1)+a(x)-2a(x)/步长的平方
a'(x)=a(x+1)-a(x-1)/2
具体的翻翻数学物理方法有关差值那部分吧。代码看看徐世良的fortran常用算法集。