# using System using System Collections Generic using System Linq using

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122``` ```using System; using System.Collections.Generic; using System.Linq; using System.Text; namespace L1 { class Seidel { Matrix mMatrix; Matrix b; int numIter = 0; public Seidel(Matrix matrix, Matrix b) { this.mMatrix = matrix; this.b = b; } public double[] getSolve(double eps) { Matrix R = getMatrixR(mMatrix), L = getMatrixL(mMatrix), D = getMatrixD(mMatrix); Matrix Q1 = -1 * (D + L).GetInverse() * R; R = getMatrixR(mMatrix); L = getMatrixL(mMatrix); D = getMatrixD(mMatrix); Matrix Q2 = (D + L).GetInverse() * b; Matrix curX = new Matrix(b), nextX = new Matrix(mMatrix.GetLength(0), 1); while (Math.Abs((nextX - curX).GetNorm()) >= eps) { ++numIter; if (numIter > 100) { throw new Exception("Решение не сходится."); } curX = nextX; nextX = Q1 * curX + Q2; } return nextX.GetColumnArray(0, 0); } private Matrix getMatrixL(Matrix original) { Matrix matrix = new Matrix(original); for (int i = 0; i < matrix.GetLength(0); ++i) { for (int j = 0; j < matrix.GetLength(1); ++j) { if (i <= j) { matrix[i, j] = 0; } } } return matrix; } private Matrix getMatrixR(Matrix original) { Matrix matrix = new Matrix(original); for (int i = 0; i < matrix.GetLength(0); ++i) { for (int j = 0; j < matrix.GetLength(1); ++j) { if (i >= j) { matrix[i, j] = 0; } } } return matrix; } private Matrix getMatrixD(Matrix original) { Matrix matrix = new Matrix(original); for (int i = 0; i < matrix.GetLength(0); ++i) { for (int j = 0; j < matrix.GetLength(1); ++j) { if (i != j) { matrix[i, j] = 0; } } } return matrix; } public int getNumIter() { return numIter; } public bool isValid() { for (int i = 0; i < mMatrix.GetLength(0); ++i) { double maxValue = -100; int ind = 0; for (int j = 0; j < mMatrix.GetLength(1); ++j) { if (mMatrix[i, j] > maxValue) { maxValue = mMatrix[i, j]; ind = j; } } if (ind != i) return false; } return true; } } } ```