Демидович Б.П., Марон І.А., «Основи обчислювальної математики», М., Наука, 1963 р., с.660
. С.А.Калоеров, «Введення в програмування на мові С + +»,: Учеб. Посібник, 1999 р., с.183
. В.М.Долгіх, «Лінійна алгебра та аналітична геометрія», 2008р., С.103
Дотаток А
Рішення сістемма лінійніх рівнянь методом Гауса-Жордана мовою Borland C + +4.5.
Дана система:
# include
# include
# include
# include ()
{(); a [11] [11], b [11], x [11], t; i, j, k, m, n;=3; [1] [1]=2; a [1] [2]=1; a [1] [3]=- 1; b [1]=2; [2] [1]=3; a [2] [2]=1; a [2] [3]=- 2; b [2]=3; [3] [1]=1; a [3] [2]=0; a [3] [3]=1; b [3]=3; << « N n»; (i=1; i <= n; i + +)
{(j=1; j <= n; j + +) << a [i] [j] << « »; << b [i] << « N»;
} eps=.000001; max; max_i; lead, a_div_lead; (k=1; k <= n; k + +)
{= 0; _i=- 1; (int i=k; i <= n; i + +)
{(fabs (a [i] [k])> max)
{= fabs (a [i] [k]); _i=i;
}
} (max_i == - 1 | | fabs (a [max_i] [k])
{<< « N n» << «Error: det=0» << « N n»;;
}=a [k] [k]; (j=k; j <= n; j + +) [k] [j] /=lead; [k] /=lead; (i=1; i <= n; i + +)
{_div_lead=a [i] [k] / a [k] [k]; (i!=k)
{(j=k; j <= n; j + +) [i] [j] -=a [k] [j] * a_div_lead; [i] -=b [k] * a_div_lead;
}
}
} << « N n»; (i=1; i <= n; i + +)
{(j=1; j <= n; j + +) << a [i] [j] << « »; << b [i] << « N»;
} ();
}
