Computes the Jacobian of f at x. fx is the
value of f at x.
# File ext/bigdecimal/lib/bigdecimal/jacobian.rb, line 76
def jacobian(f,fx,x)
n = x.size
dfdx = Array::new(n*n)
for i in 0...n do
df = dfdxi(f,fx,x,i)
for j in 0...n do
dfdx[j*n+i] = df[j]
end
end
dfdx
end