高斯消元法（02）

using System;

namespace 高斯消元法2
{
        class Program
        {
            static int n;
            static double[,] a;
            static double[] x;
            static void xiaoyuan()
            {
                for (int k = 0; k < n - 1; k++)
                {
                    for (int i = k + 1; i < n; i++)
                    {
                        double ak = a[i, k] / a[k, k];
                        for (int j = k; j < n + 1; j++)
                        {
                            a[i, j] = a[i, j] - a[k, j] * (ak);
                        }
                    }
                    Console.WriteLine("\n第{0}次消元：", k + 1);
                    for (int i = 0; i < n; i++)
                    {
                        string s = "";
                        for (int j = 0; j < n + 1; j++)
                        {
                            s += string.Format("{0,-8:f2}", a[i, j]);
                        }
                        Console.WriteLine(s);
                    }
                }
                Console.WriteLine("***********消元完成*********");
            }
            static void xiaoyuan1()
            {
                for (int k = 0; k < n - 1; k++)
                {
                    for (int i = k + 1; i < n; i++)
                    {
                        for (int j = k; j < n + 1; j++)
                        {
                            a[i, j] = a[i, j] - a[i, k] / a[k, k] * a[k, j];
                        }

                    }
                    Console.WriteLine("第{0}次消元！！！！！！！！！", k + 1);
                    for (int i = 0; i < n; i++)
                    {
                        string s = null;
                        for (int j = 0; j < n + 1; j++)
                        {
                            s += string.Format("  {0,-4}  ", a[i, j]);
                        }
                        Console.WriteLine(s);
                    }
                }
            }
            static void huidia()
            {
                for (int i = n - 1; i >= 0; i--)
                {
                    double sum = 0;
                    for (int j = i + 1; j < n; j++)
                    {
                        sum = a[i, j] * x[j];
                    }
                    x[i] = (a[i, n] - sum) / a[i, i];
                }
            }
            static void huidai()
            {
                for (int i = n - 1; i >= 0; i--)
                {
                    double sum = 0;
                    for (int j = i + 1; j < n; j++)
                    {
                        sum = a[i, j] * x[j];
                    }
                    x[i] = (a[i, n] - sum) / a[i, i];
                }
                Console.WriteLine("解为：");
                for (int i = 0; i < n; i++)
                {
                    Console.WriteLine("x{0}={1,6:f2}", i + 1, x[i]);
                }
            }
            static void input()
            {
                for (int i = 0; i < n; i++)
                {
                    string s = "";
                    for (int j = 0; j < n + 1; j++)
                    {
                        s += string.Format("{0,-6:f2}", a[i, j]);
                    }
                    Console.WriteLine(s);
                }
            }
            static void Main(string[] args)
            {
                Console.WriteLine("请输入阶数：");
                n = Convert.ToInt32(Console.ReadLine());
                a = new double[n, n + 1];
                x = new double[n];
                Console.WriteLine("请输入各行系数：");
                for (int i = 0; i < n; i++)
                {
                    string sa = Console.ReadLine();
                    string[] ss = sa.Split(' ');
                    for (int j = 0; j < n + 1; j++)
                    {
                        a[i, j] = Convert.ToDouble(ss[j]);
                    }
                }
                Console.WriteLine("矩阵系数为：");
                input();
                xiaoyuan();
                huidia();
                Console.WriteLine("方程解为：");
                for (int i = 0; i < n; i++)
                {
                    Console.WriteLine("x{0}={1,6:f2}", i + 1, x[i]);
                }
            }
        }
    
}