- 677 名前： ◆/91kCCQXBo mailto:sage [2009/07/06(月) 14:25:27 ]
- >>624 ガウス・ザイデル法の計算方法 山本昌志 2005年12月16日
akita-nct.jp/yamamoto/lecture/2005/5E/Linear_eauations/ex_GaussSeidel_html/node2.html
必要なら以下にコメントも付けとく。仮にｴｸｾﾙのPADでもいい？
#if 0
for( j=1; j<=N; j++ ){
temp += a[i][j] * x[j];
}
/* a[i][i] * newx + (temp - a[i][i] * x[i]) = b[i]; */
/* a[i][i] * newx = b[i] - (temp - a[i][i] * x[i]); */
newx = 1.0 / a[i][i] * (b[i] - (temp - a[i][i] * x[i]));
#else
for( j=1; j<=N; j++ ){
if( i==j ) continue;
temp += a[i][j] * x[j];
}
/* Xi(k+1) = 1/Aii * ( Bi - (Ai,1*X1(k) + Ai,2*X2(k) + ... + Ai,i-1*Xi-1(k) + Ai,i+1*Xi+1 + ... + Ai,n*Xn )) */
newx = (b[i] - temp)/a[i][i];
#endif
|
 |
