数理解析研究所講究録
第 1780 巻 2012 年 55-63

55

数式処理の発想を利用した数学教育の試み

四ツ谷 晶二 (Shoji Yotsutani)

龍谷大学理工学部

Fuculty of Science and Technology,Ryukoku University

1

はじめに2009 年度のこの研究集会において，数式処理システムを用いた，完全楕円積分の$\sqrt{1-k^{2}}$ と
商 $E/K$ の上下からの比較関数の構成法について報告をした．
$E(k)/K(k)$ のマクロリン展開を組み合わせた上下からの比較関数を構成するものであった．しかし，$\frac{E(k)}{K(k)}\sim\frac{1}{\log(4(1-k^{2})^{-1/2})}$

as $karrow 1-0$

$k=1$ での特異性のため応用上は不十分であった．
であり，前半部分では，まず，相加・相乗平均を用いたガウスによる高精度の近似計算法を利用した，分数幕 $(1-k^{2})^{1/2^{n}}(n=1,2, \cdots)$ を用いた比較関数を紹介する．
後半部分では，多項式の零点の個数を決定する，Sturm の定理のアイデアをわかりやすく説明するひとつの方法を紹介する．その方法は，高校の 2 次関数等の教育において別の観点や指導法を与えることを可能とする．

2

完全楕円積分の商の近似について

Cahn-Hilliard 方程式等の解の大域的分岐構造を調べる際，第 1 種完全楕円積分 $K(k)$ ,第 2 種完全楕円積分 $E(k)$ , および からなる有理式の増減，値域を調べる問題に出会う([2]). 応用上，特に K, E に関する斉次式の場合が重要である．この場合，$E/K$ を上下$k$から挟むことにより，2 変数多項式に関する種々の標準的な手法が利用でき解析可能となる．ここで，第 1 種，第 2 種完全楕円積分は次式で定義される．
$K(k):=\int_{0}^{\pi/2}\frac{d\varphi}{\sqrt{1-k^{2}\sin^{2}\varphi}}$,$E(k):=\int_{0}^{\pi/2}\sqrt{1-k^{2}\sin^{2}\varphi}d\varphi$

先に述べたように，2009 年度のこの研究集会おいて，$\sqrt 1=k^{2}$.と $E/K$ のマクロリン$k=1$ での特異性のため応用上展開を組み合わせた上下からの比較関数を報告したが，は不十分であった．

相加相乗平均を用いたガウスによる，古典的な高精度近似計算法を利用した分数幕$(1-k^{2})^{1/2^{n}}(n=1,2, \cdots)$による比較関数を紹介する．

$E(k)/K(k)arrow 0(karrow 1-O)$ であるから，$E(k)/K(k)|_{k=1}:=0$ と定義すると，まず，$E(k)/K(k)$ は閉区間 $[0,1]$ で連続な関数となる．

56
任意の $a>b>0$ に対して，算術幾何平均の漸化式$a_{l+1}= \frac{a_{\ell}+b_{\ell}}{2}$$a_{0}=a,$,$b_{0}=b,$,$b_{\ell+1}=\sqrt{a_{\ell}b_{\ell}}$$c_{l+1}= \frac{a_{\ell}-b_{\ell}}{2}(\ell=0,1,2,3, \cdots)$.

$c_{0}=\sqrt{a^{2}-b^{2}}$を考える．極限値を AGM(a,b) とかく．1818 年，Gauss は次を得た ([3],[1]).
Theorem A. $a=1,$ $b=\sqrt{1-h}$ とする．このとき，$\frac{E(\sqrt{h})}{K(\sqrt{h})}=1-\sum_{\ell=0}^{\infty}2^{\ell-1}c_{\ell}^{2}$.

この定理は，AGM(a,b) と完全楕円積分との関係を調べるなかで発見されたものである．右辺は の関数である．これより，上からの比較関数はごく自然に，下からの比較関数は若干の補正で得られる ([4]).
$h$Theorem. 関数を$\underline{g}_{n}(h):=1-\sum_{\ell=0}^{n}2^{\ell-1}c_{\ell}^{2}-2^{n-1}c_{n}^{2}$,$\overline{g}_{n}(h)_{\backslash }=1-\sum_{\ell=0}^{n}2^{\ell-1}c_{\ell}^{2}$で定義する．このとき，次の事実が成立する．

1 任意の非負の整数$n$$)$

に対して，次の評価式が成り立つ:$\underline{g}_{n}(h)\leq E(\sqrt \mathcal{H})/K($而$h\in[0,1]$$)\leq\overline{g}_{n}(h)$.

なお， ん との等号成立は $h=0,1$ の場合のみであり， ん との等号成立は $h=0$ の$\underline{g}_{n}($$\overline{g}_{n}($$)$$)$場合のみである．

2$)$$\underline{g}_{n}(h),$$\overline{g}_{n}(h)$具体的に，$\supset$$\underline{g}_{n}(h),$$E(\sqrt{h})/K(\sqrt{h})$$\overline{g}_{n}(h)(n=1,2,3)$$uni$

formly on $[0,1]$

as$narrow\infty$を示す．

,$\underline{g}_{0}(h)$$=$$1-h$ ,$\overline{g}_{0}(h)$$=$$1- \frac{h}{2}$$\underline{g}_{1}(h)$$=$$(1-h)^{1/2}$ ,$\overline{g}_{1}(h)$$=$$\frac{1}{2}-\frac{h}{4}+\frac{(1-h)^{1/2}}{2}$$\underline{g}_{2}(h)$$=$$(1-h)^{1/4}+(1-h)^{3/4}$$\overline{g}_{2}(h)$$=$$-(1-h)^{1/2}$ ,.

,$\frac{1}{4}-\frac{h}{8}-\frac{(1-h)^{1/2}}{4}$$+ \frac{(1-h)^{1/4}}{2}+\frac{(1-h)^{3/4}}{2}$.

関数 $1/K(k)$ にも，同様の高精度な近似関数が構成されることがわかる．

上記の方法を，数式処理システムと組合わせれば，いくらでも高精度に近似式を自動的に生成することができる．

57

3Sturm の定理とは話題を全くかえて，多項式の零点の個数を決定する，Sturm の定理とその証明のアイデアをわかりやすく説明し，教育に利用する試みについて説明したい．
Sturm の定理を説明するためには Sturm 列が不可欠である．このため，ユークリッドの互除法を思い出す．互除法は 2 つの自然数の最大公約数を簡単かつ高速にもとめるアルゴリズムである．互除法は多項式の共通因子のうちで最大の次数のものを求める問題にもそのまま適用することができる．
例として，多項式$f(x):=x^{3}+x^{2}-3x+1$を考える．$f(x)=0$ が重解をもつかどうを判定しよう．この判定のためには，$f(x)$ と$f’(x):=3x^{2}+2x-3$の共通因子を求めればよい．これらに対してユークリッドの互除法を適用すると$f(x)$$=$$( \frac{x}{3}+\frac{1}{9})f’(x)+(-\frac{20}{9}x+\frac{4}{3})$$f’(x)$$=$$(- \frac{27}{16}x-\frac{171}{100})(-\frac{20}{9}x+\frac{4}{3})+(-\frac{18}{25})$を得る．よって，最大共通因子は -2763/256 である．したがって，$f(x)=0$ と $f’(x)=0$$f(x)=0$ は重解を持たない．
は共通解を持たない．ゆえに，ユークリッドの互除法の計算において，余りの項をわざわざマイナスにして得られる関数列を Sturm 列という．例えば，上記の場合では$f(x)$$=$$( \frac{x}{3}+\frac{1}{9})f’(x)-(\frac{20}{9}x-\frac{4}{3})$$f’(x)$$=$$( \frac{27}{16}x+\frac{171}{100})(\frac{20}{9}x-\frac{4}{3})-\frac{18}{25}$と書き直すことができるので，Sturm 列は$f(x),$ $f’(x),$$\frac{20}{9}x$$\frac{4}{3},$$\frac{18}{25}$である．

スツルムの定理 実係数の代数方程式 $f(x)=0$ が重解をもたずとする．このとき区間 $(a, b)$ における零点の個数は$f(a)\neq 0$ かっ$f(b)\neq 0$$V(a)-V(b)$である．ただし，$V(x)$ は符号変化数である．なお，符号変化の回数を数えるとき，途中にO が現れた場合はないものと思うことにする

58
スツルムの定理とは，実係数をもつ方程式が与えられた区間 $(a, b)$ の中にいくつの実根をもつかを決定するためのものである．多項式 $f(x)$ とその導関数 $f’(x)$ にユークリッ$f(x)$ を最大公約因子で割ることによって重解ドの互除法を行い最大公約因子を求め，$f(a)\neq 0$ かつ $f(b)\neq 0$ としていい．もしとなる部分は取り除く事が出来る．更に，$f(x)$ をそれで割って$f(a)=0,$ $f(b)=0$ であれば，もともと因数 $x-a,$ $x-b$ をもつので，おけばいい．そのため多項式 $f(x)$ とその導関数 $f’(x)$ は互いに素であり，方程式 $f(x)\neq 0$を重解をもたないという状況で，さらに $f(a)\neq 0$ かつ $f(b)\neq$ であるとする．このとき前項で説明したユークリッドの互除法を実行すると最後の余は 0 ではない定数である．
例区間 $(-2,2)$ における，方程式$x^{3}+x^{2}-3x+1=0$の実数解の個数を求めてみよう．

まず，$f(x):=x^{3}+x^{2}-3x+1$とおく．前例より方程式 $f(x)=0$ は重解をもたない．また$f(-2)=3\neq 0$ ,$f(2)=7\neq 0$である．上でみたように，スツルム列は $\{f_{0}(x), fi(x), f_{2}(x), f_{3}(x), \}$ である．ただし，$f_{1}(x)=f’(x)=3x^{2}+2x-3$ ,$f_{0}(x)$$=$$f(x)$ ,$f_{2}(x)$$=$$\frac{20}{9}x-\frac{4}{3}$,$f_{3}(x)= \frac{18}{25}$.

前節で求めたスツルム列に $x=\pm 2$ を代入する．
$f_{0}(-2)$$=$

3, $f_{1}(-2)=f’(x)=5,$ $f_{2}(-2)=- \frac{52}{9},$$f_{3}(-2)= \frac{18}{25}$$f_{0}(2)$$=$7, $f_{1}(2)=f’(x)=13$ ,$f_{3}(2)= \frac{18}{25}$$f_{2}(2)=- \frac{28}{9}$,となる．よってとなり$V(-2)=2,$ $V(2)=0$である．したがって，$V(-2)-V(2)=2$を得る．以上より，区間 $(-2,2)$ における実数解の個数は 2 個である．

,

59

4

Sturm の定理の証明のアイデア
これまで，多くの本に Sturm の定理の証明が紹介されているが，抽象的に書かれていて，なかなか腋に落ちたとはいいにくいものである．高校生や大学生にも納得してもらえるような説明を探していたが，ひとつの説明法を思いっいたので以下で紹介したい．
ポイントは 1 次式，2 次式の場合にしっかり納得することが大切で，一旦そこを納得してしまえば，一般の場合も容易に想像がつくようになるということである．
多項式 $f(x)$ が 1 次式，2 次式の場合に，具体的に $V(x)$ の動きを調べてみれば定理の背後にあるものを納得できるのである．

まず，1 次式の場合を考えよう．
1 次式の場合$f(x)=px+q$とする．
$p>0$ のとき

$p<0$ のときである．表からわかるように，零点通過のときに符号変化数 $V(x)$ の値が 1 減っている．
これらの表を左から右
$x$

$x$

$x$$($x$=$-$\infty$から $x=\infty)$ の方に見ていくと次のことが分かる．

$f(x)$ と $f’(x)$ は異符号が零点通過直前，$f’(x)$ の符号不変の零点通過の付近，$f(x)$が零点通過直後， と $f’(x)$ は同符号

これらの事から $V(x)$ は，区分的に整数値，非増加，右連続関数である．さらに，$V(-\infty)-V(\infty)=$ 方程式の解の全個数であることもわかる．これらの性質から $f(x)$ が 1 次式の場合に上記のスツルムの定理が成立している．また，これらの $V(x)$ が持つ性質は，重解をもたない一般の多項式が零点を通過する際に $f(x)$と $f’(x)$ の問の符号変化数に対しても成り立っことである．

60
次に，2 次式の場合を考えよう．
2 次式の場合

2 次の係数は 1 としても一般性を失わない．
$f(x)=x^{2}+px+q$$f’(x)=2x+p$ から，平方完成によりとおく．
$f(x)= \frac{f’(x)}{4}\cdot f’(x)-\frac{D}{4}$である．ただし，判別式 $D=p^{2}-4q$
したがって $\{f(x), f’(x), D\}$ がスツルム列となる．

$D>0$ のとき，相異なる 2 実根を$\beta$$\alpha,$とおく．

$f(x)$ の零点
表より，零点通過のときに符号変化数 $V(x)$ の値が 1 減っている．
$f’(x)$ の符号は零点の前後で変化しない．
$V(x)$ の$D>0$ であるので，の前後に注目する．
値の変動は $f(x)$ と $f’(x)$ の間の符号変化数の変動そのものである．ここの部分だけの表をみると，1 次関数のときの表と同じパターンであることがわかる．すなわち $V(x)$ は．1$(x=\alpha, \beta)$減るのである．
$f’(x)$ の零点 $(x=-p/2)$ の前後に注目する．この零点においては$f(x)= \frac{f’(x)}{4}\cdot f’(x)-\frac{d}{4}$と $D>0$ より，$f(x)$ は$0$ではなく $D$ と異符号の $f(x)<0$ である．
$x=-p/2$ の前後では上から$-\cdot$の形で，”. “ の中に“–“ ,‘’$0$“, “$+$”$+$が入っている．したがって，これらの符号変化数$V(x)$は，いずれも 1 である．結局，の前後で符号に変動があるが，$f(x)$ が零点を通過することに，$V(x)$ の値は 1 減るのでの値は不変である．したがって，$f’(x)$ の零点 $(X=-p/2)$

ある．

次に $D<0$ のとき，常に $f(x)>0$ であり，$f(x)=0$ は実数解をもたない．

61
上から，$+\cdot-$の形で，””の中に，“ $-,,$ ,“$0$“,$+$”が入っている．したがって，これらの符号変化数は，いずれも 1 である．結局，$f’(x)$ の零点 $(x=-p/2)$$V(x)$の前後で符号に変動があるが，の値は不変である．

つでに本論とははずれるが，参考のために，$D=0$ のときも調べておく．常に $f(x)\geq 0$$f(x)=0$ はただ一つの重解をもつ．
であり，となっている．

高校 1 年生で，2 次関数のさまざまな問題を解かせて多くの生徒を悩ましているが上記 3 つの表から，従来とはちがった観点でそれらの問題を系統的に解くことを可能とする．

本論にもどろう．実は，1 次関数と 2 次関数で考察したことで本質は尽きている．念のために，3 次関数の例を一つ挙げてこのことを確認しよう．

3 次式の場合$f(x)=x^{3}-x$とする．Sturm 列は$\{f_{0}(x), f_{1}(x), f_{2}(x), f_{3}(x)\}$である．ただし，$f_{0}(x)=f(x)=x^{3}-x$ ,12$f_{1}(x)=f’(x)=x--$3’$f_{2}(x)=x$ ,$f_{3}(x)=1$である．これより，

1 次式，2 次式のときと同様に，零点通過のときのみにに符号変化数 $V(x)$ の値が 1 減って$fi(x),$いるのがわかる．さらに，$f_{2}(x)$の零点の前後では符号変化数は変化しない．

625

スツルムの定理の証明
ここまでスツルムの定理について説明してきたが，定理の証明をする．スツルム列が 4つの関数$\{f_{0}(x), f_{1}(x), f_{2}(x), f_{3}(x)\}$からなる場合で説明する．ただし，
ノ 0(x)$:=f(x)$$f_{1}(x);=f’(x)$$f_{2}(x)$$;=$$q_{1}(x)f_{1}(x)-f_{2}(x)$$f_{2}(x)$$;=$$q_{2}(x)f_{2}(x)-f_{3}(x)$$f_{3}(x)$は$0$でない定数は多項式である．
$V(x)$ の値の変動がある可能性があるのは $f_{0}(x),$ $fi(x),$ $f_{2}(x)$ のどれかがになると$f_{0}(x)$ が零点を通過するときだけである$V(x)$ の値を変動させるのは，きである．しかし，ことを示していく．

である．なお$q_{1}(x),$$q_{2}(x)$$0$

まず，隣り合う 2 つの関数は同時に O にはならないことを確認しておく．なぜならば，もしも同時に O になったとすると，上式を芋づる式にたどっていくとが になってしまう．これは，が でない定数であることに矛盾する．したがって，関数値が 0 となる添え字は少なくとも 1 つおきに現れる．関数値が 0 となる添え字とその前後をまとめて 1 組と考えて符号変化数を考えればいい．
$f_{3}(x)$$[f_{0}(z)\neq 0$のとき$0$$f_{3}(x)$$0$$]$$fi(z),$ $f_{2}(z),$ $f_{3}(z)$ のどれかが$0$になったときのみ $x=z$ の前後で $V(x)$ の値が変動する可能性がある．

$f_{1}(z),$ $f_{2}(z),$ $f_{3}(z)$のどれかが$0$になったときは，その添え字は少なくとも一つおきにあらわれるので，0 となる添え字とその前後をまとめて 1 組と考えてブロックに小分けして符号変化数の変動を見れば良い．
このブロックの添え字の前後では値は でなく異符号である．関数の連続性より，$x=z$の前後でこのブロック内での符号変化数は変動しない．
ブロック以外のところは，$x=z$ で値は でないので，関数の連続性よりその前後で符号は変化せず，符号変化数に影響を及ぼさない．
$x=z$ の前後で $V(x)$ の値は変動しない．
したがって，この状況では，
$0$

$0$

$[f_{0}(z)=0$ のとき$f_{2}(z)$が$0$$]$になる可能性はあるが，上の議論をそのまま繰り返すことにより，$f_{i}(z),$ $f_{2}(z)$の部分は $x=z$ の前後で符号変化数は変動しない．
$f_{0}(z),$ $fi(z)$ の符号変化数の変動が $V(x)$ の変動となるのである．
$f_{i}(z)\neq 0$

する．

であることより，$f_{0}(z),$$fi(z)$ 部分の符号変化数は，$x=z$ のところで 1 減少

63
したがって，この状況では，$x=z$ を通過するとき $V(x)$ の値は 1 減少する．
$f(x)$ の零点を通過することに $V(x)$ の値は 1 減ることが示された．
以上より，
Sturm 列が 4 個の場合を説明したが，任意の個数であっても証明は全く同じである．

本内容の詳細について興味のある方は，数学セミナーの連載記事 ([5]) をご覧下さい．

参考文献[1] Eugene Salamin, Computation of30 (1976), no. 135, 565-570.

$\pi$

using arithmetic-geometric mean. Math. Comp.

[2] S. Kosugi, Y. Morita and S. Yotsutani, Stationary solutions to the one-dimensional
Cahn-Hilliard equation: proof by the complete elliptic integrals, Discrete Contin.
Dyn. Syst., 19 $(2ppppp007),$ $609-629$ .
$\pi-\pi$ の計算[3] 竹之内脩一伊藤隆，

アルキメデスから現代まで -, 共立出版，2007.

[4] 村井-松本 四ッ谷: 日本数学会 2010 年度年会$-$

函数方程式論分科会アブストラクト．

[5] 四ッ谷-村井: 楕円関数と友達になろう 一微分方程式の解の全体像を求めて–, 数学セミナー，2011 年 4 月号 2012 年 4 月号， 本評論社．
$\sim$$H$