text_070_diff_integral

微分

高校で習う微分の公式は、以下の２つです。

• $\frac{d}{dx} c = 0$
• $\frac{d}{dx} x^n = n x^{n-1}$

sageで微分は、diff関数を使います。diff関数の使い方は以下の通りです。次数を省略すると1階微分となります。

diff(関数, 微分する変数, 次数)

Sageで微分公式を確認

sageで変数xの関数fを以下のように定義します。

x = var('x')
f = function('f', x)
	

以下に示す微分の公式をSageで実行し、f'=D[0](f)(x)の置換を実行した結果を示します。

• $(cf(x))' = cf'(x)$
• $(f(x) + g(x))' = f'(x) + g'(x)$
• $(f(x)g(x))' = f'(x)g(x) + f(x)g'(x)$
• $(f(x)/g(x))' = (f(x)'g(x) - g'(x)f(x))/g(x)^2$
• $(1/g(x))' = -g'(x)/g(x)^2$
• $(f(g(x)))'= f'(g(x))g'(x)

いろんな関数の微分

以下に主な関数を微分した結果を示します。

微分の応用

微分の応用として、接線の方程式とその法線を計算してみましょう。 関数f(x)を以下のように定義します。 $$ f(x) = x^3 - x^2 - 2*x $$ このとき、$x = x_1$での接線の式は、以下のようになります。 $$ y - y_{1} = f'(x_{1})(x - x_{1}) $$ また、法線は接線と直交することから、以下のようになります。 $$ y - y_{1} = -\frac{1}{f'(x_{1})}(x - x_{1}) $$ f(x)の微分した式をg(x)とし、x=0での接線と法線を計算し、プロットしてみます。

積分

sageでは積分は、integrate関数で計算します。

		integrate(被積分関数, 積分変数)
		また
		integrate(被積分関数, 積分変数, 積分範囲)
		とすると定積分を計算します。
		

例として、$\int x dx$を計算した後、その微分を求めています。 元のxが返されています。

いろんな関数の積分

以下に主な関数を積分した結果を示します。

定積分

次に定積分$\int_{0}^{\pi/2}sin(x) dx$の例を示します。

$sin(x)$を0から$\pi/2$の範囲でプロットすると以下のようになります。

sinの積分は、-cosですから、定積分は以下のようになります。 $$ [ -cos(x) ]_0^\pi = \{ 0 - (-1) \} = 1 $$

それでは、Sageを使って、$sin(x)$を0から$\pi/2$の範囲で積分してみましょう。

もう一つ、以下の関数を積分した例をみてみましょう。 $$ \int_{0}^{3}x^2 + 2sin(x) dx $$

先ほどと同様に、積分範囲をプロットしてみます。

積分の結果は、関数で返されるのでN(_)を使って数値に変換します。

このように関数のプロット結果から推定される概算と積分の結果をみながら、 計算結果を確認しながら、進められるのもSageの特徴です。

数値積分

すべての関数が定積分できるとは限りません、そこでSageには数値積分を行う numerical_integral関数が用意されています。

numerical_integral関数の使い方は、以下の通りです。 関数の戻り値は、積分結果と誤差のタプルが返されます。

	numerical_integral(被関分館数, 積分範囲)
	

例として、$sin(x)$を0から$\pi/2$で数値積分した結果を以下に示します。

残念ながら、integrate関数を使うと符号がマイナスの原因が不明ですが、-16が返ってきます。

次に数値積分で計算してみます。16に近い値が求まります。

Sageのような数式処理で計算する場合には、途中計算のチェックや他の手法との相互 チェックをするように心がけるとよいでしょう。