Проверьте кому не трудно код, делал вроде правильно, а все равно не так

namespace ConsoleApplication1
{
static double[,] mas =
 {
 {4.3, 0.217, 0, 0, 2.663},
 {0.1, -3.4, -0.207, 0, 2.778},
 {0, 0.09, 2.5, 0.197, 2.533},
 {0, 0, 0.08, -1.6, 1.928}
 };
 
    class CLAY
    {
      static double[,] TreygMas = new double[4, 4];
      static double[] b = new double[4];
      public static int x1 = 0, x2 = 0, x3 = 0, x4 = 0, eps = Math.Pow(10, -6);
 
      static void Main(string args)
      {
       Console.WriteLine("Массив коэфициентов СЛАУ");
       for (int i = 0; i < 4; i++)
       for (int j = 0; j < 5; j++)
      {
       Console.Write(mas[i, j]+"\t");
    }
     Console.WriteLine(); 
     Gayss_Zeidel(mas);
     Console.ReadKey();
    }
      public static void Gayss_Zeidel()
    {
       //n = 0;
      double x01 = 0, x02 = 0, x03 = 0, x04 = 0;
      int sig;
      do
      x1 = (1 / mas[0, 0]) * (mas[0, 4] - (mas[0, 1] * x02 + mas[0, 2] * x03 + mas[0, 3] * x04));
      x01 = x1;
      x2 = (1 / mas[1, 1]) * (mas[1, 4] - (mas[1, 0] * x01 + mas[1, 2] * x03 + mas[1, 3] * x04));
      x02 = x2;
      x3 = (1 / mas[2, 2]) * (mas[2, 4] - (mas[2, 0] * x01 + mas[2, 1] * x02 + mas[2, 3] * x04));
      x03 = x3;
      x4 = (1 / mas[3, 3]) * (mas[3, 4] - (mas[3, 0] * x01 + mas[3, 1] * x02 + mas[3, 2] * x03));
      x04 = x4;
      // n++;
      sig = Max(Math.Abs(x1 - x01), Math.Abs(x2 - x02), Math.Abs(x3 - x03), Math.Abs(x4 - x04));
      while (sig > eps);
      Console.WriteLine("Корни СЛАУ \n x1= {0} \n x2= {1} \n x3={2} \n x4={3} ", x1, x2, x3, x4);
      }
      public static float Max( x, y, z, f)
      {
     x= Math.Max(x, y);
     x = Math.Max(x, z);
     x = Math.Max(x, f);
    }
  }
}