szeja_kod

            N = 5;
            numericUpDown1.Value = N;
            UstawTablice();
            int blad;
            AZ = new Complex[N + 1, N + 1];
            BZ = new Complex[N + 1];
            XZ = new Complex[N + 1];

            Complex E = 10;
            double Rw = 10, R1 = 10, R3 = 10, R4 = 10;
            double C2 = 10, C3 = 10;
            double L1 = 10;


            AZ[1, 1] = 1; AZ[1, 2] = -1; AZ[1, 3] = -1; AZ[1, 4] = 0; AZ[1, 5] = 0;
            AZ[2, 1] = 0; AZ[2, 2] = 0; AZ[2, 3] = 1; AZ[2, 4] = -1; AZ[2, 5] = -1;
            AZ[3, 1] = Rw; AZ[3, 2] = new Complex(R1, L1); AZ[3, 3] = 0; AZ[3, 4] = 0; AZ[3, 5] = 0;
            AZ[4, 1] = 0; AZ[4, 2] = new Complex(R1, L1); AZ[4, 3] = new Complex(0, C2); AZ[4, 4] = new Complex(-R3, C3); AZ[4, 5] = 0;
            AZ[5, 1] = 0; AZ[5, 2] = 0; AZ[5, 3] = 0; AZ[5, 4] = new Complex(R3, -C3); AZ[5, 5] = -R4;

            BZ[1] = Complex.Zero;
            BZ[2] = Complex.Zero;
            BZ[3] = E;
            BZ[4] = Complex.Zero;
            BZ[5] = Complex.Zero;

            for (int i = 1; i <= N; i++)
            {
                for (int j = 1; j <= N; j++)
                {
                    MacierzA.Rows[i-1].Cells[j-1].Value = AZ[i,j].ToString();
                }
                WektorB.Rows[i - 1].Cells[0].Value = BZ[i].ToString();
            }
            if (radioButton1.Checked)
            {
                blad = MetodaGaussa.RozRowMacGaussa(AZ, BZ, XZ, 1e-30);
            }
            else
            {
                blad = MetodaRozkladuLU.RozRowMacCroutDoolitta(AZ, BZ, XZ, 1e-30);
            }
            if (blad == 0)
                for (int i = 1; i <= N; i++)
                    WektorX[0, i - 1].Value = XZ[i].ToString("0.000000");


            textBox1.Text = "IpK bilans: " + ((XZ[1] - XZ[2] - XZ[3]) + (XZ[3] - XZ[4] - XZ[5])).ToString();
            textBox2.Text = "IIpK bilans: " + (E - Rw*XZ[1] - (R1 + Complex.j * L1)*XZ[2]).ToString();

            double Pz = (E*XZ[1].Conjugate).Re;
            textBox3.Text = "  Pz = " + Pz.ToString();


            double Podb = (XZ[1] * XZ[1].Conjugate * Rw).Re + (XZ[2] * XZ[2].Conjugate * (R1 + Complex.j * L1)).Re + (XZ[3] * XZ[3].Conjugate * -Complex.j * C2).Re + (XZ[4] * XZ[4].Conjugate * (R3 - Complex.j * C3)).Re + (XZ[5] * XZ[5].Conjugate * R4).Re;
            textBox4.Text = "Podb = " + Podb.ToString();

            Complex Sz = new Complex(0, 0);
            Sz = E * XZ[1].Conjugate;
            textBox5.Text = "  Sz= " + Sz.ToString();

            Complex Sodb = new Complex(0, 0);
            Sodb = (XZ[1].Abs * XZ[1].Abs * Rw) + (XZ[2].Abs * XZ[2].Abs * (R1 + Complex.j * L1)) + (XZ[3].Abs * XZ[3].Abs * (-Complex.j * C2)) + (XZ[4].Abs * XZ[4].Abs * (R3 - Complex.j * C3)) + (XZ[5].Abs * XZ[5].Abs * R4);
            textBox6.Text = "Sodb = " + Sodb.ToString();