Gauss列主元法(C#实现)——计算方法
Gauss列主元法是在高斯消去法的基础上增加找主元与换行的操作,因此为了不使代码冗余,可用面向对象的继承的特点
代码实现如下:
1.定义Guass1(Gauss列主元法类)
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
namespace Line
{
class Guass1:Gauss
{
public void FindMaster(int k, ref double[,] a, int n)//找主元 换行
{
bool flag = false;
double s = a[k, k];
int m = k;
for (int i = k + 1; i <= n - 1; i++)
{
if (Math.Abs(s) < Math.Abs(a[i, k]))
{
m = i;
s = a[i, k];
flag = true;
}
}
if (flag == true)
{
for (int v = k; v <= n; v++)
{
double t = a[k, v];
a[k, v] = a[m, v];
a[m, v] = t;
}
}
}
public new void Elimination(int n, ref double[,] a)//消元
{
for (int k = 0; k <= n - 2; k++)
{
FindMaster(k, ref a, n);
for (int i = k + 1; i <= n - 1; i++)
{