概説微分積分の問題3.19の解答。exp(0.5)、cos(1)、log(1.1)を小数第二位まで求める問題。ラグランジュの剰余項Rn(x)の絶対値を評価する必要があるのだけど、評価方法を経験則的にやってるので一般化できてない。 import math def trunc(x): if x > 0: retur…

概説微分積分の問題3.21の解答。x^3 + x^2 - 3 = 0 と x^4 - x - 1 = 0 を解く問題。 def f1(x): return float(x) ** 3 + float(x) ** 2 - 3 def df1(x): return 3. * float(x) ** 2 + 2. * x def f2(x): return float(x) ** 4 - float(x) - 1 def df2(x): r…

概説微分積分の問題3.16から3.18までの解答。ソフト使って解く問題だったのでSymPyを使ってみた。機能を探すのにかなり手間取ったけど、数式を代数的に扱えて面白い。∞はooだったりする。SymPyとnumpyとscipyとの関連性はまだよくわからず。ただしPlotはpygl…