星空の夢

scilab2

scilabで非線形方程式も簡単に解けます。
exp(x)-x=2　を解くとき，まずグラフを描くには
　x=-5:0.05:2;y=exp(x)-x-2;plot(x,y);plot(x,0)
グラフを見て初期値を決める。x0=1とすれば
　deff('y=f(x)','y=exp(x)-x-2');x0=1;fsolve(x0,f)　より
　1.1461932　が得られる。
また−１とすれば - 1.8414057が得られる。
初期値を 0 とすると解は 0 と表示されるが，この場合 微分係数が 0 となるから実は解なしです。

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ではグラフの描き方がまだよくわからん･･･のにくらべるとこれでいいのかなと思ってしまいます。