[問題] 計算線性方程錯誤

[loser113 (洨大魯蛇ㄍ)]

我用 slove(A,b) 計算 A=nxn矩陣 =b 時 我的解不該出現負值

x1~xn >0

n最多到7千多

code檢查多遍無誤

發現在結尾Xn和前幾個Xk 解出數值小於10^-17  (|Xk|<10^-17)

就會有可能是負值 (像Xk=-2.357...*10^-25)

R在做計算的時候 有可能有這種誤差嗎謝謝

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 114.42.139.52
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1447599481.A.855.html
→ celestialgod: 這種誤差是有可能的
推 Wush978: 10^-17... 可能是浮點數運算的誤差
→ loser113: 謝謝
→ loser113: 有方法解決嗎
→ celestialgod: 沒有...