Scilab

連続系アルゴリズム演習「常微分方程式を、刻み幅制御を行う線形多段階法で解け」なる課題にScilabを使ってみた。数学ソフトに触るのは初めて。これ、面白いかも。

function qx = integral(px) // integral constant is 0
  coeffs = coeff(px)
  qx = poly([0, coeffs ./ (1:length(coeffs))], 'x', 'coeffs')
endfunction

4行で多項式の積分が書けた。Haskellの2行quicksort見た時と似た感覚。
変数宣言不要だからtypoが怖いとかあるけど、面白いから次の課題でも使おう。