Я хотел бы использовать модуль bigdecimal/newton
в Ruby для реализации XIRR .Я написал скрипт, чтобы попробовать его, следуя этому примеру .
Когда я запускаю код (Ruby 1.9.2 в Mac OS X 10.6), я получаю следующую ошибку:
/opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:84:in `div': wrong number of arguments(2 for 1) (ArgumentError)
from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:84:in `block in lusolve'
from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:78:in `downto'
from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/ludcmp.rb:78:in `lusolve'
from /opt/local/lib/ruby1.9/1.9.1/bigdecimal/newton.rb:60:in `nlsolve'
from gistfile1.rb:52:in `<main>'
Это из-за ошибки в моей программе или из-за чего-то в библиотеке BigDecimal
?
#!/opt/local/bin/ruby1.9
require 'bigdecimal'
require 'bigdecimal/newton'
require 'time'
include Newton
transactions = [
{amount:BigDecimal::new("-1000"), date:Time.parse("2011-01-01")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-02-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-03-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-04-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-05-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-06-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-07-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-08-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-09-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-10-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-11-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2011-12-11")},
{amount:BigDecimal::new("100"), date:Time.parse("2012-01-11")},
]
class Function
values ={ eps: "1.0e-16", one: "1.0", two: "2.0", ten: "10.0", zero: "0.0" }
values.each do |key, value|
define_method key do
BigDecimal::new value
end
end
def initialize(transactions)
@transactions = transactions
end
def values(x)
start = @transactions[0][:date]
value = []
value << @transactions.reduce(0) do |sum, t|
pwr = (t[:date] - start) / 365
r = t[:amount] / (1.0 + x[0]) ** pwr
sum + r
end
value
end
end
f = Function.new(transactions)
x = [f.zero]
n = nlsolve(f,x)
puts x
Любая помощь при отладке будет принята с благодарностью.
Спасибо!