본문내용
", k);
for(i=0;i
sum=0.;
for(j=0;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
for(j=0;j
sum = sum + pow(A[i][j],2);
}
err = pow(sum,0.5);
if(err <= esp) // 구한 오차가 허용 오차보다 적으면 실행을 중단한다.
{
printf("err=%lf \n", err);
return;
}
}
return;
}
*** Jacobi보다 빨리 수렴함을 확인할 수 있다.
for(i=0;i
for(j=0;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
}
err = pow(sum,0.5);
if(err <= esp) // 구한 오차가 허용 오차보다 적으면 실행을 중단한다.
{
printf("err=%lf \n", err);
return;
}
}
return;
}
*** Jacobi보다 빨리 수렴함을 확인할 수 있다.
소개글