四次方程式
四次方程式(よじほうていしき、quartic equation)とは、次数が 4 である代数方程式のことである。この項目では主に一変数の四次方程式を扱う。
の形で表現される。a4 で割り
x4 + A3 x3 + A2 x2 + A1 x + A 0 = 0 (A_n = \frac{a_n}{a_4})
の形にしても解は変わらないのでこの形で論じられることが多い。
一般的な四次方程式の解法は、ジェロラモ・カルダーノの弟子であるルドヴィコ・フェラーリによって発見され、カルダノの著書『アルス・マグナ』で概要が述べられた。カルダノは x, x2, x3 をそれぞれ、線分の長さ、一辺の長さが x の正方形の面積、一辺の長さが x の立方体の体積と対応させてとらえ、4次以上の方程式には意味がないと考えていたため、三次方程式と違って詳細には述べられていない。
しかし、カルダノの死後、ルネ・デカルトは著書『方法序説』の試論の一つである『幾何学』において定規とコンパスによる作図を論じ、長さ x の線分、長さ y の線分、長さ 1 の線分から長さ x y の線分が得られることを示している。これによると、長さ x の線分と長さ 1 の線分から長さ xn(n は任意の自然数)の線分の作図が可能であることが分かるため 4 次以上の方程式を解くことにも幾何学的な意味を与えることは可能であり、カルダノの捉え方は不十分であったことが分かる。
その後、四次方程式は三次方程式と同様に様々な解法が発見され、五次方程式の代数的解法の探索と合わせて詳細な研究が進められた。
の形の式は x2 を変数とする二次方程式と見ることができ、複二次方程式 (biquadratic equation)、左辺は複二次式と呼ばれる。二次方程式の解法を知っていれば簡単に解くことができる。
y  x2 と変換することで y に関する二次方程式
を得ることができ、この二次方程式を解くことによって解を求められる。
また、実数を係数とする複二次式
に対して、次のような二次式の積への因数分解もよく行われる。x2 の二次方程式とみたときの判別式
D > 0 であれば x2 について平方完成することにより
x^4 +A_2 x^2 + A_0 = \left(x^2 + {A_2 \over 2}\right)^2 - \frac \left(x^2 + \sqrt{A_0} \right)^2 - \left(2 \sqrt{A_0} - A_2 \right) x^2
と変形すれば、いずれの場合も因数分解の公式
を利用して実数を係数とする二次式の積に因数分解できる。
四次方程式は、代数学の基本定理より、高々4個の複素数解を持つ。
四次方程式 ax4 + bx3 + cx2 + dx + e  0 の判別式は
\begin{align}
\Delta=\; &256a^3e^3-192a^2bde^2-128a^2c^2e^2+144a^2cd^2e-27a^2d^4\\
&\ +144ab^2ce^2-6ab^2d^2e-80abc^2de+18abcd^3+16ac^4e\\
&\ -4ac^3d^2-27b^4e^2+18b^3cde-4b^3d^3-4b^2c^3e+b^2c^2d^2
\end{align}
によって与えられ、係数によって定まる以下の4個の定数によってさらに詳細な情報が得られる。
\begin{align}
P&=8ac-3b^2\\
R&=b^3+8a^2d-4abc\\
\Delta_0&=c^2-3bd+12ae\\
D&=64a^3e-16a^2c^2+16ab^2c-16a^2bd-3b^4
\end{align}
Δ, P, R, Δ0, D に関して、以下の事実が成立する。
Δ < 0 のとき、異なる2個の実数解と1組の共役複素数解を持つ。
Δ > 0 のとき、
P < 0 かつ D < 0 ならば、相異なる4個の実数解を持つ。
P > 0 または D > 0 ならば、2組の共役複素数解を持つ。
Δ = 0 のときにのみ、方程式は重解を持ち、
P < 0 かつ D < 0 かつ Δ0 ≠ 0 ならば、1個の実数二重解と、異なる2個の重複度 1 の実数解を持つ。
D > 0 または(P < 0 かつ(D, R のどちらかが0でない))ならば、1個の実数二重解と、1組の共役複素数解を持つ。
Δ0  0 かつ D ≠ 0ならば、1個の実数三重解と、1個の重複度 1 の実数解を持つ。
D  0 のとき、
P < 0 ならば、異なる 2個の実数二重解を持つ。
P < 0 かつ R  0 ならば、1組の共役複素数である、異なる 2個の虚数二重解を持つ。
Δ0  0 ならば、−b を実数四重解として持つ。
以上には、例えば Δ > 0 かつ P·D < 0 である場合などが記されていない。しかし、このような組み合わせは実際には存在しない。
フェラーリの解法は、一般的な四次方程式の解法のうちで最初に与えられた解法である。四次方程式
を a4 で割り
の形にする。(A_n = \frac{a_n}{a_4})
B_3 := \frac{A_3}{4}
となり、3次の項が消えた方程式が得られる。見やすいように
と書く。q  0 の時は、複二次式として解けばよいので、以後は q ≠ 0 とする。
媒介変数 u ≠ 0 を用い
\left(y^2 + \frac{p + u}{2}\right)^2 - u \left(y - \frac{q}{2u}\right)^2 = 0
と変形する。ここで上式を展開し係数を比較すると、u の三次方程式
が得られる。このような補助的な方程式を、与えられた四次方程式に関する三次分解方程式(resolvent cubic equation) という。q ≠ 0 なので、この分解方程式の解は u ≠ 0 を満たしており、この解の一つを u として取る。また、求める四次方程式は
\left\{ \left(y^2 + \frac{p + u}{2}\right) + \sqrt{u} \left(y - \frac{q}{2u}\right) \right\} \left\{ \left(y^2 + \frac{p + u}{2}\right) - \sqrt{u} \left(y - \frac{q}{2u}\right) \right\} = 0
となり、この2つの二次方程式から、四次方程式の解を求めることができる。
デカルトは、著書『方法序説』の試論の一つである『幾何学』において四次方程式
を解くために、二次式による因数分解
を仮定した方法を奨めた。係数を比較すると
が得られる。上の 3 つの式から
が得られる。
という、c1 に関する六次方程式が得られる。偶数次の項しか無いので u = c12 とでもおけば
という u に関する三次方程式が得られる。この方程式は、フェラーリの方法で得たのと同じ三次分解方程式であり、これを解くことによって、元の方程式の解が得られる。
レオンハルト・オイラーは、三次方程式や四次方程式の解法をいくつか発見した。ここに述べる方法もオイラーの方法と呼ばれる解法の一つである。
という等式を用いると、x を未知数とする四次方程式
であることが分かる。
この方程式と、3 次の項の消えた四次方程式
の係数を比べ、p, q, r から a, b, c を求めることができれば、3 次の項の消えた四次方程式の解は上にあるように 4 つ求まる。
ここで f0 = (2a)2, f1 = (2b)2, f2 = (2c)2 とおけば
となり、根と係数の関係により f0, f1, f2 は三次方程式
の解であり、これもフェラーリの方法に現れた三次分解方程式である。この三次方程式を解くことによって a, b, c が得られる。
ジョゼフ=ルイ・ラグランジュは、既に知られていた三次方程式や四次方程式の解法を、いろいろな視点から詳しく調べ上げた。ここで述べるのは、ラグランジュによるフェラーリの方法の解釈であり、現代的に言えば対称群を用いた方法である。
フェラーリの方法において、四次方程式は
y4 + p y2 + q y + r = 0
の形に変形される。この方程式の 4 つの解を r0, r1, r2, r3 とする。三次分解式を解くことで四次方程式は、 2 つの二次方程式
y^2 + \frac{p+u}{2} = \pm \sqrt{u} \left(y - \frac{q}{2u} \right)
に分解することができた。
y^2 + \frac{p+u}{2} = \sqrt{u} \left(y - \frac{q}{2u} \right)
は、元の四次方程式の 4 つの解のうちの 2 つを解とするが、これをとりあえず r0, r1 の 2 つとしたとき、
y^2 + \frac{p+u}{2} = - \sqrt{u} \left(y - \frac{q}{2u} \right)
の解は r2, r3 となり、根と係数の関係から
r_0 + r_1 = \sqrt{u}
r_2 + r_3 = - \sqrt{u}
(r0 + r1) (r2 + r3) = − u
y^2 + \frac{p+u}{2} = \sqrt{u} \left(y - \frac{q}{2u} \right)
の解を r0, r1 としたが、解の並び方はいろいろ考えられる。 rm と rn を入れ替える互換を σm,n と書けば、例えば
σ0,1 (r0 + r1) (r2 + r3) = (r0 + r1) (r2 + r3)
σ0,2 (r0 + r1) (r2 + r3) = (r2 + r1) (r0 + r3)
など、一般には異なる値を取ることになる。このように調べていくと 4 つの解の並び方は 4! = 24 通りあるが
(r0 + r1) (r2 + r3) = − u
の値は、最初の解の並べ方によって
s0 = (r0 + r1) (r2 + r3)
s1 = (r0 + r2) (r1 + r3)
s2 = (r0 + r3) (r1 + r2)
の 3 通りとなる。
例えば、互換 σ0,1 を作用させると、
σ0,1 s0 = s0
σ&0,1 s1 = s2
σ0,1 s2 = s1
となる。
一般に、互換 σm,n は s0, s1, s2 の並べ替えしかしないため s0, s1, s2 に関する基本対称式
s0 + s1 + s2
s0 s1 + s1 s2 + s2 s0
s0 s1 s2
は、互換 σm,n によって不変であり、 r0, r1, r2, r3 の基本対称式で書けることになる。
すなわち s0, s1, s2 の基本対称式は、最初に考えた四次方程式の係数 p, q, r で書ける。
u = − (r0 + r1) (r2 + r3)
は、根の並べ方によって 3 つの値 − s0, − s1, − s2 をとり、これらを解とする方程式
(u + s0) (u + s1) (u + s2) = 0
の左辺は u についての多項式として展開すると、その係数が p, q, r の多項式として書ける式である。この u に関する三次方程式こそ、フェラーリの方法で三次分解方程式として求められた方程式に他ならない。
このようにしてラグランジュは、四次方程式を解くための補助方程式である三次分解方程式の解が、元の四次方程式の解の多項式で書けることを発見し、補助方程式の次数が三次である理由を、根の置換という立場からはっきりと示した。
t_0 = \left(r_0 + r_1\right) - \left( r_2 +r_3\right)
t_1 = \left(r_0 + r_2\right) - \left( r_1 +r_3\right)
t_2 = \left(r_0 + r_3\right) - \left( r_1 +r_2\right)
とすれば、{t_0}^2, {t_1}^2, {t_2}^2 を解とする三次方程式で四次方程式を解くこともできる。ラグランジュは補助方程式の解を用いて、問題の方程式の解の公式を表現するのとは逆に、補助方程式の解を、元の方程式の解の整式(あるいは一般に有理式)として書けることが代数的に解ける理由と考え、特に
u = r_0 +i r_1 - r_2 - i r_3
の形の式、さらに一般に、n次方程式であれば 1の原始n乗根\zeta_n を用いて
u = \sum_{k=0}^{n-1} {\zeta_n}^k r_k = r_0 +\zeta_n r_1 + {\zeta_n}^2 r_2 + \cdots + {\zeta_n}^{n-2} r_{n-2} + {\zeta_n}^{n-1} r_{n-1}
の形の式の性質を詳しく調べたが、五次以上の代数方程式の代数的解法の発見には至らなかった。この形の式をラグランジュの分解式 (Lagrange resolvent) という。
五次以上の代数方程式の代数的解法の存在については、パオロ・ルフィニ、オーギュスタン=ルイ・コーシー、ニールス・アーベルらの研究がアーベル-ルフィニの定理として結実し、否定されることになるが、彼らの研究は、このようなラグランジュの研究を源流としている。
a_4x^4 + a_3x^3 + a_2x^2 + a_1x + a_0 = 0 \quad (a_4\ne 0)
の解の公式は以下の通りである:
\begin{align}
x_1=&\color{Red}{-\frac{a_3}{4 a_4 \\
&-\frac{1}{2} \sqrt{\color{Blue}{\frac{a_3^2}{4 a_4^2}+\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}+\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{2 a_2}{3 a_4}} \\
&-\frac{1}{2} \sqrt{\color{Green}{\frac{a_3^2}{2 a_4^2}-\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}-\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{-\frac{a_3^3}{a_4^3}+\frac{4 a_2 a_3}{a_4^2}-\frac{8 a_1}{a_4}{4 \sqrt{\frac{a_3^2}{4 a_4^2}+\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}+\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{2 a_2}{3 a_4}}-\frac{4 a_2}{3 a_4}} \\
\\
\end{align}
\begin{align}
x_2=&\color{Red}{-\frac{a_3}{4 a_4} \\
&-\frac{1}{2} \sqrt{\color{Blue}{\frac{a_3^2}{4 a_4^2}+\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}+\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{2 a_2}{3 a_4}} \\
&+\frac{1}{2} \sqrt{\color{Green}{\frac{a_3^2}{2 a_4^2}-\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}-\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{-\frac{a_3^3}{a_4^3}+\frac{4 a_2 a_3}{a_4^2}-\frac{8 a_1}{a_4}{4 \sqrt{\frac{a_3^2}{4 a_4^2}+\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}+\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{2 a_2}{3 a_4}}-\frac{4 a_2}{3 a_4}} \\
\\
\end{align}
\begin{align}
x_3=&\color{Red}{-\frac{a_3}{4 a_4} \\
&+\frac{1}{2} \sqrt{\color{Blue}{\frac{a_3^2}{4 a_4^2}+\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}+\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{2 a_2}{3 a_4}} \\
&-\frac{1}{2} \sqrt{\color{Green}{\frac{a_3^2}{2 a_4^2}-\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}-\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}+\frac{-\frac{a_3^3}{a_4^3}+\frac{4 a_2 a_3}{a_4^2}-\frac{8 a_1}{a_4}{4 \sqrt{\frac{a_3^2}{4 a_4^2}+\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}+\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{2 a_2}{3 a_4}}-\frac{4 a_2}{3 a_4}} \\
\\
\end{align}
\begin{align}
x_4=&\color{Red}{-\frac{a_3}{4 a_4} \\
&+\frac{1}{2} \sqrt{\color{Blue}{\frac{a_3^2}{4 a_4^2}+\frac{\sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}{3 \sqrt[3]{2} a_4}+\frac{\sqrt[3]{2} \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right)}{3 a_4 \sqrt[3]{2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4+\sqrt{\left(2 a_2^3-9 a_1 a_3 a_2-72 a_0 a_4 a_2+27 a_0 a_3^2+27 a_1^2 a_4\right){}^2-4 \left(a_2^2-3 a_1 a_3+12 a_0 a_4\right){}^3}}-\frac{2 a_2}{3 a_4}} \\