24小时热门版块排行榜    

查看: 902  |  回复: 0

springwu

铁虫 (初入文坛)

[求助] 大家帮忙看一下这个C#求解线性方程组的程序那里有问题?

程序编译没错误,可是算出来的结果不对。小弟才学习C#还请多多指教!
namespace 解线性方程组
{
    class Program
    {
        static void Main(string[] args)
        {
            int amount;//存储方程组的方程个数
            double[,] m;//二维数组用于表示存放方程系数的系数矩阵
            double[] solution;//一维数组用于存放方程组的解
            double t1, t2;
            Console.WriteLine("Input number of equation";//接收用户输入的方程个数
            amount = Convert.ToInt32(Console.ReadLine());
            m = new double[amount, amount + 1];
            solution = new double[amount];

            Console.WriteLine("Input coefficient of equation";
            for (int i = 0; i < amount; i++)
            {
                for (int j = 0; j < amount + 1; j++)
                    m[i, j] = Convert.ToDouble(Console.ReadLine());
            }

            Console.WriteLine("The coefficient of equation is:";
            Console.WriteLine("-------------------------------";
            for (int i = 0; i < amount; i++)
            {
                for (int j = 0; j < amount + 1; j++)
                    Console.Write("{0,8}", m[i, j]);
                Console.WriteLine();
            }
            for (int i = 1; i < amount; i++)
            {
                for (int j = i; j < amount; j++)
                {
                    double coef = -m[j, i - 1] / m[i - 1, i - 1];
                    for (int k = i - 1; k < amount + 1; k++)
                    {
                        m[j, k] += m[i - 1, k] * coef;
                    }
                }
            }
            t1 = m[amount - 1, amount - 1];
            solution[amount - 1] = m[amount - 1, amount] / t1;
            for (int i = amount - 1; i >= 0; i--)
            {
                t2 = 0d;
                for (int j = i + 1; j <= amount - 1; j++)
                {
                    t2 += m[i, j] * solution[j];
                }
                solution = (m[i, amount] - t2) / m[i, i];
            }
            Console.WriteLine();
            Console.WriteLine("\n The result is:";
            Console.WriteLine("-----------------------";
            for (int i = 0; i < amount; i++)
            {
                Console.WriteLine("x[{0}] = 1", i, solution);
            }
            Console.Read();
        }
        
    }

}
回复此楼
让人迷茫的原因只有一个,那就是本该拼搏的年纪,却想得太多,做得太少
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 springwu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 7/350 2026-11-02 12:46 by n0P4fnsXKbBK
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 10:02 by dcc2yOIRGGTT
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:13 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 xvJS1BYDxvBS 2026-10-31 13/650 2026-11-02 06:03 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-02 06:00 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 05:54 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 05:47 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
信息提示
请填处理意见