- 练习17可知,都在稳定域时,rk 方法精度很高,但隐式欧拉似乎不会崩溃,我记得稳定域是整个平面来着(几年前上的课了),rk 在 lambada 到-500 多就崩溃了
- 隐式欧拉的编写让我学会了 lambda 函数的写法,ode 类中的函数是
f(double t, doublt x);
fder(double t, double x); // f() 对 x 的偏导,newton 方法迭代时要用
而我以前写的 newton() 是单参数的
newton(double x);
而且迭代时实际求的函数是 其解正是$x_{i+1}$,编写过程中使用 lambda 以方便地生成迭代需要的函数及其导数,而且返回类型需是
funcion<double(double)>
其在<functional>
中定义,newton()
也需用模板,这些都是我第一次接触的
- 利用
gnuplot
和powershell
完成批量做图,在windows
下用胶水的话,感觉powershell
比python
更好