[問題] 高斯消去法的問題
大家好!
最近在寫一個高斯消去法的C++ code,
但是邏輯上已經檢查過很多次了,
還是不知道發生了什麼問題,
答案跟用手算的答案不同。
我的原始array是
a[0][0] = 1;
a[0][1] = 1;
a[0][2] = 3;
a[1][0] = 1;
a[1][1] = -1;
a[1][2] = 1;
也就是 x1 + x2 = 3
x1 - x2 = 1
用以下的程式碼去跑:
t = 0;
double factor = 0;
do
{
for(j=0; j<=n; j++)
{
a[t][j]=a[t][j]/a[t][t];
}
for(i=0; i<= n-1; i++)//row
{
if (i != t)
{
for(j=0; j<=n; j++)//column
{
factor = (a[i][t] / a[t][t]);
a[i][j] = a[i][j]-( factor * a[t][j]);
}
}
}
t++;
}
while(t<n);
拜託大家了!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.107.241
→
06/29 22:51, , 1F
06/29 22:51, 1F
推
06/30 00:47, , 2F
06/30 00:47, 2F
→
06/30 00:47, , 3F
06/30 00:47, 3F
→
06/30 01:41, , 4F
06/30 01:41, 4F
→
06/30 01:41, , 5F
06/30 01:41, 5F
推
06/30 04:04, , 6F
06/30 04:04, 6F
→
06/30 11:48, , 7F
06/30 11:48, 7F
→
06/30 17:30, , 8F
06/30 17:30, 8F
→
07/10 18:23, , 9F
07/10 18:23, 9F
→
07/10 18:23, , 10F
07/10 18:23, 10F
C_and_CPP 近期熱門文章
PTT數位生活區 即時熱門文章