|
", k);
for(i=0;i<n;i++){
sum=0.;
for(j=0;j<n;j++){
if(i!=j)
sum=sum+A[i][j]*x[j];
}
next_x[i]=(A[i][n]-sum)/A[i][i];
printf("%lf \t", next_x[i]);
x[i]=next_x[i];
}
err = 0.;
sum = 0.;
for(i=0;i<n;i++){
for(j=0;j<n;j++)
sum = sum + pow(A[i][j],2);
}
err = pow(sum,0.5);
if(err <= esp) /
|
|
메뉴펼치기
