N88-BASICで連立方程式

2022/11/3(土)
N88-BASICで連立方程式
 
例えば
次の連立3元1次方程式
 
a₁x + b₁y + c₁z = d₁ 
a₂x + b₂y + c₂z = d₂ 
a₃x + b₃y + c₃z = d₃ 
 
は、行列で
|a₁ b₁ c₁||x| = |d₁|
|a₂ b₂ c₂||y| = |d₂|
|a₃ b₃ c₃||z| = |d₃|
と書け
対角成分を1、それ以外を0にすると
|1  0  0||x| = |d₁'|
|0  1  0||y| = |d₂'|
|0  0  1||z| = |d₃'|
つまり
|x| = |d₁'|
|y| = |d₂'|
|z| = |d₃'|
となり、解くことが出来る
 
プログラムでは
|a₁ b₁ c₁ d₁|
|a₂ b₂ c₂ d₂|
|a₃ b₃ c₃ d₃|
の
1行目をa₁で割りa₁=1とし
1行目をa₂倍し2行目から引きa₂=0とし
1行目をa₃倍し3行目から引きa₃=0とする
これを、b₂,c₃についても行うことで
対角成分を1、それ以外を0にしています
 
この方法には
対角成分に0がある式は解けないという
欠点がありますので
対角成分に0が来ない様に式の順番を
手動で変えてから実行して下さい
 
 
N88-BASIC互換?VL,NL,XL-BASICと
blg～.zip(simu001.bas)は
以下のリンクからダウンロードできます

https://ulprojectmail.blogspot.com
Readme.txtを読んで遊んで下さい