ＶＢ２０１０での連立方程式の計算

質問

ＶＢ２０１０でガウスの消去法のプログラムを作りたいんですが、うまく解が求まらないのでメールしました。今回の例は未知数８です。ちなみにＶＢは初心者です。
Myarrayは、列の１－８個目までが、各未知数の係数、９個目が整数を表しています。
行は式を表していて、未知数が８個あるので８行あります。

Dim MyArray(,) As Double = {{-1, 1, -1, 1, 0, 0, 0, 0, 1},
{0, 0, 0, 1, 0, 0, 0, 0, -1},
{0, 0, 0, 0, 0, 0, 0, 1, -1},
{0, 0, 0, 0, 1, 1, 1, 1, 0},
{0, 0, 1, 0, 0, 0, -1, 0, 0},
{0, 2, 0, 0, 0, -2, 0, 0, 0},
{-6, 2, 0, 0, 0, 0, 0, 0, 0},
{6, 2, 0, 0, 0, 0, 0, 0, 0}
}

Const N = 8

Dim i, j, k, l, pivot As Integer
Dim x(N) As Double
Dim p, q, m As Double
Dim b(1, N + 1) As Double

For i = 0 To N - 1 Step 1
m = 0
pivot = i

For l = i To N - 1 Step 1 'i列の中で一番値が大きい列を選ぶ

If System.Math.Abs(MyArray(l, i)) > m Then
m = System.Math.Abs(MyArray(l, i))
pivot = l
End If

Next
'pivotがiと違えば、行の入れ替え
If pivot <> i Then
For j = 0 To 8 Step 1
b(0, j) = MyArray(i, j)
MyArray(i, j) = MyArray(pivot, j)
Myarray(pivot, j) = b(0, j)
Next
End If

Next

For k = 0 To N - 1 Step 1
p = Myarray(k, k) '対格要素を保存
MsgBox(p)
MyArray(k, k) = 1 '対格要素は１になることが分かっているので代入

For j = k + 1 To N Step 1
MyArray(k, j) = MyArray(k, j) / p
Next

For i = k + 1 To N - 1 Step 1
q = MyArray(i, k)

For j = k + 1 To N Step 1
MyArray(i, j) = MyArray(i, j) - q * MyArray(k, j)
Next

'0となることが分かってるので代入
MyArray(i, k) = 0
Next
Next

'解の計算
For i = N - 1 To 0 Step -1
x(i) = MyArray(i, N)

For j = N - 1 To i + 1 Step -1
x(i) = x(i) - MyArray(i, j) * x(j)
Next
Next

MsgBox(x(0))
MsgBox(x(1))
MsgBox(x(2))
MsgBox(x(3))
MsgBox(x(4))
MsgBox(x(5))
MsgBox(x(6))
MsgBox(x(7))

出力結果としてNAN(非数値）と出てきてしまいます。原因は、pに０が入ってしまっていることがあるからだと思いますが、これを回避する方法はないでしょうか？ご教授よろしくお願いします。

回答No.1

nag0720 · ベストアンサー

エラーの原因は行の入れ替えを最初の１回しかしていないからです。

行の入れ替えは、各行の計算ごとに行う必要があります。

For i = 0 To N - 1 Step 1
　　m = 0
　　pivot = i
　　・・・・・
　　・・・・・
Next

For k = 0 To N - 1 Step 1
　　p = Myarray(k, k) '対格要素を保存
　　・・・・・
　　・・・・・
Next

この部分を、

For k = 0 To N - 1 Step 1
　　m = 0
　　pivot = k
　　・・・・・
　　・・・・・

　　p = Myarray(k, k) '対格要素を保存
　　・・・・・
　　・・・・・
Next

とすれば大丈夫でしょう。
初めのFor文のインデックスをiからkに変えているので注意してください。

補足

解答ありがとうございます。修正プログラムfor分の終わりのNEXTが元のプログラムの２つ目のfor分の終わりのＮＥＸＴの場所でよいでしょうか？

色々やってみましたが if pivot <> i thenの中の
b(0,j) = MyArray(i,j)
の所でインデックスが配列の境界外です、ってエラーが出てしまいますね。解答の意味を勘違いしてたら申し訳ないです＞＜

回答No.2

nag0720

＞修正プログラムfor分の終わりのNEXTが元のプログラムの２つ目のfor分の終わりのＮＥＸＴの場所でよいでしょうか？

それでいいです。要するに２つのFor文を１つにまとめただけです。

＞色々やってみましたが if pivot <> i thenの中の
＞b(0,j) = MyArray(i,j)
＞の所でインデックスが配列の境界外です、ってエラーが出てしまいますね。

For文の変数iをkに変えたので、１つめのFor文の中のiをすべてkに変えてください。

お礼

よく考えたらそうですよね。不注意失礼しました＞＜

無事に動きました。本当にありがとうございました！