Octaveはsymbolicで対応しているらしい?が、Ubuntu 14.04のでは未対応みたい。Sympyは対応している。 これで楽になるかな…と思ったけど、結果を見るとうーむ。 補足すると、以下のが、 >>> t = Symbol('t') >>> f = 0.3635819 + (0.4891775 + 0.0106411) * t…

数学は強くないので、変なところがあるかも。hypot(x, y)とは、sqrt(x*x+y*y)のことだけども、xが正なら、以下のように変形できる。 sqrt(x*x+y*y) = sqrt((y/x)^2+1)*xq=sqrt(k^2+1)は、k (つまりy/x) の値が大きくなるほど+1の要素が小さくなるため、kの値…