11 解多元线性方程组
require 'matrix'
m = Matrix[[2, 3], [1, 4]]
p m.lup_decomposition.solve([2, 3]).to_a
#=> [(-1/5), (4/5)]
12 分解质因数 (其实实现比较简单,大数还是 mathematica 内建的算法比较快)
require 'mathn'
p 123232913892382.prime_division
#=> [[2, 1], [7, 1], [17, 1], [4519, 1], [114579631, 1]]