n raquo ;, r [i]);
}=0; (int i=0; i lt; n; i ++)
{(max lt; fabs (r [i]))
{= fabs (r [i]);
}
} ( n || Ax-f || =% 1.18f n raquo ;, max);
} Gauss_column (double A [n] [n], double f [n])
{max; line=0, column=0;//позиції максимального елемента
double d, x [n], r [n], Ax [n];
double A_1 [n] [n + 1];//створюємо розширену матрицю
for (int i=0; i lt; n; i ++)
{(int j=0; j lt; n; j ++)
{_ 1 [i] [j]=A [i] [j];
} _ 1 [i] [n]=f [i];
} lt; lt; endl; (int k=0; k lt; n; k ++)
{= fabs (A_1 [k] [k]);=k;=k; (int i=k; i lt; n; i ++)
{(fabs (A_1 [i] [k]) gt; max)
{= fabs (A_1 [i] [k]);
line=i;//знаходження макс елемента і його позиції
}
} (line!=k)//міняємо рядки місцями
{(int j=k; j lt; n + 1; j ++)
{(A_1 [k] [j], A_1 [line] [j]);
}
}=A_1 [k] [k]; (int j=k; j lt; n + 1; j ++)//поділ к-го рядка на макс елемент
{_ 1 [k] [j]=(double) A_1 [k] [j]/d;
} (int i=0; i lt; n; i ++)
{(i!=k)
{= A_1 [i] [k]; (int j=k; j lt; n + 1; j ++)
{_ 1 [i] [j] -=d * A_1 [k] [j];
}
}
}
} lt; lt; endl; lt; lt; Root lt; lt; endl; (int i=0; i lt; n; i ++)
{[i]=A_1 [i] [n]; (% .4f raquo ;, x [i], ); lt; lt; endl;
} (A); lt; lt; endl; (int i=0; i lt; n; i ++)
{s=0; (int j=0; j lt; n; j ++)
{+=A [i] [j] * x [j];
} [i]=s; [i]=Ax [i] -f [i];
} lt; lt; Nevazka lt; lt; endl; (int i=0; i lt; n; i ++)
{(% 1.18f n raquo ;, r [i]);
}=0; (int i=0; i lt; n; i ++)
{(max lt; fabs (r [i]))
{= fabs (r [i]);
}
} ( n || Ax-f || =% 1.18f n raquo ;, max);
} Gauss_Line (double A [n] [n], double f [n])
{P [n]; (int i=0; i lt; n; i ++) [i]=i; max;
int line=0, column=0;//позиції максимального елемента
double d, x [n], r [n], Ax [n];
double A_1 [n] [n + 1];//створюємо розширену матрицю
for (int i=0; i lt; n; i ++)
{(int j=0; j lt; n; j ++)
{_ 1 [i] [j]=A [i] [j];
} _ 1 [i] [n]=f [i];
} lt; lt; endl; (int k=0; k lt; n; k ++)
{= fabs (A_1 [k] [k]);=k;=k; (int j=k; j lt; n; j ++)
{(fabs (A_1 [k] [j]) gt; max)
{= fabs (A_1 [k] [j]);=j;
}
} (column!=k)//міняємо стовпці місцями
{(int i=0; i lt; n; i ++)
{(A_1 [i] [k], A_1 [i] [column]);
} (P [k], P [column]);
}=A_1 [k] [k]; (int j=k; j lt; n + 1; j ++)//поділ к-го рядка на макс елемент
{_ 1 [k] [j]=(double) A_1 [k] [j]/d;
} (int i=0; i lt; n; i ++)
{(i!=k)
{= A_1 [i] [k]; (int j=k; j lt; n + 1; j ++)
{_ 1 [i] [j] -=d * A_1 [k] [j];
}
}
}
} lt; lt; endl; (int i=0; i lt; n; i ++)
{[P [i]]=A_1 [i] [n];
} lt; lt; Root lt; lt; endl; (int i=0; i lt; n; i ++)
{(% .4f raquo ;, x [i], ); lt; lt; endl;
} (A); lt; lt; endl; (int i=0; i lt; n; i ++)
{s=0; (int j=0; j lt; n; j ++)
{+=A [i] [j] * x [j];
} [i]=s; [i]=Ax [i] -f [i];
} lt; lt; Nevazka lt; lt; endl; (int i=0; i lt; n; i ++)
{(% 1.18f n raquo ;, r [i]);
}=0; (int i=0; i lt; n; i ++)
{(max lt; fabs (r [i]))
{= fabs (r [i]);
}
} ( n || Ax-f || =% 1.18f n raquo ;, max);
} Determinant (double A [n] [n])
{max, p=1; s=1;
int line=0, column=0;//позиції максимального елемента
double d; (int k=0; k lt; n; k ++)
{= fabs (A [k] [k]);=k;=k; (int i=k; i lt; n; i ++)
{(int j=k; j lt; n; j ++)
{(fabs (A [i] [j]) gt; max)
{= fabs (A [i] [j]);
line=i;//знаходження макс елемента і фікс позиції=j;
}
}
} (line!=k)//міняємо рядки місцями
{(int j=k; j lt; n; j ++)
{(A [k] [j], A [line] [j]);
} *=- 1;
} (column!=k)//міняємо стовпці місцями
{(int i=k; i lt; n; i ++)
{(A [i] [k], A [i] [column]);
} *=- 1;
}=A [k] [k]; (int j=k; j lt; n; j ++)//поділ к-го рядка на макс елемент
{[k] [j]=(double) A [k] [j]/d;
}=p * d; (int i=k + 1; i lt; n; i ++)
{= A [i] [k]; (int j=k; j lt; n; j ++)
{[i] [j] -=d * A [k] [j];
}
}
} lt; lt; s * p; lt; lt; endl;
} Multipluying (double A [n] [n], double B [n] [n])
{C [n] [n]; s=0; (int i=0; i lt; n; i ++)
{(int j=0; j lt; n; j ++)
{= 0; (int k=0; k lt; n; k ++)
{= s + A [i] [k] * B [k] [j];
} [i] [j]=s;
}
} (C);
...
