Чисельне розв’язання системи лінійних алгебраїчних рівнянь за методом простих ітерацій
17
printf("\t\t%4.1f",b[i]);
printf("\n");
}
printf("\n");
do
{
Далі іде ітераційний процес:
norma=0.0;
for(i=0;i < N;i++)
{
xn[i]=-b[i];
for(j=0;j < N;j++)
{
if(i!=j)
xn[i]+=a[i][j]*x[j];
}
xn[i]/=-a[i][i];
}
for(i=0;i < N;i++)
{
Далі іде обчислення норми вектора:
if(fabs(x[i]-xn[i]) >norma)
norma=fabs(x[i]-xn[i]);