From d2d2fed4e484b5c2caea652838becd884a8225d5 Mon Sep 17 00:00:00 2001 From: shigek Date: Fri, 28 Mar 2003 05:00:21 +0000 Subject: Copied from rough/bigdecimal,documents & some sample programs added. git-svn-id: http://svn.ruby-lang.org/repos/ruby/trunk@3625 b2dd03c8-39d4-4d8f-98ff-823fe69b080e --- ext/bigdecimal/lib/linear.rb | 46 ++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 46 insertions(+) create mode 100644 ext/bigdecimal/lib/linear.rb (limited to 'ext/bigdecimal/lib/linear.rb') diff --git a/ext/bigdecimal/lib/linear.rb b/ext/bigdecimal/lib/linear.rb new file mode 100644 index 000000000..f93404fb6 --- /dev/null +++ b/ext/bigdecimal/lib/linear.rb @@ -0,0 +1,46 @@ +#!/usr/local/bin/ruby + +# +# linear.rb +# +# Solves linear equation system(A*x = b) by LU decomposition method. +# where A is a coefficient matrix,x is an answer vector,b is a constant vector. +# +require "bigdecimal" +require "ludcmp" + +include LUSolve + +def rd_order + printf("Number of equations ?") + n = gets().chomp.to_i +end + +zero = BigDecimal::new("0.0") +one = BigDecimal::new("1.0") + +while (n=rd_order())>0 + a = [] + as= [] + b = [] + printf("\nEnter coefficient matrix element A[i,j]\n"); + for i in 0...n do + for j in 0...n do + printf("A[%d,%d]? ",i,j); s = gets + a <<=BigDecimal::new(s); + as<<=BigDecimal::new(s); + end + printf("Contatant vector element b[%d] ? ",i);b<<=BigDecimal::new(gets); + end + printf "ANS=" + x = lusolve(a,b,ludecomp(a,n,zero,one),zero) + p x + printf "A*x-b\n" + for i in 0...n do + s = zero + for j in 0...n do + s = s + as[i*n+j]*x[j] + end + p s-b[i] + end +end -- cgit