36件ヒット
[1-36件を表示]
(0.009秒)
ライブラリ
-
bigdecimal
/ ludcmp (24) -
bigdecimal
/ newton (12)
検索結果
先頭3件
-
LUSolve
. # lusolve(a , b , ps , zero = 0 . 0) -> [BigDecimal] (9101.0) -
LU 分解を用いて、連立1次方程式 Ax = b の解 x を求めて返します。
...要素を Row-major order で並べて 1 次元の配列にし、
LUSolve.#ludecomp で変換したものを指定します。
@param b ベクトルを BigDecimal の配列で指定します。
@param ps LUSolve.#ludecomp の返り値を指定します。
@param zero 0.0 を表す値を......al/ludcmp'
include LUSolve
a = [['1.0', '2.0'], ['3.0', '1.0']].flatten.map(&:to_d)
# x = ['1.0', -1.0']
b = ['-1.0', '2.0'].map(&:to_d)
zero = '0.0'.to_d
one = '1.0'.to_d
# 以下の 2 行は
ps = ludecomp(a, b.size, zero, one) # a が破壊的に変更される
x = lusolve(a, b, ps, zero)
#......こう書いてもよい
# x = lusolve(a, b, ludecomp(a, b.size, zero, one), zero)
p x.map(&:to_f) #=> [1.0, -1.0]
//}... -
Newton
. # nlsolve(f , x) -> Integer (6101.0) -
ニュートン法を用いて非線形方程式 f(x) = 0 の解 x を求めます。
ニュートン法を用いて非線形方程式 f(x) = 0 の解 x を求めます。
@param f 関数を表すオブジェクトを指定します。詳細は
bigdecimal/newton をご覧ください。
@param x 探索を開始する点を数値の配列で指定します。解が複数ある場合、初
期値によって得られる解が異なります。また,初期値によっては収束
せずに RuntimeError が発生する場合があります。実行後は引
数 x に指定したオブジェクトに解が代入されます。
実行後は解を表す値が代入されています。
... -
LUSolve
. # ludecomp(a , n , zero = 0 , one = 1) -> [Integer] (3001.0) -
n 次正方行列を表す配列 a を破壊的に変更し、返り値と併せて元の行列の LU 分解を提供します。
...n 次正方行列を表す配列 a を破壊的に変更し、返り値と併せて元の行列の LU 分解を提供します。
LUSolve.#lusolve の引数を構築するために使います。
@param a 行列を BigDecimal の配列で指定します。(各要素を
Row-major order...