[완료] 간단한 루비 질문입니다.

require 'matrix'

print "Enter the dimension of the matrix:"
begin
dim=Integer(gets)
rescue ArgumentError
puts "This is not an integer!"
end

class Matrix
#this method is required due to poor definition of matrices in ruby
def Matrix.diagonal(values)
size=values.size
rows=(0..size-1).collect {
|j|
row=Array.new(size).fill(0, 0, size)
row[j]=values[j]
row
}
rows(rows, false)
end
end

def dot_product(vec1, vec2)
sum=0
0.upto(vec1.size-1) {|i| sum+=vec1[i]*vec2[i]}
return sum
end

def norm(vec) Math.sqrt(dot_product(vec, vec)) end

def normalize(vec) vec.map{|i| i/norm(vec)} end

def gen_rand_vecs(dim) Array.new(dim) {Array.new(dim) {rand}} end

def GramSchmidt(vecs)
dim=vecs.size
new_vecs=[]
new_vecs << normalize(vecs[0])
(1...dim).each do |i|
new_vec=vecs[i].clone
(0...i).each do |j|
dot=dot_product(new_vec, new_vecs[j])
(0...dim).each {|k| new_vec[k]-=new_vecs[j][k]*dot}
end
new_vecs << normalize(new_vec)
end
return new_vecs
end

def print_mat(mat)
dim=mat.rank
(0...dim).each do |i|
puts Array.new(dim) {|j| sprintf("%4f", mat[i, j])}.join("\t")
end
end

diag_mat=Matrix.diagonal(Array.new(dim) {rand})
rot_mat=Matrix.rows(GramSchmidt(gen_rand_vecs(dim)))
#puts rot_mat.det
mat=rot_mat.t*diag_mat*rot_mat

puts "Diagonal matrix:"
print_mat(diag_mat)
puts "\nRotation matrix:"
print_mat(rot_mat)
puts "\nRotated matrix:"
print_mat(mat)

혹시 나중에 참고하실 분이 계신가 하여...

루비에서 Vector와 Matrix 클래스는 부실하군요.