2014-11-02から1日間の記事一覧

山登り法

山登り法は解を近似するために使用します. 局所解に陥ってしまうので, 注意が必要です. このように, x軸を分割し, 微分値が0の時が最適解となります. 昨日の式と同じ, -x^2+4*x を用います. from sympy import * def euler(x0, y0): e =0.01 x = x0 y = y0 h…

簡単なプログラム5

研究でも使っている微分と積分についてです. 非常に難しい計算もすぐに値を返してくれます. まず微分です. 微分はdiff()関数を用います. ここでsin(x)をxで微分してみましょう. from sympy import * #sympyの呼び出し x = Symbol('x') print diff(sin(x),x) …