1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
| #! /usr/local/bin/ruby
#*********************************************
# 連立方程式の解法 ( ガウス・ジョルダン法 )
#*********************************************
#
class GaussJorden
def initialize
# 係数
@a = [
#[ 2, -3, 1, 5],
#[ 1, 1, -1, 2],
#[ 3, 5, -7, 0]
[ 1, -2, 3, -4, 5],
[-2, -5, -8, -3, -9],
[ 5, 4, 7, 1, -1],
[ 9, 7, 3, 5, 4]
]
# 次元の数
@n = @a.size
end
# 計算・結果出力
def exec
# 元の連立方程式をコンソール出力
display_equations
@n.times do |k|
# ピボット係数
p = @a[k][k]
# ピボット行を p で除算
k.upto(@n) { |j| @a[k][j] /= p.to_f }
# ピボット列の掃き出し
@n.times do |i|
next if i == k
d = @a[i][k]
k.upto(@n) { |j| @a[i][j] -= d * @a[k][j] }
end
end
# 結果出力
display_answers
rescue => e
raise
end
private
# 元の連立方程式をコンソール出力
def display_equations
0.upto(@n - 1) do |i|
0.upto(@n - 1) {|j| printf("%+dx%d ", @a[i][j], j + 1)}
printf("= %+d\n", @a[i][@n])
end
rescue => e
raise
end
# 結果出力
def display_answers
0.upto(@n - 1) {|k| printf("x%d = %f\n", k + 1, @a[k][@n])}
rescue => e
raise
end
end
if __FILE__ == $0
begin
# 計算クラスインスタンス化
obj = GaussJorden.new
# 連立方程式を解く(ガウス・ジョルダン法)
obj.exec
rescue => e
$stderr.puts "[#{e.class}] #{e.message}"
e.backtrace.each{ |tr| $stderr.puts "\t#{tr}" }
end
end
|