public static Matrix Integrate(Matrix A, comp a, int toRow, comp b, int fromRow)
{
if (toRow >= 0 && fromRow >= 0 && toRow < A.m && fromRow < A.m)
{
comp[,] bArray = new comp[A.m, A.n];
for (int j = 0; j < A.n; j++)
{
A.SetValue(a * A.GetValue(toRow, j) + b * A.GetValue(fromRow, j), toRow, j);
}
return(A);
}
return(new Matrix(new comp[0, 0]));
}
public static Matrix Exchange(Matrix A, int to, int from)
{
if (to >= 0 && from >= 0 && to < A.m && from < A.m)
{
comp[,] b = new comp[A.m, A.n];
for (int j = 0; j < A.n; j++)
{
b[to, j] = A.GetValue(from, j);
}
Matrix B = new Matrix(b);
return(A + B);
}
return(new Matrix(new comp[0, 0]));
}
public static Matrix operator *(Matrix A, Matrix B)
{
if (A.n == B.m)
{
comp[,] AB = new comp[A.m, B.n];
for (int i = 0; i < B.n; i++)
{
for (int j = 0; j < A.m; j++)
{
for (int k = 0; k < A.n; k++)
{
AB[i, j] += A.GetValue(i, k) * B.GetValue(k, j);
}
}
}
return(new Matrix(AB));
}
return(new Matrix(new comp[0, 0]));
}
public void SetValue(comp value, int i = 0, int j = 0)
{
x[i, j] = value;
}
