11/11/08 00:40:01.34
URLﾘﾝｸ(qanda.rakuten.ne.jp)にある質問なのですが
これをFORTRANでできる方はいらっしゃいますか？
あとこのプログラミングをC++でやってみたところ
const int N = 100;
const double q = 10.0, dt = 0.00001, Dm = 30.0, t0 = 2.0,
K = 1.0, pi = 3.141592, f = 3.0;
double C[N], dC[N];
double dx = q/N;
for (int i = 0; i < N; ++i) C[i] = 0; // 初期条件
for (double t = 0; t < t0; t += dt) {
C[0] = C[1]; // 境界条件1
C[N - 1] = K*sin(2*pi*f*t); // 境界条件2
for (int i = 1; i < N - 1; ++i) dC[i] = (Dm*(C[i + 1] - 2*C[i] + C[i - 1])/(dx*dx))*dt;
for (int i = 1; i < N - 1; ++i) C[i] += dC[i];
double s = 0;
double s_prev;
for (int i = 1; i < N - 1; ++i) s += C[i]*dx;
if (t) {
double ss = (s - s_prev)/dt;
cout << ss << endl; // ∂s/∂t を出力
s_prev = s;
}
といちおやってみたのですが、間違いを指摘できる方はいますか？
あとできればこのプログラムをFORTRANに変換もしていただきたいです。。。
重複と長文同時に失礼いたしました。
図はURLﾘﾝｸ(qanda.rakuten.ne.jp)にあります。
全開質問に答えていただいた方、本当にありがとうございました。