09/07/06 14:25:27
>>624 ガウス・ザイデル法の計算方法 山本昌志 2005年12月16日
URLﾘﾝｸ(akita-nct.jp)
必要なら以下にコメントも付けとく。仮にｴｸｾﾙの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