全主元高斯消元法

#include 
#include 
#define M 1000
void swap（double &adouble &b）         ///行列交换
{
    double c;
    c=a;
    a=b;
    b=c;
}

int main（）
{
    int n;
    printf（“请输入维数n:“）;
    scanf（“%d“&n）;
    for（;;）
    {
        if（n>=2&&n<=1000）
        {
            break;
        }
        printf（“维数有误，请重新输入:“）;
        scanf（“%d“&n）;
    }
    static double a[M][M]b[M];
    int ij;
    printf（“请输入各个参数:\n“）;
    for（i=0;i    {
        for（j=0;j        {
            printf（“a[%d][%d]=“ij）;
            scanf（“%lf“&a[i][j]）;
        }
        printf（“b[%d]=“i）;
        scanf（“%lf“&b[i]）;
    }
    printf（“增广矩阵为\n“）;
    for（i=0;i    {
        for（j=0;j        {
            printf（“%5.2f\t“a[i][j]）;