星空の夢

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 となるから実は解なしです。

京都にちなんだ小惑星4

　平家物語の登場人物も小惑星になっているものが多数あります。源氏では(5744)頼政，(3733)義朝，(3902)頼朝，(3178)義経，(4574)義仲，女性には(4748)常盤御前，(4200)静御前，(4896)巴御前･･･まだまだあるかもしれません。平家ではもっとあります。(4374)忠盛，(4375)清盛，(4376)重盛，(4377)惟盛，(4488)時忠，女性には(4959)二位の尼，(5242)建礼門院・・・
もちろん天文学への関係はありませんが。

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