sample4_6.cs

来自「C#的数值计算的几乎所有算法」· CS 代码 · 共 51 行

/*
 * 示例程序Sample4_6: LEquations类的求解三对角线方程组的追赶法
 */

using System;
using CSharpAlgorithm.Algorithm;

namespace CSharpAlgorithm.Sample
{
	class Class1
	{
		[STAThread]
		static void Main(string[] args)
		{
            // 系数矩阵数据
            double[] mtxDataCoef5 = {
                                        13,12,0,0,0,
                                        11,10,9,0,0,
                                        0, 8,7,6,0,
                                        0, 0,5,4,3,
                                        0, 0,0,2,1};
            // 常数矩阵数据
            double[] mtxDataConst5 = {
                                         3.0,
                                         0.0,
                                         -2.0,
                                         6.0,
                                         8.0};
		
            // 构造系数矩阵
            Matrix mtxCoef5 = new Matrix(5, mtxDataCoef5);
            // 构造常数矩阵
            Matrix mtxConst5 = new Matrix(5, 1, mtxDataConst5);

            // 构造线性方程组
            LEquations leqs5 = new LEquations(mtxCoef5, mtxConst5);
		
            // 求解三对角线方程组的追赶法
            Matrix mtxResult5 = new Matrix();
            if (leqs5.GetRootsetTriDiagonal(mtxResult5))
            {
                Console.WriteLine(mtxResult5);
            }
            else
            {
                Console.WriteLine("失败");
            }
        }
	}
}