gnuplotで

eq1

 方程式が解けてしまった。gnuplotはもっぱら関数のグラフを描くツールと思っていたが、実はややこしい非線型方程式が解けます。
図は小さくて見にくいが sin(-x)-x-4=0 の解が x=-4.96762762452019と求まったというわけ。ちょっとだけニュートン法によるくり返し文を書いていますが。多分ほとんどの方程式は解けそうです。ただしSclabやOctoveの便利な関数は効かないみたいですが。→ここはAWさんにお聞きしたい
 さらに積分などもできるでしょう、これから考えてみよう。

comments