24小时热门版块排行榜    

Znn3bq.jpeg
查看: 850  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +11 1234567wang 2026-05-17 13/650 2026-05-20 18:55 by tanfancy
[基金申请] 评审有感 +14 popular289 2026-05-18 23/1150 2026-05-20 17:39 by lincyb
[论文投稿] Sci. Bull. 悲剧经验 +6 jyang1999 2026-05-16 7/350 2026-05-20 16:39 by feng6531
[基金申请] 今年审到国自然15份,谈谈感受 +19 国自然国社科中 2026-05-17 20/1000 2026-05-20 14:14 by 仲夏夜的星星
[基金申请] 国自然评分 +3 无名者登山 2026-05-20 4/200 2026-05-20 13:29 by vito刘
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 8/400 2026-05-20 09:32 by xtlilibin
[基金申请] 河北省自然科学基金 +3 Peterchao 2026-05-18 3/150 2026-05-20 08:57 by 霸_霸
[教师之家] 上海大学实验技术岗位非升即走 +12 嘻嘻哈哈乐呵呵 2026-05-15 13/650 2026-05-20 08:34 by xli1984
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +3 光伏爱好者 2026-05-17 5/250 2026-05-19 19:13 by 光伏爱好者
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[考博] 找博士生导师 +6 小代想上岸 2026-05-15 7/350 2026-05-19 10:22 by free_fisher
[基金申请] 同样的基金本子,换个专家直接从C变A! (金币-10) +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[基金申请] 别被青基扩招骗了!26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
信息提示
请填处理意见