scilab

scilab を使うと微分方程式がいとも簡単に解けるなんて知らなかった。
y'=x-2y y(0)=1　を解くにはこれだけでいいようだ。
deff("ydot=df(x,y)", "ydot=x-2*y"); //方程式の定義
x0=0 ; y0=1 ; //初期条件
x=0:0.1:5; //解の範囲
y=ode(y0 , x0 , x , df ); //解
plot(x,y) //プロット
１階微分方程式はすべてこれでできそう。かつてFORTRANで入出力で何回も失敗したし，CやJavaではグラフの描き方がまだよくわからん･･･のにくらべるとこれでいいのかなと思ってしまいます。