24小时热门版块排行榜    

查看: 835  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +8 yangfz 2026-03-17 8/400 2026-03-18 10:43 by yjolah
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +3 1孙悟空 2026-03-17 3/150 2026-03-18 10:28 by macy2011
[考研] 288求调剂,一志愿华南理工大学071005 +3 ioodiiij 2026-03-17 3/150 2026-03-18 09:58 by 求调剂zz
[考研] 312求调剂 +7 陌宸希 2026-03-16 8/400 2026-03-18 09:54 by macy2011
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
信息提示
请填处理意见