2009年10月26日月曜日へのコメント

return 0は関数の実行結果の値をゼロにするという意味だと思います
main関数でこうしてもあまり意味がないように思います

stdio.hは標準の入出力で利用する関数を呼び出すために必要な
情報を与えていると思えば良いと思います

来週の予告

磁場中の荷電粒子の運動（ローレンツ力あり）

運動方程式（xとyの連立方程式）
m*(dVx/dt)= q*Vy*B
m*(dVy/dt)=-q*Vx*B

どうやって解くか考えましょう。

重力場中（空気抵抗あり）の物体のx-tグラフのための数値計算

#include < stdio.h >
main()
{
double t, x, dt, dx, ta, tb, v;
double f(double ,double );
int i, n;
x = 0.0 ;
t = 0.0 ;
ta = 0.0 ;
tb = 100 ;
v = 0.0 ;
n = 1000;
dt = (tb - ta) / n ;
for ( i = 0; i < n ; i++){
x = x+ v*dt ;
v = v+ f(v, t)*dt ;
t = t + dt ;
printf(" %e %e \n", t, x);
}
}
double f( double x, double t){
double y, z, alpha, g;
alpha = 0.1;
g = - 9.8;
z = -alpha * x + g ;
return z;
}

来週のメニュー

来週は
2階の微分方程式（1次の項を含む）を
電解中の荷電粒子の運動（空気抵抗あり）をモデルに
やります。

ヒント
for()
{
x = x + v*dt
v = v + (A-B*v)*dt

}

(A=qE/m, B=alpha/m)

グラフを画像ファイルとして保存する

gnuplotで表示した画像の保存
clipboardに保存（グラフ画面の左上のアイコンをクリック）
GIMPを起動しファイルメニューから”画像の生成”を選び”クリップボードから”をクリック。するとグラフが表示される。名前をつけて保存を選び、形式を”JEPG画像”にして保存。

ブログにアップロードするには、ブログの投稿画面で”画像を追加”から”参照”

計算結果を保存する

コンパイルしたプログラムを
./a.out
で端末に計算結果が表示されるのでそれを
./a.out > ファイル名
で計算結果をファイルに保存する。