C/C++の宿題片付けます 128代目

C/C++の宿題片付けます 128代目 at TECH

677: ◆/91kCCQXBo
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

次ページ

続きを表示

1を表示