物理数学I 解析学のソースを表示

<small> [[物理数学I]] > 解析学</small>	
----
==解析学==

解析学は高校までの数学の延長としてとらえることも出来るが、高校までの数学を厳密に基礎づける科目ととらえることも出来る。例えば、高校までの範囲では数列の極限や関数の連続は厳密には定義されていなかった。解析学ではこのような極限を取る手法を扱う。また、微分や積分に関するより進んだ計算も扱う。ここで学んだ手法は線形代数と並んで、より進んだ計算を行なうための基礎となるので、ここで学ぶ手法には十分習熟する必要がある。

===1変数の計算===

ここでは、1つの変数を扱う関数を用いて収束や連続性の定義を扱う。また、それらを用いて厳密に定義された手法を用いてテイラー展開やより複雑な積分を導入する。

====実数の連続性====

最初に、無理数を定義する手法を考える。高校までの範囲では、実数のうちで有理数でないものを無理数と定義した。ここで有理数とは、2つの互いに素の整数n,mを用いて、
:<math>
\frac n m
</math>
とかかれるもの全体を指す。しかし、この構成ではそもそも実数が何なのかが示されていないため、無理数というものがとらえにくいという難点がある。
ここで、実数の性質について1つの仮定をおく。
:実数が全て書かれた直線を数直線とする。この数直線上でただ1点に対する切断を考えるとき、その点はその点より小さい数の集合と大きい数の集合を作りだす。このとき、この数自身は小さい数の集合に含まれて、大きい数の集合には含まれないものとする。
この定義はデデキントの切断([[w:en:Dedekind cut]])と呼ばれる。このとき、ある実数をその数より小さい有理数の集合によって定義する。この定義は有理数と無理数の両方に対して適用できる。なぜなら、切断で選ばれた点が有理数だったときには、その点自身までの有理数の集合を選んだ有理数を表わす有理数の集合として扱えばよい。一方、切断によって選ばれた点が無理数だったときには、その切断は必ずその近くにある別の数を表わす切断とは区別される。なぜなら、ある数を選んだときその数と別の数の間には必ずある有理数が存在するからである。有理数のこの性質は有理数の[[w:稠密]]性と呼ばれ、有理数の重要な性質である。これは、どんな数でも数値として書くならその値はどんな場合でも無限小数で書くことが出来、無限小数はどれほど小さい数でも有理数で書かれ循環小数を含んでいることから確かに成立するのである。このようにして、無理数はその数より小さい有理数全体の集合によってとらえられた。

====数列の収束の定義====

ここからは、上で述べた実数の連続性を用いて、数列の収束を定義する。まずは、収束の定義を述べる。任意の(小さい)ある数<math>\epsilon</math>をとったとき、あるNが存在してn <math>>=</math> N を満たす全てのnについて
:<math>
|a _n - a| < \epsilon
</math>
が成り立つとき数列<math>a _n</math>は、定数aに収束するという。

ここで、実数の連続性は無限にある定数aに近い数がただ1つしかないということを見るために用いられている。これは、ある定数aと異なった点bは、定数aとの間に何らかの有理数を持つため、定数aと無限に近くにあることは出来ない。そのため、数列<math>|a _n-a|</math>が、定数aと選んだ点bの距離よりも小さい<math>\epsilon</math>よりも小さいという条件を満たすとき、<math>a _n</math>が収束する点は確かに点bではなく、点aであることが
保証されるのである。上の定義は高校までに行なった極限の定義に適合しているはずなので、実際に極限の計算を行なうときには、これまでに用いた結果をそのまま用いてもよい。<!-- ただし、この定義によると収束する数列の和や、積に関する結果は、この定義から直接導出することができるため、以前よりも少ない仮定で計算が進められるといえる。-->この定義を用いたとき、以下が成り立つ。

定数a,bに収束する数列<math>a _n</math>,<math>b _n</math>に対して、

(I)
:<math>
\lim (a _n + b _n) = a + b
</math>
(II)
:<math>
\lim (a _n \times  b _n) = a b
</math>
が成り立つ。

*導出

(I)について、数列<math>a _n</math>がaに収束することから、ある定数<math>\epsilon _1</math>を取ったとき、ある定数<math>N _1</math>が存在し、<math>N _1 < n</math>を満たす全てのnについて、
:<math>
|a _n - a | < \epsilon  _1
</math>
が成立する。同様に数列<math>b _n</math>がbに収束することから、ある定数<math>\epsilon _2</math>を取ったとき、ある定数<math>N _2</math>が存在し、<math>N _2 < n</math>を満たす全てのnについて、
:<math>
|b _n -b | < \epsilon  _2
</math>
が存在する。

ここで、
:<math>
a _n + b _n
</math>
について、
:<math>
N =  \textrm{max} (N _1,N _2)
</math>
としたとき、全ての<math>n>N</math>を満たす整数nに対して
:<math>
| a _n + b _n - ( a+b)|
</math>
を計算すると、この量は三角不等式を用いることで、
:<math>
< | a _n - a | + |b _n - b| 
</math>
:<math>
< \epsilon  _1 + \epsilon  _ 2
</math>
が成り立つ。しかし、<math>\epsilon _1</math>,<math>\epsilon _2</math>はNを大きく取ることでいくらでも小さくできるため、全ての<math>\epsilon</math>に対して
:<math>
 \epsilon  _1 + \epsilon  _ 2 < \epsilon
</math>
となるような整数Nが存在する。よって、
:<math>
\lim ( a _n + b _n)  =  a+b
</math>
が示された。

(II)

同様に
:<math>
a _n b _n
</math>
について、
:<math>
| a _n b _n - ab|
</math>
は、
:<math>
= | a _n ( b _n - b) + b (a _n - a) |
</math>
:<math>
\le  |a _n ( b _n - b) | +  |b| |a _n - a | 
</math>
となる。ここで、<math>n>N</math>に対しては
:<math>
 a - \epsilon _1< a _n < a + \epsilon _1
</math>
が成り立つことに注目すると、
:<math>
< ( a + \epsilon _ 1 ) \epsilon _2 + |b| \epsilon _1
</math>
が得られる。ここで、<math>\epsilon _1</math>,<math>\epsilon _2</math>はNを大きく取ることでいくらでも小さくできるため、a,bが有限のときa,bの値に関わらず上の値は限りなく小さくなる。よって、
:<math>
\lim (a _n \times  b _n) = a b
</math>
が示された。

**問題

次の数列
:<math>
n  \rightarrow \infty 
</math>
の極限値を求めよ。

(I)
:<math>
\lim ( 1 + \frac 1 n) 
</math>
(II)
:<math>
\lim ( 2 \times \frac 1 n)
</math>

**解答
上の結果である

(I)
:<math>
\lim (a _n + b _n) = a + b
</math>
(II)
:<math>
\lim (a _n  b _n) = a b
</math>
を用いればよい。ただし、定数は全てのnに対して同じ数を取る数列として扱う。

(I)
:<math>
\lim ( 1 + \frac 1 n) 
</math>
は、1は極限値1をとり
:<math>
\frac 1 n
</math>
は、極限値0を取ることから、
:<math>
\lim ( 1 + \frac 1 n)  = 1 +0 =1
</math>
となる。

(II)
:<math>
\lim ( 2 \times \frac 1 n)
</math>
について、2は、極限値2を取り、
:<math>
\frac 1 n
</math>
は極限値0を取ることから、
:<math>
\lim ( 2 \times \frac 1 n) = 2 \times 0 = 0
</math>
が成り立つ。一般に定数倍や定数の足し算は、極限値に定数倍や定数の足し算をすればよい。


次に数列の発散の定義をする。ここでも上の場合と同様無限個の数列の値がある値より大きくなることが重要である。あるNが存在してn <math>\ge</math> N を満たすすべてのnについて任意に取った(大きい)Rに対して、
:<math>
a _n > R
</math>
が成り立つとき、<math>a _n</math>はn無限大で正の無限大に発散するという。このことを
:<math>
\lim a _n = \infty
</math>
と書かれる。

*問題例
**問題
:正の無限大に発散する例
数列
:<math>
a _n = n
</math>
の場合についてこの数列が上の定義を用いたときに正の無限大に発散することを示せ。

**解答
ここでも、Nの選び方が重要である。ここでは、あるRに対して
:<math>
N \ge R
</math>
と選べばよい。この場合、どのような(大きい)Rを取ったとしても
:<math>
N \ge R
</math>
を満たすような整数Nを選ぶと、それ以降の全てのnについて
:<math>
a _n = n  \ge R = N 
</math>
が成り立つ。値Rはいくらでも大きくできるので、このことは数列の発散の条件を満たしている。よって、数列
:<math>
a _n = n
</math>
はn無限大で正の無限大へと発散する。


同じ様にして、 あるNが存在してn <math>\ge</math> N を満たすすべてのnについて任意に取った(小さい)Rに対して、
:<math>
a _n < R
</math>
が成り立つとき、<math>a _n</math>はn無限大で負の無限大に発散するという。このことは
:<math>
\lim a _n = - \infty
</math>
と書かれる。

このうちのいずれにも当てはまらない場合もある。例えば、次の場合は数列はどの値に収束することもないため、数列は極限値を持たない。


*問題例
**問題
:ある有限値に収束しない場合
:<math>
a _n = (-1)^n
</math>
が上の定義のいずれも満たさないことを示し、この数列が収束も発散もしないことを導出せよ。


**解答
このとき、非常に大きなNを取ったとしても、そのNから先の全てのnについて<math>a _n</math>がきわめてaに近い値に留まるようなaは存在しない。例えば、a = 1と取ったとすると、ある値kにおいて
:<math>
a _n - a = 0
</math>
となり、両者は非常に近くなる。しかし、n=k+1においては既に、その値は-1となり、
:<math>
|a _n - a| = 2
</math>
となり、任意に小さい数<math>\epsilon</math>に対してより小さい数であり続けることはできない。これはどれほど大きなkをとっても、もしくはa = -1 もしくはそれ以外の量を選んでも同じである。よって、この数列はn無限大である値に収束することは無い。一方、この数列は1と-1しか値を取らないため、どのような数よりも大きくなるような数列ではない。よって、この数列は正負の無限大に発散することもない。よって、この数列は収束も発散もしないことが示された。

====連続の定義====

ある区間<math>I</math>において定義された関数<math>f</math>が<math>a \in I</math>で連続とは、<br/>
どんな<math>\epsilon >0</math>についても,ある<math>\delta>0</math>が存在して<br/>
<math>|x - a| <\delta</math>を満たす全ての<math>x(\in I)</math> について
:<math>
| f (x) - f(a)| < \epsilon
</math>
が成り立つことである。<br/>
区間Ｉの全ての点で連続のとき、関数fはＩ上で連続であると呼ぶ。

====複数回微分の定義====

n回微分を
<math>
f^{(n)} = (f^{(n-1)})'
</math>
で定義する。


====テイラー展開====

=====テイラー級数の定義=====

ある関数
f(x)について、fが定義された全ての実数について
:<math>
f(x) = f(a) + f'(a)(x -a ) + \frac 1 2 f''(a) (x - a)^2 + ... + 
\frac 1 {n!} f ^{(n)} (\xi ) (x -a)^n
</math>
が成り立つ。(<math>\xi</math>はaとxの間にある,ある実数。)これを発見者にちなんで[[w:テイラー級数]]と呼ぶ。これは複雑な関数をべき級数という比較的分かり易い関数で近似することが出来るということを表わす定理である。

=====テイラー展開の定義=====

上で述べたテイラー級数はn次までのべき級数によって展開したが、ある性質のいい関数については最後のややこしい項からの寄与が無限に小さくなり、単にその項をよりわかりやすい無限和で置き換えることが出来る。このときテイラー級数は
:<math>
f(x) = f(a) + f'(a)(x -a ) + \frac 1 2 f''(a) (x - a)^2 + \cdots + 
\frac 1 {n!} f ^{(n)} (a ) (x -a)^n + \cdots 
</math>
と書き換えられるが、これを[[w:テイラー展開]]と呼ぶ。テイラー展開は短く
:<math>
f(x) = \sum  _ {n=0} ^{\infty} \frac 1 {n!} f^{(n)} (a) (x-a) ^n
</math>
と書くことが出来る。

=====重要なテイラー展開の例=====

======eのx乗の例====== 
:<math>
f(x) = e^x
</math>
に対してx=0のまわりでのテイラー展開を導出する。
:<math>
f(x) = f'(x)=f''(x)= \cdots = f^{(n)}(x) = \cdots = e^x
</math>
であることを用いると、テイラー展開の定義の式で
:<math>
f^{(n)}(a)= f^{(n)}(0)=1
</math>
が得られる。
よって、<math>e^x</math>のx=0のまわりでのテイラー展開は、
:<math>
e^x = \sum  _n \frac 1 {n!} x^n
</math>
となる。

======1+xのa乗の例 (aは実数)======

:<math>
(1+x)^a
</math>
についてテイラー展開を考える。実際には、aが整数の場合にはこの値は通常のべき級数展開に一致する。例えば、
:<math>
(1+x)^2
</math>
をx=0のまわりでテイラー展開すると、
:<math>
f(0) = 1
</math>
:<math>
f'(0) = 2
</math>
:<math>
f''(0) = 2
</math> 
となり、
2次の代数式であるので3階以降の微分は0になることを考慮すると、そのテイラー展開は、
:<math>
\begin{matrix}
(1+x)^2 &= 1 + \frac 1 {1!} + 2 \cdot x + \frac 1 {2!} 2 x^2\\
&= 1 + 2x + x^2
\end{matrix}
</math>
となり、確かに通常の展開と一致する。

aが整数でない場合にはこの展開は無限に続く。この展開の係数をaが整数の場合の2項定理の拡張として、
:<math>
\begin{pmatrix}
a\\
n
\end{pmatrix}
</math>
と定義し、2項定数と呼ぶ場合がある。ここでaは<math>(1+x)^a</math>のaであり、nはxについてのn次の項を表わす。この係数を用いると、このテイラー展開は、
:<math>
(1+x)^a = \sum  _{n=0}^{\infty} 
\begin{pmatrix}
a\\
n
\end{pmatrix}
x^n
</math>
と書くことが出来る。例えば、a= 1/2では、x=0のまわりの展開について
:<math>
(1+x)^{1/2}
</math>
について、
:<math>
f(0) = 1
</math>
:<math>
f'(0) =\frac 1 2
</math>
:<math>
f''(0) = -\frac 1 4
</math>
が得られることから、2項目までのテイラー展開として、
:<math>
\begin{matrix}
(1+x)^{1/2}&= 1 + \frac 1 2 x + \frac 1 2 \cdot ( - \frac 1 4 ) x^2 + \cdots\\
&=  1 + \frac 1 2 x  - \frac 1 8  x^2 + \cdots
\end{matrix}
</math>
が得られる。もちろん根気があればどこまででも値を得ることが出来る。よって、
:<math>
\begin{pmatrix}
1/2\\
0
\end{pmatrix}
=1
</math>
:<math>
\begin{pmatrix}
1/2\\
1
\end{pmatrix}
=1/2
</math>
:<math>
\begin{pmatrix}
1/2\\
2
\end{pmatrix}
=-1/8
</math>
が得られる。


======sin x, cos xのテイラー展開======

<math>\sin x</math>と<math>\cos x</math>は微分によって互いに移り変わるのでそのテイラー展開は同時に扱うことが出来る。詳しく計算すると、x = 0のまわりでの展開は
:<math>
\sin x = \sum  _{n=1} ^{\infty} (-1)^n\frac 1 {(2n-1)!} x^{2n-1}
</math>
:<math>
\cos x = \sum  _{n=0} ^{\infty} (-1)^n\frac 1 {(2n)!} x^{2n}
</math>
を得ることが出来る。このとき、この値と、
:<math>
e^{ix}
</math>
のテイラー展開の値を比較した場合、
:<math>
e^{ix} = \cos x + i\sin x
</math>
の関係が示唆される。この関係は発見者の名にちなんで[[w:オイラーの公式]]と呼ばれる。この公式の正当化は複素関数論を使わないとうまくいかないようなのでこの稿の範囲を超えるが、[[物理数学II]]以降で扱われる予定である。<!-- しかし、今後はこの関係を断わり無く使うことがあるので注意して欲しい。
実際、この関係が無いと古典力学でいう単振動の方程式が
いきなり解けなくなるのでこの関係は非常に重要である。-->オイラーの公式を用いると、三角関数を指数関数を用いて表すことができる。具体的には、
:<math>
\cos x = \frac {e^{ix} + e^{-ix}}2, \sin x = \frac{e^{ix} - e^{-ix}}{2i}, \tan x =\frac 1 i  \frac{e^{ix} - e^{-ix}}{e^{ix} + e^{-ix}}
</math>
が成り立つ。
*問題例
*:上の結果を確かめよ。
*コーヒーブレイク
*:'''[[2のi乗の求め方]]'''
======テイラー展開を用いた関数の極限の計算======

テイラー展開を用いて極限を求めることが出来ることがある。例えば、
<math>x 
\rightarrow 0</math>で、
:<math>
\frac {e^x -1} x = \frac {x + x^2/2 + \cdots } x = 1 + x/2 \rightarrow 1
</math>
となる。

====ロピタルの定理====

aを実数または<math>\plusmn \infin</math>として

<math>
\lim _{x \rightarrow a}f(x) = \lim _{x \rightarrow a}g(x) = 0
</math>

または

<math>
\lim _{x \rightarrow a}f(x) = \lim _{x \rightarrow a}g(x) = \infin
</math>

となる微分可能な関数について

<math>
\lim  _{x \rightarrow a} \frac { f(x)}{g(x)} = \lim  _{x \rightarrow a} \frac { f'(x)}{g'(x)} 
</math>
例えば、
<math>
\lim _{x\rightarrow 0} \frac {\sin x} x = \lim _{x\rightarrow 0} \frac {\cos x} 1 = 1 
</math>
となる。

====積分の定義====

ある区間を考え、区間を細かく分割する。ここで、ある関数fに対して、分けられた区間でもっとも大きい部分をとり、区間の広さをかけて、足し合わせたものをその関数の上積分と呼ぶ。同様にもっとも小さい部分を取り足し合わせたものを関数の下積分と呼ぶ。上積分と下積分が一致するとき、それをその関数の積分と呼び、fを積分可能と呼ぶ。

Note:連続な関数は積分可能である。

*上積分と下積分が一致しない例
例えば関数
<math>
f (x) = 
\begin{cases}
1 ~ \textrm{(x is rational)} \\ 
0 ~ \textrm{(x is irrational)}
\end{cases}
</math>
について区間<math>0<x<1</math>で考えたとき、どんな小さい区間を使って<math>0<x<1</math>を分割したとしても有理数の稠密性により、上積分は1,下積分は0となる。よってfは積分可能でない。

====双曲線関数====
[[w:双曲線関数]]は三角関数と関係の深い一連の関数群である。これらは積分を行うための変数変換で使うことがあるので、ここで導入する。双曲線関数は次の3つの関数である。
:<math>
\sinh x = \frac {e^x - e^{-x}} 2
</math> 
:<math>
\cosh x = \frac {e^x + e^{-x}} 2
</math> 
:<math>
\tanh x = \frac {\sinh x} {\cosh x}
</math>
を双曲線関数と呼ぶ。
:それぞれのグラフ
これらは関係式
:<math>
\cosh ^2 x - \sinh ^2 x = 1
</math>
を満たすが、<math>x^2 -y^2 = 1</math>が双曲線の関数表示であることから、この関数は双曲線関数と呼ばれる。上の式は三角関数の対応物である<math>
\cos ^2 x + \sin ^2 x = 1</math>に類似しているが、この結果は偶然ではない。上のオイラー公式を使った三角関数の式を見ると、
<!--
更に
<math>
e ^{ix} = \cos x + i \sin x
</math>
を用いると、
 <math>e^{x}</math> から解析接続? -->
<math>
\sin ix = i\sinh x, \cos ix = \cosh x, \tan ix = i\tanh x
</math>
が得られる。この式を<math>\cos ^2 x + \sin ^2 x = 1</math>でx=izとしたものに代入すると、<math>\cosh ^2 x - \sinh ^2 x = 1</math>の関係が得られる。

====三角関数の逆関数====

<math>
\sin^{-1} x 
</math>
を<math>\sin x</math>の逆関数とする。これは多価関数であるので通常
<math>
-\pi  < y <  \pi
</math>
の範囲を選んで用いる。
<!--  複素解析によると無限多価は\ln から来る。 -->
同様に
<math>
\tan ^{-1} x 
</math>
も
<math>
-\pi  < y <  \pi
</math>
の範囲を選んで用いる。一方
<math>
\cos ^{-1} x 
</math>
は
<math>
0  < y <  2\pi
</math>
の範囲を選んで用いる。

====様々な積分====

:<math>
\int dx \frac 1 {\sqrt{1-x^2}}  = \textrm{Arcsin} x
</math> 
:<math>
\int dx \frac 1 { 1+x^2} = \textrm{Arctan} x
</math>
が得られる。

*導出

まず、
:<math>
\int dx \frac 1 {\sqrt{1-x^2}}  = \textrm{Arcsin} x
</math> 
を導出する。<math>y= \sin x</math>とする。このとき、
:<math>
\frac{d{y}}{d{x}} = \cos x
</math>
:<math>
\begin{matrix}
\frac{d{x}}{d{y}} &= \frac 1 {\cos x }\\
&= \frac 1 {\sqrt { 1 - (\sin x)^2} }\\
&= \frac 1 {\sqrt { 1 - y^2} }\\
\end{matrix}
</math>
よって、
:<math>
x = \sin ^{-1} y
</math>
と合わせると、
:<math>
\frac{d{x}}{d{y}} = \frac{d(\sin ^{-1} y)}{dy} = \frac 1 {\sqrt { 1 - y^2} }
</math>
となり、2番目の式と、3番目の式をyで積分することで求める式を得る。

<math>y = \tan x </math>とおく。
:<math>
\frac{d{y}}{d{x}} = \frac 1 {\cos ^2 x}
</math>
より、
:<math>
\begin{matrix}
\frac{d{x}}{d{y}} &= \cos ^2 x\\
&= \frac 1 {1+\tan^2 x } \\
&= \frac 1 {1+y^2 } \\
\end{matrix}
</math>
となる。よって、
:<math>
\frac{d{x}}{d{y}} = \frac{d(\tan^{-1}y)}{dy} = \frac 1 {1+y^2}
</math>
が得られた。この式の2番目の式と、3番目の式をyで積分することで求める式を得る。

====有理関数の積分====
[[w:有理関数]]の積分

有理関数は必ず[[w:初等関数]]を用いて積分できる。

<!-- (
*TODO
初等関数の定義
) -->

*説明

有理関数の積分は
:<math>
\int dx \frac {P(x)}{Q(x)}
</math>
の形に書くことが出来る。(P,Qはxの整式。)ここで、次のような手順を実行する。

**上の式の次数が下の式の次数より高かったら、上の式を下の式で割る。
このことによって、被積分関数の分母の次数は、上の式の分子の次数より低くなる。割ることであまった部分は必ず、分数でない形になるので(普通の数やx,<math>x^2</math>などになる。)積分できる。

**次に、分母を因数分解する。
代数式は必ず複素数の範囲で因数分解できることが知られているので、([[w:代数学の基本定理]])
分母は必ず(x-a)の積の形に書ける。ここで、元々の被積分関数が実数だったとすると、<!-- (この場合はそのようにしている。複素数だった場合は物理数学IIの範囲となる。) -->因数分解された式は、必ず、<math>(x-a)(x-a^*)</math>の形になっているはずである。(*は複素共役)これらの2因数をかけ合わせることにすると、結局これらの式の分母は、1次式か2次式の積で書ける。

**次に、得られた分母を使って[[w:部分分数分解]]を行なう。例えば、
:<math>
\frac 1 {(x^2-1) }
</math>
については、
:<math>
\begin{matrix}
=&\frac 1 {(x-1)(x+1) } \\
=&\frac {1} 2 (\frac 1 {x-1 }-  \frac 1 {x+1})
\end{matrix}
</math>
が得られる。

**これまでの手順を経た結果、得られた式の形は、
:<math>
\int \frac 1 {x-a}
</math>
か、
:<math>
\int \frac {x-b} {cx^2 + dx +e}
</math>
が得られることが分かる。これらは共に初等関数の範囲で積分可能である。実際、上の式は
:<math>
\int \frac 1 {x-a}
</math>
:<math>
= \ln |x-a|
</math>
を満たすことが分かる。下の式については、まず、分母を平方完成すると、分母は、
:<math>
c(x-\alpha) ^2 + \beta
</math>
の形になるが、ここで
:<math>
y=x-\alpha
</math>
の置き換えをすると、元々の積分は、
:<math>
\int \frac {y+f} {y^2 +g}
</math>
となる。ここで、このうちの第1項は、
:<math>
\begin{matrix}
\int \frac y {y^2+g} dy\\
=& \int \frac 1 2 \frac 1 {y^2+g} dy^2\\
=&  \frac 1 2 \ln (y^2 +g ) \\
\end{matrix}
</math>
が得られ、積分できることが分かる。次に、第2項については
:<math>
z = \sqrt {g} y
</math>
の置き換えをすると、定数因子を除いて、
<!-- (
*TODO
定数因子の計算
) -->
:<math>
\int \frac {1} {z^2 +1}
</math>
となるが、この積分の結果はこのページの上の方で見た通り、
:<math>
= \tan^{-1} z
</math>
となる。

よって、全ての有理関数は、初等関数の範囲で積分できることが分かった。


*計算例
計算例として、
:<math>
\int dx \frac {x+3} {(x+2) (x^2+1) }
</math>
を実際に計算してみる。
<!-- あらかじめ、[[Maxima]]を使って、計算してみた。
 
 integrate((x+3)/((x^2+1)*(x+2)),x);

答は、
 
                                        3 + x
 (%i2)                   INTEGRATE(----------------, x)
                                         2
                                   (1 + x ) (2 + x)
                            2
                       LOG(x  + 1)   LOG(x + 2)   7 ATAN(x)
 (%o2)               - ----------- + ---------- + ---------
                           10            5            5

読みづらいかも知れないが、さしあたり計算値はわかるものと期待する。
-->
計算を行なうときにはまず、分子の次数が分母の次数よりも低いことを確認する。次に、部分分数分解を行なうが、このときには、
:<math>
\frac A {x+2} + \frac {Bx+C} { x^2+1} 
</math>
とおいて計算すればよい。ここで、分母を通分すると、分子は、
:<math>
(A+B) x^2 + (2B+C) x + A+2C
</math>
が得られるが、これは元々の式の分子である
:<math>
x+3
</math>
と一致していなくてはならない。よって、
:<math>
\begin{matrix}
A+B = 0\\
A+2C = 3\\
2B+C = 1
\end{matrix}
</math>
が得られる。これを解くと、
:<math>
A= \frac 1 5, B = - \frac 1 5, C = \frac 7 5
</math>
が得られる。元の積分は
:<math>
\int dx (\frac 1 5 \cdot \frac 1 {x+2} + \frac {\frac{-1}5 x + \frac 7 5  }{x^2 +1} )
</math>
に帰着するが、これらの項ははそれぞれ初等関数の範囲で積分できる。実際に積分を行なうと、
:<math>
\int dx (\frac 1 5 \cdot \frac 1 {x+2} + \frac {\frac{-1}5 x + \frac 7 5  }{x^2 +1} )
</math>
:<math>
= \frac 1 5 \ln (x+2)  -\frac 1 {10} \ln (x^2+1) + \frac 7 5 \tan^{-1} x
</math>
が得られ、上で得た値と一致する。

====無理数を含んだ積分====

関数が有理数だけで書かれない場合、その式はもはや積分が出来るとは限らない。簡単に積分が実行できる場合を挙げる。すぐに積分の仕方が見当たらない場合、それが定積分であったら、数値的に求めることを考えることも必要である。
*<math>
\sqrt{1-x^2}
</math> の場合。
:<math>
\sqrt{1-x^2}
</math>
で書かれる無理式が含まれ、それ以外の無理式が含まれない場合には、
:<math>
x =  \sin t ~\textrm{or}~ x = \cos t
</math>
の置き換えをすることで、この式を三角関数の積分に置き換えることが出来る。三角関数の積分は、後に述べる通り有理関数の積分に帰着させることが出来るので、この積分は解析的に実行できる。

*<math>
\sqrt{1+x^2}
</math> の場合。
:<math>
\sqrt{1+x^2}
</math>
で書かれる無理式が含まれ、それ以外の無理式が含まれない場合には、
:<math>
x = \sinh t
</math>
の置き換えをすることで、この式を三角関数の積分に置き換えることが出来る。
(
:<math>
\cosh ^2 x  = 1+ \sinh^2 x
</math>
の関係を用いて、根号を消すことが出来る。
)

====三角関数を含んだ積分====

三角関数
:<math>
\sin x, \cos x, \tan x
</math> 
だけを含んだ積分については、
:<math>
\tan (x/2) = t
</math>
の置き換えをすることで、これを有理関数の積分に帰着させることができる。実際、
:<math>
\begin{matrix}
\sin x =& 2 \sin \frac x 2 \cos \frac x 2\\
=& 2 \tan \frac x 2 \cos^2 \frac x 2\\
=&  \frac {2 \tan \frac x 2} {1+\tan ^2 \frac x 2}\\
=& \frac {2 t} {1+t^2} 
\end{matrix}
</math>

:<math>
\begin{matrix}
\cos x &= 2 \cos^2 \frac x 2 - 1\\
&= 2 \frac 1 {1+t^2} - 1\\
&= 2 \frac {1-t^2} {1+t^2} 
\end{matrix}
</math>

:<math>
\begin{matrix}
\tan x &= \frac {\sin x} {\cos x }\\
&= \frac {2t} {1-t^2}
\end{matrix}
</math>
さらに、
:<math>
dx = d(2\tan^{-1} t) = 2 \frac {dt} {1+t^2}
</math>
となり、確かにtについての有理関数に帰着することが分かる。よって、三角関数だけの関数は初等関数の範囲で積分され得ることが分かった。

<!--
*TODO 周期が<math>\pi</math>であるとき
-->

*問題例

**問題
(I)
:<math>{{x+2}\over{x^2+5\,x+4}}</math>
(II)
:<math>{{x^2+3\,x+3}\over{2\,x^3-x+1}}</math>
(III)
:<math>{{\sin x}\over{\cos ^2x}}</math>
(IV)
:<math>
\int dx \frac 1 {1+\sin x}
</math>
をそれぞれ積分せよ。

**解答
(I)
:<math>
\frac {x+2} {x^2+5x+4}
</math>
:<math>
=\frac {x+2} {(x+4)(x+1)}
</math>
ここで、この式が
:<math>
\frac A{x+4} + \frac B {x+1}
</math>
に等しいとすると、
:<math>
\frac {x+2} {(x+4)(x+1)}
=\frac 1 {(x+4)(x+1)} (A(x+1) + B(x+4))
</math>
両辺が等しいことから、
:<math>
A + B = 1, A+4B = 2
</math>
となり、
:<math>
B = \frac13, A = \frac23
</math>
が得られる。よって始めの式について、
:<math>
\frac {x+2} {x^2+5x+4}
</math>
:<math>
= \frac 13 \frac 1 {x+1} + \frac 23 \frac 1{x+4}
</math>
が得られる。この関数をxで積分すると
:<math>
\frac 13 \ln (x+1) + \frac 23 \ln (x+4)
</math>
が得られる。
(II)
:<math>
\frac {x^2 + 3x+3}{2x^3 - x +1}
</math>
:<math>
=\frac{x^2 + 3x+3}{(x+1)(2x^2 -2x +1)}
</math>
ここで、この式が
:<math>
\frac A {x+1} + \frac {Bx +C }{2x^2 -2x + 1}
</math>
に等しいと仮定すると、両辺の分母を比較することで、
:<math>
x^2 + 3x +3 = A(2x^2 -2x +1)+ (Bx+c)(x+1)
</math>
となり、
:<math>
2A+B = 1, -2A + B + C = 3, A+C = 3
</math>
が得られる。これを解くと、
:<math>
A= \frac 15, B = \frac 3 5, C = \frac {14}5
</math>
が得られる。よって元の式は、
:<math>
\frac 15 \frac 1 {x+1}+ \frac {\frac 3 5x + \frac {14}5}{2x^2-2x+1}
</math>
となる。更にこの式の第2項について、項の分子を
:<math>
\frac 35(2x^2-2x+1)' \frac 14 + \frac 3 {20} \cdot 2 + \frac {14}5
</math>
:<math>
= \frac 3{20}(2x^2-2x+1)' + \frac {31} {10}
</math>
と書き換えられる事に注目すると、元の式は
:<math>
\frac 15 \frac 1 {x+1}+ \frac {\frac 3 {20}(2x^2-2x+1)' }{2x^2-2x+1}
+ \frac {31}{10} \frac 1{2x^2-2x+1}
</math>
となる。ここで、この式の1, 2項については、簡単に積分できて、
:<math>
\frac 15 \ln(x+1) + \frac 3{20} \ln (2x^2-2x+1)
</math>
が得られる。最後に第3項については、
:<math>
\frac 1{2x^2-2x+1} = \frac 2 {(2x-1)^2 + 1}
</math>
が成り立つことに注目すると、<math>t = 2x -1, dt = 2dx</math>の置き換えをして、
:<math>
\int dx 
\frac 2 {(2x-1)^2 + 1}
</math>
:<math>
= \int dt \frac 1{t^2 +1}
</math>
:<math>
= \textrm{Arctan} t = \textrm{Arctan} (2x-1)
</math>
が得られる。よって、全体をまとめると積分値として
:<math>
\frac 15 \ln(x+1) + \frac 3{20} \ln (2x^2-2x+1)+\frac{31}{10}\textrm{Arctan} (2x-1)
</math>
が得られる。

(III)
:<math>
\int \frac {\sin x}{\cos^2 x} dx
</math>
:<math>
= - \int 
\frac {1}{\cos^2 x} d(\cos x)
</math>
:<math>
= \frac 1 {\cos x}
</math>

(IV)
:<math>
t = \tan (\frac x 2)
</math>
としたとき、
:<math>
dx = \frac {2dt}{1+t^2}, \sin x = \frac {2t}{1+t^2}
</math>
となることを考慮すると、
:<math>
\int dx \frac 1 {1+\sin x}
</math>
:<math>
= \int \frac 1 {1+ \frac {2t}{1+t^2}}\frac {2dt}{1+t^2}
</math>
:<math>
= \int dt \frac 2 {(1+t)^2} 
</math>
:<math>
= - \frac 2 {1+t} = \frac {-2}{1+\tan (\frac x 2)}
</math>
となる。別の方法として、<!-- 結果から逆算したのだが ... 。 -->
:<math>
\frac 1 {1+\sin x} = \frac 1 {1+2\sin \frac x 2 \cos \frac x 2}
</math>
:<math>
= \frac 1 {(\sin \frac x 2 + \cos \frac x 2)^2}
</math>
:<math>
= \frac 1{1+\tan \frac x 2}\frac 1 {\cos ^2 \frac x 2}
</math>
:<math>
= \frac d {dx}(\frac 1 {1 + \tan \frac x 2}) (-2)
</math>
となるので、両辺を積分して結果を得てもよい。

===多変数関数の微積分===

====偏微分====

多変数で定義された関数fがあるときのある変数のみを対象にした微分、例えば
<math>
\lim_{h\rightarrow 0} \frac {f(x _1 + h, ... ,x _n) - f(x _1, ...,x _n) } h
</math>
を
<math>
f_x
</math>
や
<math>
\frac{\partial{f}}{\partial{{x _1}}}
</math>
、
<math>
(\frac{\partial{f}}{\partial{{x _1}}})_{x_2, x_3...}
</math>
と書き{{Ruby|偏微分|へんびぶん}}と呼ぶ。

====多変数関数の最大最小値====

====偏微分による計算====
多変数関数ではあらゆる独立変数による偏微分がすべて0になる点で、関数が最大値または最小値を取ることが期待される。

例えば
<math>
f = x^2 + y^2
</math>
では、
<math>
\frac{\partial{f}}{\partial{x}} = 2x
</math>
<math>
\frac{\partial{f}}{\partial{y}} = 2y
</math>
であるので、
<math>
x = 0 , y = 0
</math>
で、極大値または極小値を取ることが期待される。

====2変数関数の極値====
2変数関数<math>f(x,y)</math>において、点<math>(a,b)</math>で<math>f_x(a,b)=f_y(a,b)=0</math>とする。判別式'''D'''を
:<math>D=f_{xx}(a,b)f_{yy}(a,b)-f_{xy}(a,b)^2</math>

と定義する。

<math>D>0</math>のとき
:<math>f_{xx}(a,b)>0</math>ならば、関数<math>f(x,y)</math>は点<math>(a,b)</math>で極小値をとり
:<math>f_{xx}(a,b)<0</math>ならば、関数<math>f(x,y)</math>は点<math>(a,b)</math>で極大値をとる。

<math>D<0</math>のときは、極値はとらない。

====全微分====
2変数関数<math>f(x,y)</math>における全微分は
:<math>df=\frac{\partial{f}}{\partial{x}}dx+\frac{\partial{f}}{\partial{y}} dy </math>
と定義される。例えば<math>f(x,y)=x^2+y^2</math>における全微分は
:<math>df=2xdx+2ydy</math>
となる。同様に''n''変数関数<math>f(x_1,x_2,\cdots,x_n)</math>における全微分は
:<math>df=\frac{\partial{f}}{\partial{x_1}}dx_1+\frac{\partial{f}}{\partial{x_2}} dx_2 + \cdots + \frac{\partial{f}}{\partial{x_n}}dx_n</math>
と定義される。

====ヘッセ行列====

ヘッセ行列は2階偏微分によって作られた行列
<math>
H = \left[
\frac{\partial^2{f}}{\partial x_ix_j}(P)
\right]
</math>である。

点Pを、
<math>
\frac{\partial f}{\partial x_1}(P)=
\frac{\partial f}{\partial x_2}(P)=
\cdots
\frac{\partial f}{\partial x_n}(P)=0
</math>
なる点([[w:臨界点]])とする。ヘッセ行列のPにおける固有値が全て正であれば、関数は点Pで極小値を持ち、全て負であれば、点Pで極大値を持つ。どちらでもないなら点Pは[[w:鞍点]]である。

例えば、
<math>
f = x^2 + y^2
</math>
について、臨界点(0,0)におけるヘッセ行列は、
<math>
H = 
\begin{pmatrix}
\frac{\partial^2{f}}{\partial{x}^2} & \frac {\partial^2 f}{\partial x\partial y}\\
\frac {\partial^2 f}{\partial y\partial x} & \frac{\partial^2{f}}{\partial{y}^2} 
\end{pmatrix}
</math>
<math>
=
\begin{pmatrix}
2 & 0\\
0&2
\end{pmatrix}
</math>
となる。固有値は2であるので、点(0,0)はfの極小値である。

====陰関数定理====


:<math>
F(x,y) = 0
</math>
の形で表わされる関数があるとき、
:<math>
\frac{\partial{F}}{\partial{y}}
</math>
が存在したとすると、この関数は
:<math>
y = y(x)
</math>
の形に(局所的には)表わすことが出来る。このとき、
:<math>
y' = - \frac { \frac{\partial{F}}{\partial{x}}}{\frac{\partial{F}}{\partial{y}}}
</math>
が成り立つ。

右辺の形は少し妙に見えるかも知れないが例えば、
:<math>
F(x,y) = ax+by
</math>
(a,bは定数)について考えてみると、上の式は、
:<math>
y' = - \frac a b
</math>
となっており、通常の仕方で見たyの傾きと一致している。

この定理は結局のところどんな複雑な曲線でも、ある点のすぐ近くに限れば、それはほとんど直線と同じ様になっているということを述べている。

====Lagrangeの未定定数法====

F(x,y) = 0の形の条件が課せられた中で、
:<math>
z = f(x,y) 
</math>
の最大値を求める問題を考える。このとき
:<math>
g = f + \lambda F
</math>
で新しい関数gを定義し、
(<math>\lambda</math>はある定数)
:<math>
\frac{\partial{g}}{\partial{x}} = \frac{\partial{g}}{\partial{y}} = \frac{\partial{g}}{\partial{{\lambda}}} = 0
</math>
で与えられる<math>x,y,\lambda</math>を計算する。得られた点が極大か極小値を取る点である。

*計算例
:<math>
z = x^2+y^2, F= x+y-1
</math>
として、この方法を適用してみる。極値は、(図を書いてみると)
:<math>
x= y = \frac 1 2
</math>
で現われると期待される。
<!-- (実際にはグラフ中の最小値。?) -->
この式の場合は、
:<math>
y= -x +1 
</math>
を代入することで答を得ることもできる。平方完成した形は
:<math>
z = 2(x-\frac 1 2)^2 + \frac 1 2
</math>
であり、確かに
:<math>
x= y = \frac 1 2
</math>
で極値を取ることが分かる。未定定数法を用いると

:<math>
g= f +\lambda F = x^2+y^2 + \lambda (x+y-1)
</math>
が得られる。
ここで、
:<math>
\frac{\partial{g}}{\partial{\lambda}} = x+y-1 = 0
</math>
:<math>
\frac{\partial{g}}{\partial{x}} = 2x +\lambda = 0
</math>
:<math>
\frac{\partial{g}}{\partial{y}} = 2y +\lambda = 0
</math>
が得られるが、これはx,y,<math>\lambda</math>についての連立1次方程式となっている。これを解くと、答は、
:<math>
\lambda = 1, x = y = \frac 1 2
</math>
となり、確かに正確な値と一致する。

====多重積分====

複数の文字について積分を行なうときこれを多重積分と呼ぶ。例えば、
<math>
\iint f(x,y) dx dy 
</math>

====累次積分====

<math>
\iint f(x,y) dx dy 
</math>
は、
<math>
\iint f(x,y) dx dy 
</math>
<math> 
= \int dy (\int f dx ) = \int dx (\int f dy ) 
</math>
で書き変えられる。

====ガウス積分====
ガウス積分<math>\int_{-\infty}^{\infty} e^{-x^2}dx=\sqrt{\pi}</math>の導出。

:<math>I=\int_{-\infty}^{\infty} e^{-x^2}dx=\int_{-\infty}^{\infty} e^{-y^2}dy</math>とおくと、

:<math>I^2=\int_{-\infty}^{\infty} e^{-x^2}dx \int_{-\infty}^{\infty} e^{-y^2}dy</math>
::<math>=\iint_{D} e^{-x^2-y^2}dxdy</math>（<math>D=\left\{(x,y)|-\infty<x<\infty,-\infty<y<\infty \right\}</math>）
::<math>=\iint_{D'} e^{-r^2} r dr d\theta</math> （直交座標から極座標に変換、<math>x=r\cos \theta,y=r\sin\theta,dxdy\to rdrd\theta</math>、<math>D'=\left\{(r,\theta)|0\le r <\infty,0\le\theta\le 2\pi\right\}</math>）

::<math>=\int_{0}^{\infty}re^{-r^2}dr \int_{0}^{2\pi}d\theta</math>
::<math>=\frac{1}{2}\times 2\pi</math> 
::<math>=\pi</math>
:<math>I=\int_{-\infty}^{\infty} e^{-x^2}dx=\sqrt{\pi}</math>

また:<math>\int_{-\infty}^{\infty} e^{-ax^2}dx</math>の積分は
:<math>\int_{-\infty}^{\infty} e^{-ax^2}dx=\int_{-\infty}^{\infty} e^{-t^2}\frac{1}{\sqrt{a}}dt</math>（<math>\sqrt{a}x=t</math>と置いて置換積分）
::<math>=\frac{1}{\sqrt{a}}\int_{-\infty}^{\infty} e^{-t^2}dt</math>
::<math>=\sqrt{\frac{\pi}{a}}</math>

となる。

====ガンマ関数====
ガンマ関数は<math>\Gamma(t)=\int_{0}^{\infty}x^{t-1}e^{-x}dx </math>で定義される関数である。

;ガンマ関数の基本公式
:<math>\Gamma(t+1)=t\Gamma(t)</math>
:<math>\Gamma(n+1)=n!</math>（''n''は自然数）
:<math>\Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}</math>

====ベータ関数====
ベータ関数は<math>\Beta(p,q)=\int_{0}^{1} x^{p-1}(1-x)^{q-1}dx</math>で定義される関数である。

;ベータ関数の基本公式
:<math>\Beta(q,p)=\Beta(p,q)</math>
:<math>\Beta(p+1,q)=\frac{p}{q}\Beta(p,q+1)</math>
:<math>\Beta(p,q+1)=\frac{q}{p}\Beta(p+1,q)</math>
:<math>\Beta(p+1,q)=\frac{p}{p+q}\Beta(p,q)</math>
:<math>\Beta(p,q+1)=\frac{q}{p+q}\Beta(p,q)</math>
:<math>\Beta(p,q)=2\int_{0}^{\frac{\pi}{2}} \sin^{2p-1} \theta \cos^{2q-1} \theta d\theta</math>
:<math>\int_{0}^{\frac{\pi}{2}} \sin^a \theta \cos^b \theta d\theta = \frac{1}{2}\Beta\left(\frac{a+1}{2},\frac{b+1}{2}\right)</math>
:<math>\Beta(p,q)=\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)}</math>

===数列の収束===

====ダランベールの収束判定法====
<!--
<math>
 \frac {a _{n}}{ a _{n-1} } < 1 
</math>

を満たす数列は収束する。(ダランベールの収束判定法)
(?)
-->



====Σ{1/(n^α)}の収束発散====

<math>
\sum _{n=1}^{\infty} \frac 1 {n^{\alpha}}
</math>
は、<math>\alpha =< 1</math>のとき発散し、<math>\alpha > 1</math>のとき収束する。


TODO
* 実数の連続性 
* 複雑な積分 
* 数列の収束
* 偏微分は交換可能
[[カテゴリ:解析学|ふつりすうかくいち かいせきかく]]