よりよいプレゼンテーションのために（for Japanese, for English）

自分で執筆したいのですが如何せん時間がないため、いくつか興味深いリンクを張っておきます。（勝手なため問題とあらばご連絡下さい）尚、これは人に読んでもらうためではなく、自分のためのメモであることも付け加えておきます。

1. プレゼンテーション入門


2. 私家版　学会発表マニュアル


3. 研究発表マニュアル［口頭発表編］

Steepest Descent method ( simple version )

\documentclass{jsarticle}
\usepackage{ascmac}
\usepackage{amsmath, amssymb}
\usepackage[dvipdfm]{graphicx,color}
%
\parindent=0pt
%
\def\mathvc#1{\mbox{\boldmath $#1$}}
\begin{document}
\section{Steepest Descent method (SD method, 最急降下法)}

\subsection{一般的SD法}
$n$次元ベクトル$\mathvc x$についての関数$f(\mathvc x)$について最小化することを考える。
\begin{equation}
\mathrm{grad}f(\mathvc x_k)\ \left(=\ \nabla f(\mathvc x_k)\right)\ =\ \mathvc 0\ \mbox{（収斂条件）}
\end{equation}
また解は次のように更新される。
\begin{equation}
\mathvc x_{k+1}=\mathvc x_k-\alpha_k\nabla f(\mathvc x_k)\ \mbox{（直線探索）}
\end{equation}
ここで$\alpha_k$は$f(\mathvc x_{k+1})$が最小化されるように決定される。一般的に最急降下法はこれを繰り返すことで最良の$\mathvc x$を探査する。

　この方法の場合、収束が遅いという問題が生じる。この問題を克服するため2次近似を行うニュートン法の改良が現在の主流になっている。



\subsection{2次近似の下でのニュートン法に基づく改良}
　関数$f(\mathvc x)$について$\mathvc x_k$の下でTaylor展開を行う。
\begin{equation}
f(\mathvc x)=f(\mathvc x_k)+(\mathvc x-\mathvc x_k)^T\nabla f(\mathvc x_k)+\frac12(\mathvc x-\mathvc x_k)^T\mathvc F(\mathvc x_k)(\mathvc x-\mathvc x_k)
\end{equation}
ここで$\mathvc F(\mathvc x_k)$は$f(\mathvc x)$について$\mathvc x=\mathvc x_k$におけるHessian matrixである($\mathvc F(\mathvc x_k)=\nabla\left(\nabla f(\mathvc x)\right)^T |_{\mathvc x=\mathvc x_k}$)。$f(\mathvc x)$の第一変分は次のように求められる。
\begin{eqnarray}
\delta f&=&\delta x^T \nabla f(\mathvc x_k)+\frac12\delta\mathvc x^T\mathvc F(\mathvc x_k)(\mathvc x-\mathvc x_k)+\frac12\underbrace{(\mathvc x-\mathvc x_k)^T\mathvc F(\mathvc x_k)\delta\mathvc x}_{\left[\delta\mathvc x^T\mathvc F(\mathvc x_k)(\mathvc x-\mathvc x_k)\right]^\ast} \nonumber \\
&=&\delta\mathvc x^T\left[\nabla f(\mathvc x_k)+\mathvc F(\mathvc x_k)(\mathvc x-\mathvc x_k)\right] \qquad (\because\mbox{実空間での場合を考える})
\end{eqnarray}

したがって最適な$\mathvc x_{k+1}$は次の通り求められる。
\begin{equation}
\mathvc x_{k+1}=\mathvc x_k-\left[\mathvc F(\mathvc x_k)\right]^{-1}\nabla f(\mathvc x_k)
\end{equation}

したがって2次で十分近似可能な場合には、一回の操作で最適解を得ることが可能となる。

　$f(\mathvc x)$が十分滑らかでかつ$\mathvc F(\mathvc x)$が正定値となる場合には、$\mathvc x_0$を十分最適解$\mathvc x$に近く選ぶことにより、次式が成立する。
\begin{equation}
\left|\mathvc x_{k+1}-\mathvc x\right|\leq M\left|\mathvc x_k-\mathvc x\right|^2
\end{equation}
この不等式のように、右辺のべき乗が2で成立する場合には、2次収束するといい、誤差が小さくなると加速度的に収束することが言える。一般的な最急降下法の場合1次収束であり、最適解近傍では急速に収束することが分かる。

　この方法は非常に強力な反面、Hessian matrixの逆行列を計算する必要があり全体として計算速度は速くならない。

\end{document}

Partial Differential Equation #01

\documentclass{jsarticle}
\usepackage{ascmac}
\usepackage{amsmath}
\usepackage{enumerate}
%
\parindent=0pt
\begin{document}
\section{2階線形偏微分方程式}
\subsection{標準形}
$n\ (n\ge2)$個の独立変数を$x_i\ (i=1,2,\cdots,n)$とし、$C^2$級の未知関数を$u(x_i, i=1,2,\cdots,n)$とする。2階線形偏微分方程式の一般表示を下に示す。
\begin{equation}
\sum_{i=1}^n\sum_{j=1}^n A_{ij}\frac{\partial^2 u}{\partial x_i\partial x_j}+\sum_i B_i\frac{\partial u}{\partial x_i}+Cu=r(x_1,x_2,\cdots,x_n)
\end{equation}

\subsubsection{分類}
2階線形偏微分方程式は、楕円形、放物形、双曲形の3つのパターンに分類される。行列$A_{ij}$の固有値が総て同符号である場合を楕円形、少なくとも1個の固有値が0である場合を放物形、$m$個の正の固有値と$n-m$個の負の固有値がある場合を双曲形という。$n=2$の場合に具体的に述べると次のようになる。

$x=x_1, y=x_2$としたとき次のように表される。
\begin{eqnarray}
L[u]&\equiv&a\frac{\partial^2 u}{\partial x^2}+2b\frac{\partial^2 u}{\partial x\partial y}+c\frac{\partial^2 u}{\partial y^2}+e\frac{\partial u}{\partial x}+f\frac{\partial u}{\partial y}+gu\\
L[u]&\equiv&\varphi(x,y)
\end{eqnarray}

\begin{boxnote}
この場合の$A_{ij}$は次のように表される。
\begin{equation}
A_{ij}=\left(
\begin{array}{cc}
a&b\\
b&c
\end{array}\right)
\end{equation}

この場合の固有値は永年方程式$\lambda^2-\mathrm{tr}[A_{ij}]\lambda+\mathrm{det}[A_{ij}]=0$より次のように求めることができる。
$$\frac12\left[1\pm\frac{\sqrt{(a-c)^2+4b^2}}{a+c}\right]$$

固有値の同符号については定数項部分(この場合$\mathrm{det}[A_{ij}]$である)を確認すればいいことは周知である。
\end{boxnote}

このとき$a,b,c,e,f,g$は定数で$a^2+b^2+c^2\not=0$とする。上の方程式は次のように分類される。
\begin{enumerate}[(a)]
\item $b^2-ac<0$のとき楕円形
\item $b^2-ac=0$のとき放物形
\item $b^2-ac>0$のとき双曲形
\end{enumerate}

\subsubsection{一次変換}
ここでは方程式$L[u]=\varphi(x,y)$の標準形を次のように改める。
\begin{equation}
a\frac{\partial^2 u}{\partial x^2}+2b\frac{\partial^2 u}{\partial x\partial y}+c\frac{\partial^2 u}{\partial y^2}=q(x,y,u,u_x,u_y)
\end{equation}

次の一次変換を導入する。
\begin{equation}
\xi=\alpha x+y,\hspace{2zw}\eta=\beta x+y\hspace{2zw}(\alpha\not=\beta)
\end{equation}

\begin{eqnarray}
\frac{\partial}{\partial x}&=&\frac{\partial\xi}{\partial x}\frac{\partial}{\partial\xi}+\frac{\partial\eta}{\partial x}\frac{\partial}{\partial\eta}=\alpha\frac{\partial}{\partial\xi}+\beta\frac{\partial}{\partial\eta}\\
\frac{\partial}{\partial y}&=&\frac{\partial}{\partial\xi}+\frac{\partial}{\partial\eta}
\end{eqnarray}
を導入すると各偏導関数は次のように改められる。
\begin{eqnarray}
\frac{\partial^2}{\partial x^2}&=&\left(\alpha\frac{\partial}{\partial\xi}+\beta\frac{\partial}{\partial\eta}\right)^2 \nonumber\\
&=&\alpha^2\frac{\partial^2}{\partial\xi^2}+2\alpha\beta\frac{\partial^2}{\partial\xi\partial\eta}+\beta^2\frac{\partial^2}{\partial\eta^2}\\
\frac{\partial^2}{\partial x\partial y}&=&\left(\alpha\frac{\partial}{\partial\xi}+\beta\frac{\partial}{\partial\eta}\right)\left(\frac{\partial}{\partial\xi}+\frac{\partial}{\partial\eta}\right) \nonumber\\
&=&\alpha\frac{\partial^2}{\partial\xi^2}+(\alpha+\beta)\frac{\partial^2}{\partial\xi\partial\eta}+\beta\frac{\partial^2}{\partial\eta^2}\\
\frac{\partial^2}{\partial y^2}&=&\left(\frac{\partial}{\partial\xi}+\frac{\partial}{\partial\eta}\right)^2 \nonumber\\
&=&\frac{\partial^2}{\partial\xi^2}+2\frac{\partial^2}{\partial\xi\partial\eta}+\frac{\partial^2}{\partial\eta^2}
\end{eqnarray}
したがって$v(\xi,\eta)=u(x,y)$としたとき一次変換の結果は次のように表される。
\begin{equation}
a_1\frac{\partial^2 v}{\partial\xi^2}+2b_1\frac{\partial^2 v}{\partial\xi\partial\eta}+c_1\frac{\partial^2 v}{\partial\eta^2}=q_1(\xi,\eta,v,v_\xi,v_\eta)
\end{equation}
ここで$a_1,b_1,c_1$は次に示すとおりである。
\begin{eqnarray*}
a_1&=&a\alpha^2+2b\alpha+c\\
b_1&=&a\alpha\beta+b(\alpha+\beta)+c\\
c_1&=&a\beta^2+2b\beta+c
\end{eqnarray*}


\subsubsection{楕円形の場合}
$\alpha,\beta$の値に二次方程式$at^2+2bt+c=0\ (a\not=0)$の解を当てる。すなわち、$b^2-ac<0$であることからも$\alpha,\beta$は互いに共役な複素数と分かる。同時に以下の内容に帰着する。
\begin{eqnarray}
a_1=c_1=0&,&\hspace{2zw} b_1=2(ac-b^2)/a\not=0\\
2b_1\frac{\partial^2 v}{\partial\xi\partial\eta}&=&q_1(\xi,\eta,v,v_\xi,v_\eta)
\end{eqnarray}

このとき次の一次変換を導入する。
\begin{equation}
\sigma=\xi+\eta,\hspace{2zw}\tau=i(\xi-\eta)
\end{equation}

\begin{eqnarray}
\frac{\partial}{\partial\xi}&=&\frac{\partial\sigma}{\partial\xi}\frac{\partial}{\partial\sigma}+\frac{\partial\tau}{\partial\xi}\frac{\partial}{\partial\tau}=\frac{\partial}{\partial\sigma}+i\frac{\partial}{\partial\tau}\\
\frac{\partial}{\partial\eta}&=&\frac{\partial}{\partial\sigma}-i\frac{\partial}{\partial\tau}\\
\frac{\partial^2}{\partial\xi\partial\eta}&=&\frac{\partial^2}{\partial\sigma^2}+\frac{\partial^2}{\partial\tau^2}
\end{eqnarray}

したがって、$w(\sigma,\tau)=v(\xi,\eta)$としたとき一次変換の結果は次のように表される。
\begin{equation}
\frac{\partial^2w}{\partial\sigma^2}+\frac{\partial^2w}{\partial\tau^2}=q_2(\sigma,\tau,w,w_\sigma,w_\tau)
\end{equation}
すなわち楕円形の一般系は次のように表すことができる。
\begin{equation}
\frac{\partial^2u}{\partial x^2}+\frac{\partial^2u}{\partial y^2}=p(x,y,u,u_x,u_y)
\end{equation}


\subsubsection{放物形の場合}
$\alpha,\beta$の値に二次方程式$at^2+2bt+c=0\ (a\not=0)$の解を当てる。すなわち、$b^2-ac=0$であることからも重解（$\alpha=\beta$）であるため、$\beta$は$\alpha$に独立の任意の複素数と改める。この場合次のことが分かる。

\begin{eqnarray}
a_1=0,c_1\not=0&,&\hspace{2zw} b_1=2(ac-b^2)/a=0\\
c_1\frac{\partial^2 v}{\partial\eta^2}&=&q_1(\xi,\eta,v,v_\xi,v_\eta)
\end{eqnarray}

すなわち放物形の一般系は次のように表すことができる。
\begin{equation}
\frac{\partial^2u}{\partial x^2}=p_1(x,y,u,u_x,u_y),\ \mbox{または}\ %
\frac{\partial^2u}{\partial y^2}=p_2(x,y,u,u_x,u_y)
\end{equation}


\subsubsection{双曲形の場合}
$\alpha,\beta$の値に二次方程式$at^2+2bt+c=0\ (a\not=0)$の解を当てる。すなわち、$b^2-ac>0$であることからも$\alpha,\beta$は異なる実数と分かる。同時に以下の内容に帰着する。
\begin{eqnarray}
a_1=c_1=0&,&\hspace{2zw} b_1=2(ac-b^2)/a\not=0\\
2b_1\frac{\partial^2 v}{\partial\xi\partial\eta}&=&q_1(\xi,\eta,v,v_\xi,v_\eta)
\end{eqnarray}
楕円形に類似することはわかりやすい。このとき次の一次変換を導入する。
\begin{equation}
\sigma=\xi+\eta,\hspace{2zw}\tau=\xi-\eta
\end{equation}

\begin{eqnarray}
\frac{\partial}{\partial\xi}&=&\frac{\partial\sigma}{\partial\xi}\frac{\partial}{\partial\sigma}+\frac{\partial\tau}{\partial\xi}\frac{\partial}{\partial\tau}=\frac{\partial}{\partial\sigma}+\frac{\partial}{\partial\tau}\\
\frac{\partial}{\partial\eta}&=&\frac{\partial}{\partial\sigma}-\frac{\partial}{\partial\tau}\\
\frac{\partial^2}{\partial\xi\partial\eta}&=&\frac{\partial^2}{\partial\sigma^2}-\frac{\partial^2}{\partial\tau^2}
\end{eqnarray}

したがって、$w(\sigma,\tau)=v(\xi,\eta)$としたとき一次変換の結果は次のように表される。
\begin{equation}
\frac{\partial^2w}{\partial\sigma^2}-\frac{\partial^2w}{\partial\tau^2}=q_2(\sigma,\tau,w,w_\sigma,w_\tau)
\end{equation}
すなわち双曲形の一般系は次のように表すことができる。
\begin{equation}
\frac{\partial^2u}{\partial x^2}-\frac{\partial^2u}{\partial y^2}=p(x,y,u,u_x,u_y)
\end{equation}

\subsubsection{演習問題}
次の偏微分方程式の標準形を求めよ。
\begin{enumerate}[(a)]
\item $u_{xx}+u_{xy}=u_x$

{\bf 解答例}　2次方程式$\lambda^2+\lambda=0$の解は0,-1と分かる（双曲形）。\\
次の一次変換を考える。
$$\xi=y,\hspace{2zw}\eta=-x+y$$
\begin{eqnarray*}
\frac{\partial}{\partial x}&=&\frac{\partial\xi}{\partial x}\frac{\partial}{\partial\xi}+\frac{\partial\eta}{\partial x}\frac{\partial}{\partial\eta}=-\frac{\partial}{\partial\eta}\\
\frac{\partial}{\partial y}&=&\frac{\partial}{\partial\xi}+\frac{\partial}{\partial\eta}\\
\frac{\partial^2}{\partial x^2}&=&\frac{\partial^2}{\partial\eta^2}\\
\frac{\partial^2}{\partial x\partial y}&=&-\frac{\partial^2}{\partial\xi\partial\eta}-\frac{\partial^2}{\partial\eta^2}
\end{eqnarray*}

であるため、$v(\xi,\eta)=u(x,y)$としたとき一次変換の結果は次のようになる。
$$v_{\xi\eta}=v_\eta$$

次の一次変換を考える。
$$\sigma=\xi+\eta,\hspace{2zw}\tau=\xi-\eta$$
\begin{eqnarray*}
\frac{\partial}{\partial\xi}&=&\frac{\partial\sigma}{\partial\xi}\frac{\partial}{\partial\sigma}+\frac{\partial\tau}{\partial\xi}\frac{\partial}{\partial\tau}=\frac{\partial}{\partial\sigma}+\frac{\partial}{\partial\tau}\\
\frac{\partial}{\partial\eta}&=&\frac{\partial}{\partial\sigma}-\frac{\partial}{\partial\tau}\\
\frac{\partial^2}{\partial\xi\partial\eta}&=&\frac{\partial^2}{\partial\sigma^2}-\frac{\partial^2}{\partial\tau^2}
\end{eqnarray*}
であるため、$w(\sigma,\tau)=v(\xi,\eta)$としたとき一次変換の結果は次のようになる。
$$w_{\sigma\sigma}-w_{\tau\tau}=w_\sigma-w_\tau$$
{\bf 解答終了}

\item $u_{xx}+2u_{xy}+u_{yy}=u_y$

{\bf 解答例}　二次方程式$\lambda^2+2\lambda+1=0$の解は$-1$(重解)と分かる（放物形）。\\
次の一次変換を考える。
$$\xi=-x+y,\hspace{2zw}\eta=y$$
\begin{eqnarray*}
\frac{\partial}{\partial x}&=&\frac{\partial\xi}{\partial x}\frac{\partial}{\partial\xi}+\frac{\partial\eta}{\partial x}\frac{\partial}{\partial\eta}=-\frac{\partial}{\partial\xi}\\
\frac{\partial}{\partial y}&=&\frac{\partial}{\partial\xi}+\frac{\partial}{\partial\eta}\\
\frac{\partial^2}{\partial x^2}&=&\frac{\partial^2}{\partial\xi^2}\\
\frac{\partial^2}{\partial x\partial y}&=&-\frac{\partial^2}{\partial\xi^2}-\frac{\partial^2}{\partial\xi\partial\eta}\\
\frac{\partial^2}{\partial y^2}&=&\frac{\partial^2}{\partial\xi^2}+2\frac{\partial^2}{\partial\xi\partial\eta}+\frac{\partial^2}{\partial\eta^2}
\end{eqnarray*}
であるため、$v(\xi,\eta)=u(x,y)$としたとき一次変換の結果は次のようになる。
$$v_{\eta\eta}=v_\xi+v_\eta$$
{\bf 解答終了}

\end{enumerate}

\newpage

\subsection{3次元直交曲線座標におけるLaplaceの演算子}
ここでは応用上重要な円柱座標系、極座標系に対するLaplaceの演算子について纏める。
\subsubsection{円柱座標系}
円柱座標系は次のように与えられる。
\begin{eqnarray}
x=\rho\cos \phi,\ y=\rho\sin \phi,\ z=z\\
(0\le\rho\le +\infty,\ 0\le\phi\le 2\pi,\ -\infty\le z\le +\infty) \nonumber
\end{eqnarray}
\begin{equation}
\rho=\sqrt{x^2+y^2},\hspace{2zw}\tan\phi=\frac{y}{x}
\end{equation}

\begin{eqnarray*}
\frac{\partial\rho}{\partial x}=\frac{x}{\rho}&,\hspace{2zw}&\frac{\partial\rho}{\partial y}=\frac{y}{\rho}\\
\frac{\partial^2\rho}{\partial x^2}&=&\frac{\partial}{\partial x}\left(\frac{x}{\rho}\right)=\frac1{\rho}-\frac{x}{\rho^2}\frac{x}{\rho}\\
&=&\frac{y^2}{\rho^3}\\
\frac{\partial^2\rho}{\partial y^2}&=&\frac{\partial}{\partial y}\left(\frac{y}{\rho}\right)=\frac1{\rho}-\frac{y}{\rho^2}\frac{y}{\rho}\\
&=&\frac{x^2}{\rho^3}
\end{eqnarray*}

\begin{eqnarray*}
\frac{\partial\phi}{\partial x}=\frac{\partial}{\partial x}\left(\tan^{-1}\frac{y}{x}\right)&=&-\frac{y}{x^2}\frac{1}{1+(y/x)^2}\\
&=&-\frac{y}{\rho^2}\\
\frac{\partial\phi}{\partial y}=\frac{\partial}{\partial y}\left(\tan^{-1}\frac{y}{x}\right)&=&\frac1{x}\frac{1}{1+(y/x)^2}\\
&=&\frac{x}{\rho^2}\\
\end{eqnarray*}
\begin{eqnarray*}
\frac{\partial^2\phi}{\partial x^2}=\frac{\partial}{\partial x}\left(-\frac{y}{\rho^2}\right)&=&2\frac{y}{\rho^3}\frac{x}{\rho}\\
&=&\frac{2xy}{\rho^4}\\
\frac{\partial^2\phi}{\partial y^2}=\frac{\partial}{\partial y}\left(\frac{x}{\rho^2}\right)&=&-2\frac{x}{\rho^3}\frac{y}{\rho}\\
&=&-\frac{2xy}{\rho^4}
\end{eqnarray*}

したがって、
\begin{equation}
\frac{\partial u}{\partial x}=\frac{\partial\rho}{\partial x}\frac{\partial u}{\partial\rho}+\frac{\partial\phi}{\partial x}\frac{\partial u}{\partial\phi}+\frac{\partial z}{\partial x}\frac{\partial u}{\partial z}=\frac{x}{\rho}\frac{\partial u}{\partial\rho}-\frac{y}{\rho^2}\frac{\partial u}{\partial\phi}
\end{equation}
\begin{eqnarray}
\frac{\partial^2u}{\partial x^2}&=&\frac{\partial}{\partial x}\left(\frac{\partial\rho}{\partial x}\frac{\partial}{\partial\rho}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}+\frac{\partial z}{\partial x}\frac{\partial}{\partial z}\right) u \nonumber \\
&=&\frac{\partial}{\partial x}\left(\frac{\partial\rho}{\partial x}\frac{\partial}{\partial\rho}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}\right) u \nonumber \\
&=&\frac{\partial^2\rho}{\partial x^2}\frac{\partial u}{\partial\rho}+\frac{\partial^2\phi}{\partial x^2}\frac{\partial u}{\partial\phi}+\frac{\partial\rho}{\partial x}\left(\frac{\partial\rho}{\partial x}\frac{\partial}{\partial\rho}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}\right)\frac{\partial u}{\partial\rho}+\frac{\partial\phi}{\partial x}\left(\frac{\partial\rho}{\partial x}\frac{\partial}{\partial\rho}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}\right)\frac{\partial u}{\partial\phi} \nonumber \\
&=&\left(\frac{\partial\rho}{\partial x}\right)^2\frac{\partial^2u}{\partial\rho^2}+2\frac{\partial\rho}{\partial x}\frac{\partial\phi}{\partial x}\frac{\partial^2 u}{\partial\rho\partial\phi}+\left(\frac{\partial\phi}{\partial x}\right)^2\frac{\partial^2u}{\partial\phi^2}+\frac{\partial^2\rho}{\partial x^2}\frac{\partial u}{\partial\rho}+\frac{\partial^2\phi}{\partial x^2}\frac{\partial u}{\partial\phi} \nonumber \\
&=&\frac{x^2}{\rho^2}\frac{\partial^2u}{\partial\rho}-\frac{2xy}{\rho^3}\frac{\partial^2 u}{\partial\rho\partial\phi}+\frac{y^2}{\rho^4}\frac{\partial^2u}{\partial\phi^2}+\frac{y^2}{\rho^3}\frac{\partial u}{\partial\rho}+\frac{2xy}{\rho^4}\frac{\partial u}{\partial\phi}
\end{eqnarray}

\begin{equation}
\frac{\partial u}{\partial y}=\frac{\partial\rho}{\partial y}\frac{\partial u}{\partial\rho}+\frac{\partial\phi}{\partial y}\frac{\partial u}{\partial\phi}+\frac{\partial z}{\partial y}\frac{\partial u}{\partial z}=\frac{y}{\rho}\frac{\partial u}{\partial\rho}+\frac{x}{\rho^2}\frac{\partial u}{\partial\rho}
\end{equation}

\begin{eqnarray}
\frac{\partial^2 u}{\partial y^2}&=&\frac{\partial}{\partial y}\left(\frac{\partial\rho}{\partial y}\frac{\partial}{\partial\rho}+\frac{\partial\phi}{\partial y}\frac{\partial}{\partial\phi}\right) u \nonumber \\
&=&\frac{\partial^2\rho}{\partial y^2}\frac{\partial u}{\partial\rho}+\frac{\partial^2\phi}{\partial y^2}\frac{\partial u}{\partial\phi}+\frac{\partial\rho}{\partial y}\left(\frac{\partial\rho}{\partial y}\frac{\partial}{\partial\rho}+\frac{\partial\phi}{\partial y}\frac{\partial}{\partial\phi}\right)\frac{\partial u}{\partial\rho}+\frac{\partial\phi}{\partial y}\left(\frac{\partial\rho}{\partial y}\frac{\partial}{\partial\rho}+\frac{\partial\phi}{\partial y}\frac{\partial}{\partial\phi}\right)\frac{\partial u}{\partial\phi} \nonumber \\
&=&\left(\frac{\partial\rho}{\partial y}\right)^2\frac{\partial^2u}{\partial\rho^2}+2\frac{\partial\rho}{\partial y}\frac{\partial\phi}{\partial y}\frac{\partial^2 u}{\partial\rho\partial\phi}+\left(\frac{\partial\phi}{\partial y}\right)^2\frac{\partial^2u}{\partial\phi^2}+\frac{\partial^2\rho}{\partial y^2}\frac{\partial u}{\partial\rho}+\frac{\partial^2\phi}{\partial y^2}\frac{\partial u}{\partial\phi} \nonumber \\
&=&\frac{y^2}{\rho^2}\frac{\partial^2 u}{\partial\rho^2}+\frac{2xy}{\rho^3}\frac{\partial^2 u}{\partial\rho\partial\phi}+\frac{x^2}{\rho^4}\frac{\partial^2 u}{\partial\phi^2}+\frac{x^2}{\rho^3}\frac{\partial u}{\partial\rho}-\frac{2xy}{\rho^4}\frac{\partial u}{\partial\phi}
\end{eqnarray}

したがってLaplaceの演算子は次のように表すことができる。

\begin{equation}
\Delta\equiv\frac{\partial^2 u}{\partial\rho^2}+\frac1{\rho}\frac{\partial u}{\partial\rho}+\frac1{\rho^2}\frac{\partial^2 u}{\partial\phi^2}+\frac{\partial^2 u}{\partial z^2}
\end{equation}

\newpage

\subsubsection{極座標系}
極座標系は次のように与えられる。
\begin{eqnarray}
x=r\cos \phi\sin \theta,\ y=r\sin \phi\sin \theta,\ z=r\cos \theta\\
(0\le r\le +\infty,\ 0\le\phi\le 2\pi,-\pi\le \theta\le +\pi) \nonumber
\end{eqnarray}
\begin{equation}
r=\sqrt{x^2+y^2+z^2},\hspace{2zw}\tan\phi=\frac{y}{x},\hspace{2zw}\cos \theta=\frac{z}{r}
\end{equation}

\begin{eqnarray*}
\frac{\partial r}{\partial x}=\frac{x}{r}&,\hspace{2zw}&\frac{\partial r}{\partial y}=\frac{y}{r},\hspace{2zw}\frac{\partial z}{\partial r}=\frac{z}{r}\\
\frac{\partial^2 r}{\partial x^2}&=&\frac{\partial}{\partial x}\left(\frac{x}{r}\right)=\frac1{r}-\frac{x}{r^2}\frac{x}{r}\\
&=&\frac{y^2+z^2}{r^3}\\
\frac{\partial^2 r}{\partial y^2}&=&\frac{\partial}{\partial y}\left(\frac{y}{r}\right)=\frac1{r}-\frac{y}{r^2}\frac{y}{r}\\
&=&\frac{x^2+z^2}{r^3}\\
\frac{\partial^2 r}{\partial z^2}&=&\frac{\partial}{\partial z}\left(\frac{z}{r}\right)=\frac1{r}-\frac{z}{r^2}\frac{z}{r}\\
&=&\frac{x^2+y^2}{r^3}
\end{eqnarray*}



\begin{eqnarray*}
\frac{\partial\phi}{\partial x}&=&\frac{\partial}{\partial x}\left(\tan^{-1}{\frac{y}{x}}\right)=-\frac{y}{x^2}\frac{1}{1+(y/x)^2}\\
&=&-\frac{y}{x^2+y^2}\\
\frac{\partial\phi}{\partial y}&=&\frac{\partial}{\partial y}\left(\tan^{-1}{\frac{y}{x}}\right)=\frac1{x}\frac{1}{1+(y/x)^2}\\
&=&\frac{x}{x^2+y^2}\\
\frac{\partial\phi}{\partial z}&=&0
\end{eqnarray*}


\begin{eqnarray*}
\frac{\partial\phi^2}{\partial x^2}&=&\frac{\partial}{\partial x}\left(-\frac{y}{x^2+y^2}\right)=\frac{2xy}{(x^2+y^2)^2}\\
\frac{\partial\phi^2}{\partial y^2}&=&\frac{\partial}{\partial y}\left(\frac{x}{x^2+y^2}\right)=-\frac{2xy}{(x^2+y^2)^2}
\end{eqnarray*}



\begin{eqnarray*}
\frac{\partial\theta}{\partial x}&=&\frac{\partial}{\partial x}\left(\cos^{-1}{\frac{z}{r}}\right)=-\frac{z}{r^2}\frac{x}{r}\cdot -\frac{1}{\sqrt{1-\left(\frac{z}{r}\right)^2}}\\
&=&\frac{xz}{r^2\sqrt{x^2+y^2}}\\
\frac{\partial\theta}{\partial y}&=&\frac{\partial}{\partial y}\left(\cos^{-1}{\frac{z}{r}}\right)=-\frac{z}{r^2}\frac{y}{r}\cdot -\frac{1}{\sqrt{1-\left(\frac{z}{r}\right)^2}}\\
&=&\frac{yz}{r^2\sqrt{x^2+y^2}}\\
\frac{\partial\theta}{\partial z}&=&\frac{\partial}{\partial z}\left(\cos^{-1}{\frac{z}{r}}\right)=\left(\frac1{r}-\frac{z}{r^2}\frac{z}{r}\right)\cdot -\frac1{\sqrt{1-\left(\frac{z}{r}\right)^2}}\\
&=&-\frac{\sqrt{x^2+y^2}}{r^2}
\end{eqnarray*}

\begin{eqnarray*}
\frac{\partial^2\theta}{\partial x^2}&=&\frac{\partial}{\partial x}\left(\frac{xz}{r^2\sqrt{x^2+y^2}}\right)=\frac{z}{r^2\sqrt{x^2+y^2}}-\frac{x}{r}\frac{2xz}{r^3\sqrt{x^2+y^2}}-\frac{x^2z}{r^2(x^2+y^2)^\frac32}\\
&=&\frac{(y^4-x^2y^2+y^2z^2-2x^4)z}{r^4(x^2+y^2)^\frac32}\\
\frac{\partial^2\theta}{\partial y^2}&=&\frac{\partial}{\partial y}\left(\frac{yz}{r^2\sqrt{x^2+y^2}}\right)=\frac{z}{r^2\sqrt{x^2+y^2}}-\frac{y}{r}\frac{2yz}{r^3\sqrt{x^2+y^2}}-\frac{y^2z}{r^2(x^2+y^2)^\frac32}\\
&=&\frac{(x^4-x^2y^2+x^2z^2-2y^4)z}{r^4(x^2+y^2)^\frac32}\\
\frac{\partial^2\theta}{\partial z^2}&=&\frac{\partial}{\partial z}\left(-\frac{\sqrt{x^2+y^2}}{r^2}\right)=\frac{2z\sqrt{x^2+y^2}}{r^4}
\end{eqnarray*}


したがって、

\begin{eqnarray}
\frac{\partial u}{\partial x}&=&\frac{\partial r}{\partial x}\frac{\partial u}{\partial r}+\frac{\partial\phi}{\partial x}\frac{\partial u}{\partial\phi}+\frac{\partial\theta}{\partial x}\frac{\partial u}{\partial\theta} \nonumber \\
&=&\frac{x}{r}\frac{\partial u}{\partial r}-\frac{y}{x^2+y^2}\frac{\partial u}{\partial\phi}+\frac{xz}{r^2\sqrt{x^2+y^2}}\frac{\partial u}{\partial\theta}\\
\frac{\partial u}{\partial x}&=&\frac{\partial r}{\partial y}\frac{\partial u}{\partial r}+\frac{\partial\phi}{\partial y}\frac{\partial u}{\partial\phi}+\frac{\partial\theta}{\partial y}\frac{\partial u}{\partial\theta} \nonumber \\
&=&\frac{y}{r}\frac{\partial u}{\partial r}+\frac{x}{x^2+y^2}\frac{\partial u}{\partial\phi}+\frac{yz}{r^2\sqrt{x^2+y^2}}\frac{\partial u}{\partial\theta}\\
\frac{\partial u}{\partial x}&=&\frac{\partial r}{\partial z}\frac{\partial u}{\partial r}+\frac{\partial\phi}{\partial z}\frac{\partial u}{\partial\phi}+\frac{\partial\theta}{\partial z}\frac{\partial u}{\partial\theta} \nonumber \\
&=&\frac{z}{r}\frac{\partial u}{\partial z}-\frac{\sqrt{x^2+y^2}}{r^2}\frac{\partial u}{\partial\theta}
\end{eqnarray}

\begin{eqnarray}
\frac{\partial^2 u}{\partial x^2}&=&\frac{\partial}{\partial x}\left(\frac{\partial r}{\partial x}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial x}\frac{\partial}{\partial\theta}\right)u \nonumber \\
&=&\frac{\partial^2 r}{\partial x^2}\frac{\partial u}{\partial r}+\frac{\partial^2\phi}{\partial x^2}\frac{\partial u}{\partial\phi}+\frac{\partial^2\theta}{\partial x^2}\frac{\partial u}{\partial\theta}+\frac{\partial r}{\partial x}\left(\frac{\partial r}{\partial x}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial x}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial r} \nonumber \\
&&{}+\frac{\partial\phi}{\partial x}\left(\frac{\partial r}{\partial x}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial x}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial\phi}+\frac{\partial\theta}{\partial x}\left(\frac{\partial r}{\partial x}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial x}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial x}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial\theta} \nonumber \\
&=&\left(\frac{\partial r}{\partial x}\right)^2\frac{\partial^2 u}{\partial r^2}+\left(\frac{\partial\phi}{\partial x}\right)^2\frac{\partial^2 u}{\partial\phi^2}+\left(\frac{\partial\theta}{\partial x}\right)^2\frac{\partial^2 u}{\partial\theta^2} \nonumber \\
&&{}+2\frac{\partial r}{\partial x}\frac{\partial\phi}{\partial x}\frac{\partial^2 u}{\partial r\partial\phi}+2\frac{\partial\phi}{\partial x}\frac{\partial\theta}{\partial x}\frac{\partial^2 u}{\partial\phi\partial\theta}\frac{\partial^2 u}{\partial\phi\partial\theta}+2\frac{\partial\theta}{\partial x}\frac{\partial r}{\partial x}\frac{\partial^2 u}{\partial\theta\partial r} \nonumber \\
&&{}+\frac{\partial^2 r}{\partial x^2}\frac{\partial u}{\partial r}+\frac{\partial^2\phi}{\partial x^2}\frac{\partial u}{\partial\phi}+\frac{\partial^2\theta}{\partial x^2}\frac{\partial u}{\partial\theta}\\
&=&\frac{x^2}{r^2}\frac{\partial^2 u}{\partial r^2}+\frac{y^2}{(x^2+y^2)^2}\frac{\partial^2 u}{\partial\phi^2}+\frac{x^2z^2}{r^4(x^2+y^2)}\frac{\partial^2 u}{\partial\theta^2} \nonumber \\
&&{}-\frac{2xy}{r(x^2+y^2)}\frac{\partial^2 u}{\partial r\partial\phi}-\frac{2xyz}{r^2(x^2+y^2)^\frac32}\frac{\partial^2 u}{\partial\phi\partial\theta}+\frac{2x^2z}{r^3\sqrt{x^2+y^2}}\frac{\partial^2 u}{\partial\theta\partial r} \nonumber \\
&&{}+\frac{y^2+z^2}{r^3}\frac{\partial u}{\partial r}+\frac{2xy}{(x^2+y^2)^2}\frac{\partial u}{\partial\phi}+\frac{(y^4-x^2y^2+y^2z^2-2x^4)z}{r^4(x^2+y^2)^\frac32}\frac{\partial u}{\partial\theta}
\end{eqnarray}


\begin{eqnarray}
\frac{\partial^2 u}{\partial y^2}&=&\frac{\partial}{\partial y}\left(\frac{\partial r}{\partial y}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial y}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial y}\frac{\partial}{\partial\theta}\right)u \nonumber \\
&=&\frac{\partial^2 r}{\partial y^2}\frac{\partial u}{\partial r}+\frac{\partial^2\phi}{\partial y^2}\frac{\partial u}{\partial\phi}+\frac{\partial^2\theta}{\partial y^2}\frac{\partial u}{\partial\theta}+\frac{\partial r}{\partial y}\left(\frac{\partial r}{\partial y}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial y}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial y}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial r} \nonumber \\
&&{}+\frac{\partial\phi}{\partial y}\left(\frac{\partial r}{\partial y}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial y}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial y}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial\phi}+\frac{\partial\theta}{\partial y}\left(\frac{\partial r}{\partial y}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial y}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial y}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial\theta} \nonumber \\
&=&\left(\frac{\partial r}{\partial y}\right)^2\frac{\partial^2 u}{\partial r^2}+\left(\frac{\partial\phi}{\partial y}\right)^2\frac{\partial^2 u}{\partial\phi^2}+\left(\frac{\partial\theta}{\partial y}\right)^2\frac{\partial^2 u}{\partial\theta^2} \nonumber \\
&&{}+2\frac{\partial r}{\partial y}\frac{\partial\phi}{\partial y}\frac{\partial^2 u}{\partial r\partial\phi}+2\frac{\partial\phi}{\partial y}\frac{\partial\theta}{\partial y}\frac{\partial^2 u}{\partial\phi\partial\theta}\frac{\partial^2 u}{\partial\phi\partial\theta}+2\frac{\partial\theta}{\partial y}\frac{\partial r}{\partial y}\frac{\partial^2 u}{\partial\theta\partial r} \nonumber \\
&&{}+\frac{\partial^2 r}{\partial y^2}\frac{\partial u}{\partial r}+\frac{\partial^2\phi}{\partial y^2}\frac{\partial u}{\partial\phi}+\frac{\partial^2\theta}{\partial y^2}\frac{\partial u}{\partial\theta}\\
&=&\frac{y^2}{r^2}\frac{\partial^2 u}{\partial r^2}+\frac{x^2}{(x^2+y^2)^2}\frac{\partial^2 u}{\partial\phi^2}+\frac{y^2z^2}{r^4(x^2+y^2)}\frac{\partial^2 u}{\partial\theta^2} \nonumber \\
&&{}+\frac{2xy}{r(x^2+y^2)}\frac{\partial^2 u}{\partial r\partial\phi}+\frac{2xyz}{r^2(x^2+y^2)^\frac32}\frac{\partial^2 u}{\partial\phi\partial\theta}+\frac{2y^2z}{r^3\sqrt{x^2+y^2}}\frac{\partial^2 u}{\partial\theta\partial r} \nonumber \\
&&{}+\frac{z^2+x^2}{r^3}\frac{\partial u}{\partial r}-\frac{2xy}{(x^2+y^2)^2}\frac{\partial u}{\partial\phi}+\frac{(x^4-x^2y^2+x^2z^2-2y^4)z}{r^4(x^2+y^2)^\frac32}\frac{\partial u}{\partial\theta}
\end{eqnarray}

\begin{eqnarray}
\frac{\partial^2 u}{\partial z^2}&=&\frac{\partial}{\partial z}\left(\frac{\partial r}{\partial z}\frac{\partial}{\partial r}+\frac{\partial\phi}{\partial z}\frac{\partial}{\partial\phi}+\frac{\partial\theta}{\partial z}\frac{\partial}{\partial\theta}\right)u \nonumber \\
&=&\frac{\partial^2 r}{\partial z^2}\frac{\partial u}{\partial r}+\frac{\partial^2\theta}{\partial z^2}\frac{\partial u}{\partial\theta}+\frac{\partial r}{\partial z}\left(\frac{\partial r}{\partial z}\frac{\partial}{\partial r}+\frac{\partial\theta}{\partial z}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial r} \nonumber \\
&&{}+\frac{\partial\phi}{\partial z}\left(\frac{\partial r}{\partial z}\frac{\partial}{\partial r}+\frac{\partial\theta}{\partial z}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial\phi}+\frac{\partial\theta}{\partial z}\left(\frac{\partial r}{\partial z}\frac{\partial}{\partial r}+\frac{\partial\theta}{\partial z}\frac{\partial}{\partial\theta}\right)\frac{\partial u}{\partial\theta} \nonumber \\
&=&\left(\frac{\partial r}{\partial z}\right)^2\frac{\partial^2 u}{\partial r^2}+\left(\frac{\partial\theta}{\partial z}\right)^2\frac{\partial^2 u}{\partial\theta^2}+2\frac{\partial\theta}{\partial z}\frac{\partial r}{\partial z}\frac{\partial^2 u}{\partial\theta\partial r}+\frac{\partial^2 r}{\partial z^2}\frac{\partial u}{\partial r}+\frac{\partial^2\theta}{\partial z^2}\frac{\partial u}{\partial\theta} \\
&=&\frac{z^2}{r^2}\frac{\partial^2 u}{\partial^2 r}+\frac{x^2+y^2}{r^4}\frac{\partial^2 u}{\partial\theta^2}-\frac{2z\sqrt{x^2+y^2}}{r^3}\frac{\partial^2 u}{\partial r\partial\theta}+\frac{x^2+y^2}{r^3}\frac{\partial u}{\partial r}+\frac{2\sqrt{x^2+y^2}}{r^3}\frac{\partial u}{\partial\theta}
\end{eqnarray}



\newpage

\begin{eqnarray}
\Delta&=&\frac{\partial^2}{\partial r^2}+\frac1{x^2+y^2}\frac{\partial^2}{\partial\phi^2}+\frac1{r^2}\frac{\partial^2}{\partial\theta^2}+\frac2{r}\frac{\partial}{\partial r}+\frac{z^3}{r^4\sqrt{x^2+y^2}}\frac{\partial}{\partial\theta} \nonumber \\
&=&\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\frac{\partial^2}{\partial\theta^2}+\frac{\cos \theta}{r^2\sin\theta}\frac{\partial}{\partial\theta}+\frac1{r^2\sin^2 \theta}\frac{\partial^2}{\partial\phi^2} \nonumber \\
&=&\frac1{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial}{\partial r}\right)+\frac1{r^2}\left[\frac1{\sin \theta}\frac{\partial}{\partial\theta}\left(\sin \theta\frac{\partial}{\partial\theta}\right)+\frac1{\sin^2\theta}\frac{\partial^2}{\partial\phi^2}\right]
\end{eqnarray}


\subsection{変数分離法}

2階線形偏微分方程式の解を求める方法のひとつ変数分離法について、実際にLaplaceの方程式について円柱座標系や極座標系を導入した変数分離法を用いて求める。

\subsubsection{円柱座標系を用いた場合}

円柱座標系におけるLaplaceの方程式は以下の通りである。
\begin{equation}
\frac{\partial^2 u}{\partial\rho^2}+\frac1{\rho}\frac{\partial u}{\partial\rho}+\frac1{\rho^2}\frac{\partial^2 u}{\partial\phi^2}+\frac{\partial^2 u}{\partial z^2}=0
\end{equation}

$u=u(\rho, \phi, z)$が以下の通りに変数分離できると仮定する。
\[
u=R(\rho)\varPhi(\phi)Z(z)
\]
このときLaplaceの方程式は以下の通りに式変形することができる。
\begin{eqnarray*}
\frac{\partial^2}{\partial\rho^2}R(\rho)\varPhi(\phi)Z(z)+\frac1{\rho}\frac{\partial}{\partial\rho}R(\rho)\varPhi(\phi)Z(z)+\frac1{\rho^2}\frac{\partial^2}{\partial\phi^2}R(\rho)\varPhi(\phi)Z(z)+\frac{\partial^2}{\partial z^2}R(\rho)\varPhi(\phi)Z(z)&=&0\\
\varPhi(\phi)Z(z)\frac{\partial^2}{\partial\rho^2}R(\rho)+\frac1{\rho}\varPhi(\phi)Z(z)\frac{\partial}{\partial\rho}R(\rho)+\frac1{\rho^2}R(\rho)Z(z)\frac{\partial^2}{\partial\phi^2}\varPhi(\phi)+R(\rho)\varPhi(\phi)\frac{\partial^2}{\partial z^2}Z(z)&=&0\\
\frac1{R(\rho)}\frac{\partial^2}{\partial\rho^2}R(\rho)+\frac1{R(\rho)}\frac1{\rho}\frac{\partial}{\partial\rho}R(\rho)+\frac1{\varPhi(\phi)}\frac1{\rho^2}\frac{\partial^2}{\partial\phi^2}\varPhi(\phi)+\frac1{Z(z)}\frac{\partial^2}{\partial z^2}Z(z)&=&0\\
\frac{R''}{R}+\frac1{\rho}\frac{R'}{R}+\frac{1}{\rho^2}\frac{\varPhi''}{\varPhi}+\frac{Z''}{Z}=0
\end{eqnarray*}

\begin{equation*}
\rho^2\left[\frac{R''}{R}+\frac1{\rho}\frac{R'}{R}+\frac{Z''}{Z}\right]=-\frac{\varPhi''}{\varPhi}
\end{equation*}
上式において右辺は$\phi$のみの関数であり、左辺は$\rho$と$z$のみの関数である。両辺の値が$\rho,\phi,z$に従わずに一定である（この値を$m^2$）ことより、次のように二分できることが分かる。
\begin{eqnarray*}
\varPhi''+m^2\varPhi&=&0\\
\frac{R''}{R}+\frac1{\rho}\frac{R'}{R}-\frac{m^2}{\rho^2}&=&-\frac{Z''}{Z}
\end{eqnarray*}
後者の式もまた二分することができる（この値を$-a^2$とおく）。したがって各変数ごとの２階常微分方程式に分けて考えることができることがわかる。またこのとき導入した$m,a$は分離定数と呼ばれる。
\begin{eqnarray}
\varPhi''+m^2\varPhi&=&0\\
Z''&=&a^2 Z\\
\rho^2R''+\rho R'+(a^2\rho^2-m^2)R&=&0
\end{eqnarray}

$\varPhi$の基本解は$\cos{m\phi}$、$\sin{m\phi}$であり、$Z$の基本解は$\cosh{az}$、$\sinh{az}$である。また$R$については変数$a\rho$（ただし$a$は定数）についてBesselの微分方程式の解であることが分かる。よって基本解は$J_m(a\rho)$、$Y_m(a\rho)$である。


\vspace{2em}

円柱座標系を用いた波動の方程式は以下のように表される。
\begin{equation}
\frac{\partial^2u}{\partial t^2}=c^2\left(\frac{\partial^2u}{\partial\rho^2}+\frac1{\rho}\frac{\partial u}{\partial\rho}+\frac1{\rho^2}\frac{\partial^2u}{\partial\phi^2}+\frac{\partial^2u}{\partial z^2}\right)
\end{equation}

$u=u(\rho, \phi, z, t)$が以下の通りに変数分離できると仮定する。
\[
u=T(t)R(\rho)\varPhi(\phi)Z(z)
\]

このとき波動の方程式は以下のように変形することができる。
\begin{equation}
\frac{T''}{T}=c^2\left\{\frac{R''}{R}+\frac1{\rho}\frac{R'}{R}+\frac1{\rho^2}\frac{\varPhi''}{\varPhi}+\frac{Z''}{Z}\right\}
\end{equation}

分離定数を$-c^2k^2$とおくと方程式は次のように二分できる。
\begin{eqnarray*}
T''&=&-c^2k^2T\\
\frac{Z''}{Z}&=&-k^2-\frac{R''}{R}-\frac1{\rho}\frac{R'}{R}-\frac1{\rho^2}\frac{\varPhi''}{\varPhi}
\end{eqnarray*}

後者の分離定数を$-a^2$とおくと再び二分することができる。
\begin{eqnarray*}
Z''&=&-a^2Z\\
\rho^2\frac{R''}{R}+\rho\frac{R'}{R}+\rho^2(k^2-a^2)&=&-\frac{\varPhi''}{\varPhi}
\end{eqnarray*}

後者の分離定数を$m^2$とおくと再び二分することができる。
\begin{eqnarray*}
\varPhi''&=&-m^2\varPhi\\
\rho^2R''+\rho R'+\left\{\rho^2(k^2-a^2)-m^2\right\}R&=&0
\end{eqnarray*}

$T$の基本解は、$\cos{ckt}$、$\sin{ckt}$であり、$Z$の基本解は$\cos{az}$、$\sin{az}$であり、$\varPhi$の基本解は$\cos{m\phi}$、$\sin{m\phi}$である。また$R$については変数$\sqrt{k^2-a^2}\rho$（ただし$a,k$ともに定数）についてBesselの微分方程式の解であることが分かる。よって基本解は$J_m(\sqrt{k^2-a^2}\rho)$、$Y_m(\sqrt{k^2-a^2}\rho)$である。




\begin{boxnote}
\subsection*{Besselの微分方程式について}
特殊関数（$\Gamma$関数、Legendre関数など）に属するものである。$x$を実変数、$\nu$を実定数としたとき、
\vspace{-8pt}
\begin{equation}
x^2\frac{\mathrm d^2y}{\mathrm dx^2}+x\frac{\mathrm dy}{\mathrm dx}+(x^2-\nu^2)y=0
\end{equation}
を$\nu$次のBesselの微分方程式といい、Bessel関数はこの一般解全般をいう。

\subsection*{Bessel関数の種類}
\begin{enumerate}[(i)]
\item 第一種のBessel関数 $J_\nu(x)$\\
級数解表現で定義される。
\vspace{-8pt}
\begin{equation}
J_\nu(x)=\left(\frac{x}2\right)^\nu\sum_{l=0}^{+\infty}\frac{(-1)^l}{l!\Gamma(\nu+l+1)}\left(\frac{x}2\right)^{2l}
\end{equation}
$\nu\not=$整数のとき、$J_\nu(x)$と$J_{-\nu}(x)$とは独立であるため$y=C_1J_\nu(x)+C_2J_{-\nu}(x)\ \mbox{($\nu\not=$整数)}$はBesselの微分方程式の一般解である。

\item 第二種のBessel関数（別名、$\nu$次のNeumann関数） $Y_\nu(x)$\\
第一種のBessel関数を用いて定義される（特殊解の一種）。
\vspace{-8pt}
\begin{equation}
Y_\nu(x)=\frac1{\sin{\nu\pi}}(\cos{\nu\pi}J_\nu(x)-J_{-\nu}(x))\ \mbox{($\nu\not=$整数)}
\end{equation}
$\nu\to n$（整数）の際にも$Y_n(x)$として定義される。

\item 第三種のBessel関数 $H_\nu^{(1)}(x)$：第一種のHankel関数、$H_\nu^{(2)}(x)$：第二種のHankel関数
\begin{equation}
H_\nu^{(1)}(x)=J_\nu(x)+iY_\nu(x),\hspace{2zw}H_\nu^{(2)}(x)=J_\nu(x)-iY_\nu(x)
\end{equation}
\end{enumerate}

\subsection*{Bessel関数の特徴}
\begin{enumerate}[(a)]
\item Bessel展開というものがある。

$[0,1]$で定義された区分的に滑らかな関数$f(x)$は次のようにBessel展開できる
\vspace{-8pt}
\begin{equation}
f(x)=\sum_{k=0}^{+\infty}A_kJ_n(\lambda_kx)
\end{equation}
ここで展開係数$A_k$は次のように表される。
\vspace{-8pt}
\begin{equation}
A_k=\frac2{\left[J_{n+1}(\lambda_k)\right]^2}\int_0^1f(x)J_n(\lambda_kx)x\,\mathrm dx
\end{equation}

\item 球Bessel関数(spherical Bessel function)\\
次のように定義される。
\vspace{-8pt}
\begin{eqnarray}
j_n(x)&=&\sqrt\frac{\pi}{2x}J_{n+\frac12}(x)\\
&=&(-1)^nx^n\left(\frac1{x}\frac{\mathrm d}{\mathrm dx}\right)^n\frac{\sin x}{x}
\end{eqnarray}
解析的に表すことができるため、Besselの微分方程式の解として一般的によく多用される。

\end{enumerate}
\end{boxnote}



\subsubsection{極座標を用いた場合}
極座標系におけるLaplaceの演算子について次式の通りである。
\begin{equation}
\Delta u\equiv \frac1{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial u}{\partial r}\right)+\frac1{r^2}\left[\frac1{\sin \theta}\frac{\partial}{\partial\theta}\left(\sin \theta\frac{\partial u}{\partial\theta}\right)+\frac1{\sin^2\theta}\frac{\partial^2 u}{\partial\phi^2}\right]
\end{equation}

$u=u(r, \theta, \phi)$が以下の通りに変数分離できると仮定する。
\[
u=R(r)\varTheta(\theta)\varPhi(\phi)
\]

このとき次のように纏めることができる。
\begin{equation}
\frac{\Delta u}{u}=\frac{R''}{R}+\frac2{r}\frac{R'}{R}+\frac1{r^2}\left[\frac{\varTheta''}{\varTheta}+\frac{\cos \theta}{\sin \theta}\frac{\varTheta'}{\varTheta}+\frac1{\sin^2\theta}\frac{\varPhi''}{\varPhi}\right]
\end{equation}

$\mu=\cos \theta$とおき、$\varTheta(\theta)=M(\mu)$とした場合次のように変形することができる。
\begin{eqnarray*}
\varTheta'&=&\frac{\mathrm d\varTheta}{\mathrm d\theta}\\
&=&\frac{\mathrm d\varTheta}{\mathrm d\mu}\frac{\mathrm d\mu}{\mathrm d\theta}\ =\ \frac{\mathrm dM}{\mathrm d\mu}(-\sin \theta)\\
&=&-\sin \theta M'\\
\varTheta''&=&-\cos \theta M'-\sin \theta \frac{\mathrm dM'}{\mathrm d\theta}\\
&=&-\cos \theta M'-\sin \theta \frac{\mathrm dM'}{\mathrm d\mu}\frac{\mathrm d\mu}{\mathrm d\theta}\\
&=&\sin^2\theta M''-\cos \theta M'
\end{eqnarray*}

\begin{eqnarray}
\frac{\varTheta''}{\varTheta}+\frac{\cos \theta}{\sin \theta}\frac{\varTheta'}{\varTheta}+\frac1{\sin^2\theta}\frac{\varPhi''}{\varPhi}&=&\frac{\sin^2\theta M''-\cos \theta M'}{M}+\frac{\cos \theta}{\sin \theta}\frac{-\sin \theta M'}{M}+\frac1{\sin^2\theta}\frac{\varPhi''}{\varPhi} \nonumber \\
&=&\sin^2\theta\frac{M''}{M}-2\cos \theta\frac{M'}{M}+\frac1{\sin^2\theta}\frac{\varPhi''}{\varPhi} \nonumber \\
&=&(1-\mu^2)\frac{M''}{M}-2\mu\frac{M'}{M}+\frac1{1-\mu^2}\frac{\varPhi''}{\varPhi}
\end{eqnarray}


極座標系においてLaplaceの方程式は次式のとおりに表現される。
\begin{equation}
r^2\frac{R''}{R}+2r\frac{R'}{R}+(1-\mu^2)\frac{M''}{M}-2\mu\frac{M'}{M}+\frac1{1-\mu^2}\frac{\varPhi''}{\varPhi}=0
\end{equation}

上式は次のように変形することができる。
\begin{equation}
(1-\mu^2)\left[r^2\frac{R''}{R}+2r\frac{R'}{R}+(1-\mu^2)\frac{M''}{M}-2\mu\frac{M'}{M}\right]=-\frac{\varPhi''}{\varPhi}
\end{equation}

分離定数を$m^2$（$m$は整数）とおくと次のように二分される。
\begin{eqnarray*}
\varPhi''&=&-m^2\varPhi\\
-r^2\frac{R''}{R}-2r\frac{R'}{R}&=&(1-\mu^2)\frac{M''}{M}-2\mu\frac{M'}{M}-\frac{m^2}{1-\mu^2}
\end{eqnarray*}

後者の分離定数を$-n(n+1)$（$n$は0または自然数）とおくと次のように二分される。
\begin{eqnarray*}
r^2R''+2rR'+n(n+1)R&=&0\\
(1-\mu^2)M''-2\mu M'+\left[n(n+1)-\frac{m^2}{1-\mu^2}\right]M&=&0
\end{eqnarray*}

$\varPhi$の基本解は$\cos{m\phi}$と$\sin{m\phi}$であり、$R$の基本解は$R^n$と$R^{-(n+1)}$である。$M$についての２階常微分方程式はLegendreの陪微分方程式である。したがって、$M$の基本解は$P_n^m(\mu)$と$Q_n^m(\mu)$といえる。

ここで$(\theta,\phi)$についてのある偏微分方程式をしめす。
\begin{equation}
\frac1{\sin \theta}\frac{\partial}{\partial\theta}\left(\sin \theta\frac{\partial Y_n}{\partial\theta}\right)+\frac1{\sin^2\theta}\frac{\partial^2 Y_n}{\partial\phi^2}+n(n+1)Y_n=0
\end{equation}
上式の解$Y_n(\theta,\phi)$は$n$次の球面調和関数という。具体的には以下の$2n+1$になる。
\[
P_n(\cos \theta),\ \cos{m\phi}P_n^m(\cos \theta)\ (n=1,2,3,\ldots,n),\ \sin{m\phi}P_n^m(\cos \theta)\ (n=1,2,3,\ldots,n)
\]


\vspace{2em}

極座標を用いた熱伝導の方程式は以下の通りである。
\begin{equation}
\frac{\partial u}{\partial t}=k\left[\frac1{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial u}{\partial r}\right)+\frac1{r^2\sin \theta}\frac{\partial}{\partial\theta}\left(\sin \theta\frac{\partial u}{\partial\theta}\right)+\frac1{r^2\sin^2\theta}\frac{\partial^2 u}{\partial\phi^2}\right]
\end{equation}

$u=u(r, \mu, \phi, t)$が以下の通りに変数分離できるものと仮定する。
\[
u=T(t)R(r)M(\mu)\varPhi(\phi)
\]
このとき熱伝導の方程式は以下の通りに変形することができる。

\begin{equation}
\frac{T'}{T}=k\left[\frac1{r^2R}\frac{\mathrm d}{\mathrm dr}(r^2R')+\frac1{r^2M}\frac{\mathrm d}{\mathrm d\mu}\left\{(1-\mu^2)M'\right\}+\frac1{r^2(1-\mu^2)}\frac{\varPhi''}{\varPhi}\right]
\end{equation}

分離定数を$-k\kappa^2$とおくと次のように二分される。
\begin{eqnarray*}
T'&=&-k\kappa^2T\\
(1-\mu^2)\left[\frac1{R}\frac{\mathrm d}{\mathrm dr}(r^2R')+\frac1{M}\frac{\mathrm d}{\mathrm d\mu}\left\{(1-\mu^2)M'\right\}+r^2\kappa^2\right]&=&-\frac{\varPhi''}{\varPhi}
\end{eqnarray*}

後者の分離定数を$m^2$（$m$は整数）としたとき次のように二分される。
\begin{eqnarray*}
\varPhi''&=&-m^2\varPhi\\
-\frac1{R}\frac{\mathrm d}{\mathrm dr}(r^2R')-r^2\kappa^2&=&\frac1{M}\frac{\mathrm d}{\mathrm d\mu}\left\{(1-\mu^2)M'\right\}-\frac{m^2}{1-\mu^2}
\end{eqnarray*}

後者の分離定数を$-n(n+1)$（$n$は0または自然数）としたとき次のように二分される。
\begin{eqnarray*}
\frac{\mathrm d}{\mathrm dr}(r^2R')+\left\{r^2\kappa^2-n(n+1)\right\}R&=&0\\
\frac{\mathrm d}{\mathrm d\mu}\left\{(1-\mu^2)M'\right\}+\left\{n(n+1)-\frac{m^2}{1-\mu^2}\right\}M&=&0
\end{eqnarray*}
前者はBesselの微分方程式であるため$R$の基本解は$j_n(\kappa r)$、$n_n(\kappa r)$(球Neumman関数と呼ぶ)であり、後者はLegendreの陪微分方程式であるため$M$の基本解は$P_n^m(\mu)$、$Q_n^m(\mu)$である。また$T$の基本解は$e^{-k\kappa t}$であり、$\varPhi$の基本解は$\cos{m\phi}$、$\sin{m\phi}$である。




\begin{boxnote}
\subsection*{Legendreの陪微分方程式について}
Legendreの陪微分方程式は以下の通りである。
\vspace{-8pt}
\begin{equation}
(1-x^2)\frac{\mathrm d^2y}{\mathrm dx^2}-2x\frac{\mathrm dy}{\mathrm dx}+\left\{n(n+1)-\frac{m^2}{1-x^2}\right\}y=0
\end{equation}
ここで$m,n$は0または自然数として$m\leq n$を満たすものとする。

\subsection*{Legendre関数について}
最初にLegendre関数について纏める。Legendre関数は以下の微分方程式の解で定義される。
\vspace{-8pt}
\begin{equation}
(1-x^2)\frac{\mathrm d^2y}{\mathrm dx^2}-2x\frac{\mathrm dy}{\mathrm dx}+n(n+1)y=0
\end{equation}

この方程式の級数解は以下の通りである。

\begin{enumerate}[(a)]
\item 第一種のLegendre関数
\vspace{-8pt}
\begin{equation}
P_n(x)=\frac{(2n)!}{2^n(n!)^2}\left\{x^n-\frac{n(n-1)}{2(2n-1)}x^{n-2}+\frac{n(n-1)(n-2)(n-3)}{2\cdot 4(2n-1)(2n-3)}x^{n-4}+\ldots\right\}
\end{equation}
\item 第二種のLegendre関数
\vspace{-8pt}
\begin{equation}
Q_n(x)=\frac{2^n(n!)^2}{(2n+1)!}\left\{x^n-\frac{(n+1)(n+2)}{2(2n+3)}\frac1{x^{n+1}}+\frac{(n+1)(n+2)(n+3)(n+4)}{2\cdot 4(2n+3)(2n+5)}\frac1{x^{n+3}}+\ldots\right\}
\end{equation}
\end{enumerate}
Rodriguesの公式を用いると考えやすい。
\vspace{-8pt}
\begin{equation}
P_n(x)=\frac1{2^nn!}\frac{\mathrm d^n}{\mathrm dx^n}(x^2-1)^n
\end{equation}

\subsection*{Legendreの陪関数}
さてLegendreの陪関数であるが、Legendreの陪微分方程式の解として定義されるが、次の変数変換を行うことでLgendreの微分方程式に変形することができる。
\vspace{-8pt}
\[
y=(1-x^2)^\frac{m}2\frac{\mathrm d^mw}{\mathrm dx^m}
\]
Legendreの陪関数は以下の二種類がある。
\[
P_n^m(x)=(1-x^2)^\frac{m}2\frac{\mathrm d^m P_n(x)}{\mathrm dx^m}\hspace{2zw}Q_n^m(x)=(1-x^2)^\frac{m}2\frac{\mathrm d^m Q_n(x)}{\mathrm dx^m}
\]
それぞれ、第一種のLegendreの陪関数、第二種のLegendreの陪関数と呼ぶ。
\end{boxnote}

\newpage

\subsubsection{演習問題}
極座標系を用いた波動方程式は以下の通りである。
\begin{equation}
\frac{\partial^2 u}{\partial t^2}=c^2\left[\frac1{r^2}\frac{\partial}{\partial r}\left(r^2\frac{\partial u}{\partial r}\right)+\frac1{r^2\sin \theta}\frac{\partial}{\partial\theta}\left(\sin \theta\frac{\partial u}{\partial\theta}\right)+\frac1{r^2\sin^2\theta}\frac{\partial^2 u}{\partial\phi^2}\right]
\end{equation}
$u=u(r, \mu, \phi, t)$が以下の通りに変数分離できるものと仮定する。
\[
u=T(t)R(r)M(\mu)\varPhi(\phi)
\]
このとき波動方程式は以下の通りに変形することができる。
\begin{equation}
\frac{T''}{T}=c^2\left[\frac1{r^2R}\frac{\mathrm d}{\mathrm dr}(r^2R')+\frac1{r^2M}\frac{\mathrm d}{\mathrm d\mu}\left\{(1-\mu^2)M'\right\}+\frac1{r^2(1-\mu^2)}\frac{\varPhi''}{\varPhi}\right]
\end{equation}分離定数を$-c^2\kappa^2$としたとき次のように二分される。
\begin{eqnarray*}
T''&=&-c^2\kappa^2T\\
(1-\mu^2)\left[\frac1{R}\frac{\mathrm d}{\mathrm dr}(r^2R')+\frac1{M}\frac{\mathrm d}{\mathrm d\mu}\left\{(1-\mu^2)M'\right\}+r^2\kappa^2\right]&=&-\frac{\varPhi''}{\varPhi}
\end{eqnarray*}

あとは熱伝導の方程式を解いたときと同様である。

したがって、$R$の基本解は$j_n(\kappa r)$、$n_n(\kappa r)$であり、$M$の基本解は$P_n^m(\mu)$、$Q_n^m(\mu)$であり、$\varPhi$の基本解は$\cos{m\phi}$、$\sin{m\phi}$であり、$T$の基本解は$\cos{c\kappa t}$、$\sin{c\kappa t}$である。



\subsection{参考文献}
培風館、堀口剛・海老澤丕道・福井芳彦　共著、応用数学講義［初版］、2000

\end{document}

Virial theorem about hydrogen atom (Simple ver.)

\documentclass{jsarticle}
\parindent=0pt
\pagestyle{empty}
\begin{document}
\section{Gamma equation}
\begin{equation}
\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}\,dt
\end{equation}
\begin{equation}
\int_0^\infty t^{z-1}e^{-\alpha t}\,dt=\frac{\Gamma(z)}{\alpha^z}
\end{equation}

\section{spherical harmonics}

\begin{eqnarray}
Y_{0,0}&=&\sqrt{\frac1{4\pi}}\\
Y_{1,0}&=&\sqrt{\frac3{4\pi}}\cos \theta\\
Y_{1,\pm1}&=&\mp\sqrt{\frac3{8\pi}}\sin \theta e^{\pm i\varphi}\\
Y_{2,0}&=&\sqrt{\frac5{16\pi}}(3\cos^2 \theta -1)\\
Y_{2,\pm1}&=&\mp\sqrt{\frac{15}{8\pi}}\sin \theta\cos \theta e^{\pm i\varphi}\\
Y_{2,\pm2}&=&\mp\sqrt{\frac{15}{32\pi}}\sin^2 \theta e^{\pm i2\varphi}
\end{eqnarray}

\section{H\_atom}
\subsection{$n=1, l=0$\hspace{1zw}(1s)}
\begin{equation}
\psi_\mathrm{1s}=2e^{-r}Y_{0,0}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{1s}|\psi_\mathrm{1s}\right>&=&4\int e^{-2r}|Y_{0,0}|^2\,dV\\
&=&4\int r^2e^{-2r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&4\frac{\Gamma(3)}{2^3}\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{1s}\left|-\frac12\nabla^2\right|\psi_\mathrm{1s}\right>\\
&=&-2\int e^{-r}Y_{0,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)e^{-r}Y_{0,0}\,dV\\
&=&-2\int\left(r^2-2r\right)e^{-2r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&-2\left[\frac{\Gamma(3)}{2^3}-2\frac{\Gamma(2)}{2^2}\right]\\
&=&-2\left(\frac14-\frac12\right)\\
&=&\frac12
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{1s}\left|-\frac1{r}\right|\psi_\mathrm{1s}\right>\\
&=&-\frac1{\pi}\int e^{-r}\frac1{r}e^{-r}\,dV\\
&=&-4\int re^{-2r}\,dr\int\!\!\!\int |Y_{0,0}|^2\,d\theta d\varphi\\
&=&-4\frac{\Gamma(2)}{2^2}\\
&=&-1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac12-1\\
&=&-\frac12
\end{eqnarray*}

\newpage

\subsection{$n=2,l=0$\hspace{1zw}(2s)}
\begin{equation}
\psi_\mathrm{2s}=\frac1{\sqrt 8}(2-r)e^{-\frac{r}2}Y_{0,0}
\end{equation}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2s}|\psi_\mathrm{2s}\right>&=&\frac18\int (2-r)^2e^{-r}|Y_{0,0}|^2\,dV\\
&=&\frac18\int (r^4-4r^3+4r^2)e^{-r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&\frac18\left[\Gamma(5)-4\Gamma(4)+4\Gamma(3)\right]\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2s}\left|-\frac12\nabla^2\right|\psi_\mathrm{2s}\right>\\
&=&-\frac1{16}\int (2-r)e^{-\frac{r}2}Y_{0,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)(2-r)e^{-\frac{r}2}Y_{0,0}\,dV\\
&=&-\frac1{16}\int (2-r)\left(-\frac14r^3+\frac52r^2-4r\right)e^{-r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin\theta\,d\theta d\varphi\\
&=&-\frac1{16}\int \left(\frac14 r^4-3 r^3+9r^2-8r\right)e^{-r}\,dr\\
&=&-\frac1{16}\left[\frac14\Gamma(5)-3\Gamma(4)+9\Gamma(3)-8\Gamma(2)\right]\\
&=&-\frac{6-18+18-8}{16}\\
&=&\frac18
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2s}\left|-\frac1{r}\right|\psi_\mathrm{2s}\right>\\
&=&-\frac18\int \frac1{r}(2-r)^2e^{-r}|Y_{0,0}|^2\,dV\\
&=&-\frac18\int (r^3-4r^2+4r)e^{-r}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&-\frac18\left[\Gamma(4)-4\Gamma(3)+4\Gamma(2)\right]\\
&=&-\frac{6-8+4}{8}\\
&=&-\frac14
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}

\newpage

\subsection{$n=2,l=1$\hspace{1zw}(2p)}
\begin{equation}
\psi_\mathrm{2p}=\frac1{\sqrt{24}}re^{-\frac{r}2}Y_{1,m}
\end{equation}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2p}|\psi_\mathrm{2p}\right>&=&\frac1{24}\int r^2e^{-r}|Y_{1,m}|^2,dV\\
&=&\frac1{24}\int r^4e^{-r}\,dr\int\!\!\!\int |Y_{1,m}|^2\,d\theta d\varphi\\
&=&\frac1{24}\Gamma(5)\\
&=&\frac1{24}\times 24\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p}\left|-\frac12\nabla^2\right|\psi_\mathrm{2p}\right>\\
&=&-\frac1{48}\int re^{-\frac{r}2}Y_{1,m}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)re^{-\frac{r}2}Y_{1,m}\,dV\\
&=&-\frac1{48}\int \left(\frac14r^4-2r^3\right)e^{-r}\,dr \int\!\!\!\int|Y_{1,m}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{48}\left[\frac14\Gamma(5)-2\Gamma(4)\right]\\
&=&-\frac{6-12}{48}\\
&=&\frac1{8}
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p}\left|-\frac1{r}\right|\psi_\mathrm{2p}\right>\\
&=&-\frac1{24}\int re^{-r}|Y_{1,m}|^2\,dV\\
&=&-\frac1{24}\int r^3e^{-r}\,dr\int\!\!\!\int |Y_{1,m}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{24}\Gamma(4)\\
&=&-\frac14
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}

\newpage

\subsection{$n=3,l=0$\hspace{1zw}(3s)}
\begin{equation}
\psi_\mathrm{3s}=\frac2{\sqrt{6561\times 3}}\left(27-18r+2r^2\right)e^{-\frac{r}3}Y_{0,0}
\end{equation}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3s}|\psi_\mathrm{3s}\right>&=&\frac4{6561\times 3}\int \left(27-18r+2r^2\right)^2e^{-\frac{2r}3}|Y_{0,0}|^2\,dV\\
&=&\frac4{6561\times 3}\int \left(729-972r+432r^2-72r^3+4r^4\right)r^2 e^{-\frac{2r}3}\,dr\int\!\!\!\int |Y_{0,0}|^2\,d\theta d\varphi\\
&=&\frac4{6561\times 3}\left[729\Gamma(3)\left(\frac32\right)^3-972\Gamma(4)\left(\frac32\right)^4+432\Gamma(5)\left(\frac32\right)^5-72\Gamma(6)\left(\frac32\right)^6+4\Gamma(7)\left(\frac32\right)^7\right]\\
&=&\frac4{3^9}\left(\frac{3^9}4-\frac{3^{10}}2+3^9\times 4-3^9\times 5+\frac52\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3s}\left|-\frac12\nabla^2\right|\psi_\mathrm{3s}\right>\\
&=&-\frac2{6561\times 3}\int \left(27-18r+2r^2\right)e^{-\frac{r}3}Y_{0,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)\left(27-18r+2r^2\right)e^{-\frac{r}3}Y_{0,0}\,dV\\
&=&-\frac2{6561\times 3}\int \left(27-18r+2r^2\right)\left(\frac29r^2-6r+39-\frac{54}{r}\right)e^{-\frac23r} r^2\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&-\frac2{6561\times 3}\int \left(\frac49r^6-16r^5+192r^4-4\times3^5r^3+25\times 3^4r^2-2\times 3^6r\right)e^{-\frac23r}\,dr\\
&=&-\frac2{6561\times 3}\left[\frac49\Gamma(7)\left(\frac32\right)^7-16\Gamma(6)\left(\frac32\right)^6+192\Gamma(5)\left(\frac32\right)^5\right.\\
&&{}-\left.4\times3^5\Gamma(4)\left(\frac32\right)^4+25\times 3^4\Gamma(3)\left(\frac32\right)^3-2\times 3^6\Gamma(2)\left(\frac32\right)^2\right]\\
&=&-\frac2{6561\times 3}\left(\frac52\times 3^7-10\times 3^7+16\times 3^7-\frac12\times 3^{10}+\frac{25}4\times 3^7-\frac12\times 3^8\right)\\
&=&\frac2{6561\times 3}\times\frac14\times 3^7\\
&=&\frac1{18}
\end{eqnarray*}

\newpage
\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3s}\left|-\frac1{r}\right|\psi_\mathrm{3s}\right>\\
&=&-\frac4{6561\times 3}\int \frac1{r}\left(27-18r+2r^2\right)^2e^{-\frac23r}|Y_{0,0}|^2\,dV\\
&=&-\frac4{6561\times 3}\int \left(729r-972r^2+432r^3-72r^4+4r^5\right)e^{-\frac{2r}3}\,dr\int\!\!\!\int |Y_{0,0}|^2\sin \theta\,d\theta d\varphi\\
&=&\frac4{6561\times 3}\left[729\Gamma(2)\left(\frac32\right)^2-972\Gamma(3)\left(\frac32\right)^3+432\Gamma(4)\left(\frac32\right)^4-72\Gamma(5)\left(\frac32\right)^5+4\Gamma(6)\left(\frac32\right)^6\right]\\
&=&\frac4{3^9}\left(\frac{3^8}4-3^8+2\times 3^8-2\times 3^8+\frac52\times 3^7\right)\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{$n=3,l=1$\hspace{1zw}{3p}}
\begin{equation}
\psi_\mathrm{3p}=\frac4{\sqrt{6561\times 6}}r\left(6-r\right)e^{-\frac{r}3}Y_{1,m}
\end{equation}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3p}|\psi_\mathrm{3p}\right>&=&\frac{16}{6561\times 6}\int r^2\left(6-r\right)^2e^{-\frac23r}|Y_{1,m}|^2\,dV\\
&=&\frac{16}{6561\times 6}\int (r^6-12r^5+36r^4)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,m}|^2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 6}\left[\Gamma(7)\left(\frac32\right)^7-12\Gamma(6)\left(\frac32\right)^6+36\Gamma(5)\left(\frac32\right)^5\right]\\
&=&\frac{16}{6561\times 6}\left(\frac58\times 3^9-\frac52\times 3^8+3^8\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p}\left|-\frac12\nabla^2\right|\psi_\mathrm{3p}\right>\\
&=&-\frac8{6561\times 6}\int r\left(6-r\right)e^{-\frac{r}3}Y_{1,m}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r\left(6-r\right)e^{-\frac{r}3}Y_{1,m}\,dV\\
&=&-\frac8{6561\times 6}\int (6r-r^2)\left(-\frac19r^4+\frac83r^3-12r^2\right)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,m}|^2\sin \theta \,d\theta\,d\varphi\\
&=&-\frac8{6561\times 6}\int\left(\frac19r^6-\frac{10}3r^5+28r^4-72r^3\right)e^{-\frac23r}\,dr\\
&=&-\frac8{6561\times 6}\left[\frac19\Gamma(7)\left(\frac32\right)^7-\frac{10}3\Gamma(6)\left(\frac32\right)^6+28\Gamma(5)\left(\frac32\right)^5-72\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac8{6561\times 6}\left(\frac58\times 3^7-\frac{25}4\times 3^6+7\times 3^6-3^7\right)\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p}\left|-\frac1{r}\right|\psi_\mathrm{3p}\right>\\
&=&-\frac{16}{6561\times 6}\int r^3(6-r)^2e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,m}|^2\,d\theta\,d\varphi\\
&=&-\frac{16}{6561\times 6}\int (r^5-12r^4+36r^3)e^{-\frac23r}\,dr\\
&=&-\frac{16}{6561\times 6}\left[\Gamma(6)\left(\frac32\right)^6-12\Gamma(5)\left(\frac32\right)^5+36\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac{16}{6561\times 6}\left(\frac58\times 3^7-3^7+\frac12\times 3^7 \right)\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{$n=3,l=2$\hspace{1zw}(3d)}
\begin{equation}
\psi_\mathrm{3d}=\frac{4}{\sqrt{6561\times 30}}r^2e^{-\frac{r}3}Y_{2,m}
\end{equation}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3d}|\psi_\mathrm{3d}\right>&=&\frac{16}{6561\times 30}\int r^4e^{-\frac23r}|Y_{2,m}|^2\,dV\\
&=&\frac{16}{6561\times 30}\int r^6e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,m}|^2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 30}\Gamma(7)\left(\frac32\right)^7\\
&=&\frac{16}{6561\times 30}\left(\frac58\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d}\left|-\frac12\nabla^2\right|\psi_\mathrm{3d}\right>\\
&=&-\frac8{6561\times 30}\int r^2e^{-\frac{r}3}Y_{2,m}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r^2e^{-\frac{r}3}Y_{2,m}\,dV\\
&=&-\frac8{6561\times 30}\int \left(\frac19r^6-2r^5\right)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,m}|\,d\theta d\varphi\\
&=&-\frac8{6561\times 30}\left[\frac19\Gamma(7)\left(\frac32\right)^7-2\Gamma(6)\left(\frac32\right)^6\right]\\
&=&-\frac8{6561\times 30}\left(\frac58\times 3^7-\frac54\times 3^7\right)\\
&=&\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d}\left|-\frac1{r}\right|\psi_\mathrm{3d}\right>\\
&=&-\frac{16}{6561\times 30}\int r^3e^{-\frac23r}|Y_{2,m}|^2\,dV\\
&=&-\frac{16}{6561\times 30}\int r^5e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,m}|^2\sin \theta\,d\theta d\varphi\\
&=&-\frac{16}{6561\times 30}\Gamma(6)\left(\frac32\right)^6\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\end{document}

Virial theorem about Hydrogen atom

\documentclass{jsarticle}
\parindent=0pt
\pagestyle{empty}
\begin{document}
\section{Gamma equation}
\begin{equation}
\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}\,dt
\end{equation}
\begin{equation}
\int_0^\infty t^{z-1}e^{-\alpha t}\,dt=\frac{\Gamma(z)}{\alpha^z}
\end{equation}

\section{spherical harmonics}

\begin{eqnarray}
Y_{0,0}&=&\sqrt{\frac1{4\pi}}\\
Y_{1,0}&=&\sqrt{\frac3{4\pi}}\cos \theta\\
Y_{1,\pm1}&=&\mp\sqrt{\frac3{8\pi}}\sin \theta e^{\pm i\varphi}\\
Y_{2,0}&=&\sqrt{\frac5{16\pi}}(3\cos^2 \theta -1)\\
Y_{2,\pm1}&=&\mp\sqrt{\frac{15}{8\pi}}\sin \theta\cos \theta e^{\pm i\varphi}\\
Y_{2,\pm2}&=&\mp\sqrt{\frac{15}{32\pi}}\sin^2 \theta e^{\pm i2\varphi}
\end{eqnarray}

\section{H\_atom}
\subsection{1s}
\begin{equation}
\psi_\mathrm{1s}=\frac1{\sqrt{\pi}}e^{-r}
\end{equation}
\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{1s}|\psi_\mathrm{1s}\right>&=&\frac1{\pi}\int e^{-2r}\,dV\\
&=&4\int r^2e^{-2r}\,dr\\
&=&4\frac{\Gamma(3)}{2^3}\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{1s}\left|-\frac12\nabla^2\right|\psi_\mathrm{1s}\right>\\
&=&-\frac1{2\pi}\int e^{-r}\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)e^{-r}\,dV\\
&=&2\int e^{-r}\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}\right)e^{-r}r^2\,dr\\
&=&-2\int e^{-r}\left(1-\frac{2}{r}\right)e^{-r}r^2\,dr\\
&=&-2\int (r^2-2r)e^{-2r}\,dr\\
&=&-2\left[\frac{\Gamma(3)}{2^3}-2\frac{\Gamma(2)}{2^2}\right]\\
&=&-2\left(\frac14-\frac12\right)\\
&=&\frac12
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{1s}\left|-\frac1{r}\right|\psi_\mathrm{1s}\right>\\
&=&-\frac1{\pi}\int e^{-r}\frac1{r}e^{-r}\,dV\\
&=&-4\int \frac1{r}e^{-2r}r^2\,dr\\
&=&-4\int r e^{-2r}\,dr\\
&=&-4\frac{\Gamma(2)}{2^2}\\
&=&-1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac12-1\\
&=&-\frac12
\end{eqnarray*}

\newpage

\subsection{2s}
\begin{equation}
\psi_\mathrm{2s}=\frac1{\sqrt{32\pi}}(2-r)e^{-\frac{r}2}
\end{equation}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2s}|\psi_\mathrm{2s}\right>&=&\frac1{32\pi}\int (2-r)^2e^{-r}\,dV\\
&=&\frac18\int (r^4-4r^3+4r^2)e^{-r}\,dr\\
&=&\frac18\left[\Gamma(5)-4\Gamma(4)+4\Gamma(3)\right]\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2s}\left|-\frac12\nabla^2\right|\psi_\mathrm{2s}\right>\\
&=&-\frac1{64\pi}\int (2-r)e^{-\frac{r}2}\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)(2-r)e^{-\frac{r}2}\,dV\\
&=&-\frac1{16}\int r^2(2-r)e^{-\frac{r}2}\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}\right)(2-r)e^{-\frac{r}2}\,dr\\
&=&-\frac1{16}\int r^2(2-r)e^{-\frac{r}2}\left(1+\frac{2-r}4-\frac2{r}-\frac{2-r}{r}\right)e^{-\frac{r}2}\,dr\\
&=&-\frac1{16}\int \left(\frac14 r^4-3 r^3+9r^2-8r\right)e^{-r}\,dr\\
&=&-\frac1{16}\left[\frac14\Gamma(5)-3\Gamma(4)+9\Gamma(3)-8\Gamma(2)\right]\\
&=&-\frac{6-18+18-8}{16}\\
&=&\frac18
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2s}\left|-\frac1{r}\right|\psi_\mathrm{2s}\right>\\
&=&-\frac1{32\pi}\int (2-r)e^{-\frac{r}2}\frac1{r}(2-r)e^{-\frac{r}2}\,dV\\
&=&-\frac18\int \frac{(2-r)^2}{r}e^{-r}r^2\,dr\\
&=&-\frac18\int (r^3-4r^2+4r)e^{-r}\,dr\\
&=&-\frac18\left[\Gamma(4)-4\Gamma(3)+4\Gamma(2)\right]\\
&=&-\frac{6-8+4}{8}\\
&=&-\frac14
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}

\newpage

\subsection{2px}
\begin{eqnarray}
\psi_\mathrm{2p_x}&=&\frac1{\sqrt{32\pi}}re^{-\frac{r}2}\sin \theta\cos \varphi\\
&=&\frac1{\sqrt{24}}re^{-\frac{r}2}\frac1{\sqrt 2}\left(-Y_{1,+1}+Y_{1,-1}\right)
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2p_x}|\psi_\mathrm{2p_x}\right>&=&\frac1{32\pi}\int r^2e^{-r}\sin^2 \theta \cos^2 \varphi\,dV\\
&=&\frac1{32\pi}\int r^4e^{-r}\,dr\int_0^{\pi}\sin^3 \theta\,d\theta\int_0^{2\pi}\cos^2\varphi\,d\varphi\\
&=&\frac1{32\pi}\Gamma(5)\left[-\cos \theta+\frac13\cos^3 \theta\right]_0^\pi\left[\frac{1+\cos{2\varphi}}2\right]_0^{2\pi}\\
&=&\frac1{32\pi}\times 24\times \frac43 \times \pi\\
&=&1
\end{eqnarray*}


\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p_x}\left|-\frac12\nabla^2\right|\psi_\mathrm{2p_x}\right>\\
&=&-\frac1{96}\int re^{-\frac{r}2}(-Y_{1,+1}+Y_{1,-1})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)re^{-\frac{r}2}(-Y_{1,+1}+Y_{1,-1})\,dV\\
&=&-\frac1{96}\int (\frac14r^4-2r^3)e^{-r}\,dr \int\!\!\!\int|-Y_{1,+1}+Y_{1,-1}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{48}\left[\frac14\Gamma(5)-2\Gamma(4)\right]\\
&=&-\frac{6-12}{48}\\
&=&\frac1{8}
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p_x}\left|-\frac1{r}\right|\psi_\mathrm{2p_x}\right>\\
&=&-\frac1{48}\int re^{-r}|-Y_{1,+1}+Y_{1,-1}|^2\,dV\\
&=&-\frac1{48}\int r^3e^{-r}\,dr\int\!\!\!\int |-Y_{1,+1}+Y_{1,-1}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{24}\Gamma(4)\\
&=&-\frac14
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}

\newpage

\subsection{2py}
\begin{eqnarray}
\psi_\mathrm{2p_y}&=&\frac1{\sqrt{32\pi}}re^{-\frac{r}2}\sin \theta\sin \varphi\\
&=&\frac1{\sqrt{24}}re^{-\frac{r}2}\frac{i}{\sqrt 2}\left(Y_{1,+1}+Y_{1,-1}\right)
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2p_y}|\psi_\mathrm{2p_y}\right>&=&\frac1{32\pi}\int r^2e^{-r}\sin^2 \theta \sin^2 \varphi\,dV\\
&=&\frac1{32\pi}\int r^4e^{-r}\,dr\int_0^{\pi}\sin^3 \theta\,d\theta\int_0^{2\pi}\sin^2\varphi\,d\varphi\\
&=&\frac1{32\pi}\Gamma(5)\left[-\cos \theta+\frac13\cos^3 \theta\right]_0^\pi\left[\frac{1-\cos{2\varphi}}2\right]_0^{2\pi}\\
&=&\frac1{32\pi}\times 24\times \frac43 \times \pi\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p_y}\left|-\frac12\nabla^2\right|\psi_\mathrm{2p_y}\right>\\
&=&-\frac1{96}\int re^{-\frac{r}2}(Y_{1,+1}+Y_{1,-1})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)re^{-\frac{r}2}(Y_{1,+1}+Y_{1,-1})\,dV\\
&=&-\frac1{96}\int (\frac14r^4-2r^3)e^{-r}\,dr \int\!\!\!\int|Y_{1,+1}+Y_{1,-1}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{48}\left[\frac14\Gamma(5)-2\Gamma(4)\right]\\
&=&-\frac{6-12}{48}\\
&=&\frac1{8}
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p_y}\left|-\frac1{r}\right|\psi_\mathrm{2p_y}\right>\\
&=&-\frac1{48}\int re^{-r}|Y_{1,+1}+Y_{1,-1}|^2\,dV\\
&=&-\frac1{24}\int r^3e^{-r}\,dr\int\!\!\!\int \frac{|Y_{1,+1}+Y_{1,-1}|^2}2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{24}\Gamma(4)\\
&=&-\frac14
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}

\newpage

\subsection{2pz}
\begin{eqnarray}
\psi_\mathrm{2p_z}&=&\frac1{\sqrt{32\pi}}re^{-\frac{r}2}\cos \theta\\
&=&\frac1{\sqrt{24}}re^{-\frac{r}2}Y_{1,0}
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{2p_z}|\psi_\mathrm{2p_z}\right>&=&\frac1{32\pi}\int r^2e^{-r}\cos^2 \theta\,dV\\
&=&\frac1{32\pi}\int r^4e^{-r}\,dr\int_0^{\pi}\cos^2 \theta\sin \theta\,d\theta\int_0^{2\pi} \,d\varphi\\
&=&\frac1{32\pi}\Gamma(5)\left[-\frac13\cos^3 \theta\right]_0^\pi \times 2\pi\\
&=&\frac1{32\pi}\times 24\times \frac23 \times 2\pi\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p_z}\left|-\frac12\nabla^2\right|\psi_\mathrm{2p_z}\right>\\
&=&-\frac1{48}\int re^{-\frac{r}2}Y_{1,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)re^{-\frac{r}2}Y_{1,0}\,dV\\
&=&-\frac1{48}\int (\frac14r^4-2r^3)e^{-r}\,dr \int\!\!\!\int|Y_{1,0}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{48}\left[\frac14\Gamma(5)-2\Gamma(4)\right]\\
&=&-\frac{6-12}{48}\\
&=&\frac1{8}
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{2p_z}\left|-\frac1{r}\right|\psi_\mathrm{2p_z}\right>\\
&=&-\frac1{24}\int re^{-r}|Y_{1,0}|^2\,dV\\
&=&-\frac1{24}\int r^3e^{-r}\,dr\int\!\!\!\int |Y_{1,0}|^2\sin \theta\,d\theta\,d\varphi\\
&=&-\frac1{24}\Gamma(4)\\
&=&-\frac14
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac18-\frac14\\
&=&-\frac18
\end{eqnarray*}

\newpage

\subsection{3s}
\begin{equation}
\psi_\mathrm{3s}=\frac1{\sqrt{6561\times 3\pi}}\left(27-18r+2r^2\right)e^{-\frac{r}3}
\end{equation}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3s}|\psi_\mathrm{3s}\right>&=&\frac1{6561\times 3\pi}\int \left(27-18r+2r^2\right)^2e^{-\frac{2r}3}\,dV\\
&=&\frac4{6561\times 3}\int \left(729-972r+432r^2-72r^3+4r^4\right)r^2 e^{-\frac{2r}3}\,dr\\
&=&\frac4{6561\times 3}\left[729\Gamma(3)\left(\frac32\right)^3-972\Gamma(4)\left(\frac32\right)^4+432\Gamma(5)\left(\frac32\right)^5-72\Gamma(6)\left(\frac32\right)^6+4\Gamma(7)\left(\frac32\right)^7\right]\\
&=&\frac4{3^9}\left(\frac{3^9}4-\frac{3^{10}}2+3^9\times 4-3^9\times 5+\frac52\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3s}\left|-\frac12\nabla^2\right|\psi_\mathrm{3s}\right>\\
&=&-\frac1{6561\times 6\pi}\int \left(27-18r+2r^2\right)e^{-\frac{r}3}\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)\left(27-18r+2r^2\right)e^{-\frac{r}3}\,dV\\
&=&-\frac2{6561\times 3}\int \left(27-18r+2r^2\right)\left(\frac29r^2-6r+39-\frac{54}{r}\right)e^{-\frac23r} r^2\,dr\\
&=&-\frac2{6561\times 3}\int \left(\frac49r^6-16r^5+192r^4-4\times3^5r^3+25\times 3^4r^2-2\times 3^6r\right)e^{-\frac23r}\,dr\\
&=&-\frac2{6561\times 3}\left[\frac49\Gamma(7)\left(\frac32\right)^7-16\Gamma(6)\left(\frac32\right)^6+192\Gamma(5)\left(\frac32\right)^5\right.\\
&&{}-\left.4\times3^5\Gamma(4)\left(\frac32\right)^4+25\times 3^4\Gamma(3)\left(\frac32\right)^3-2\times 3^6\Gamma(2)\left(\frac32\right)^2\right]\\
&=&-\frac2{6561\times 3}\left(\frac52\times 3^7-10\times 3^7+16\times 3^7-\frac12\times 3^{10}+\frac{25}4\times 3^7-\frac12\times 3^8\right)\\
&=&\frac2{6561\times 3}\times\frac14\times 3^7\\
&=&\frac1{18}
\end{eqnarray*}

\newpage
\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3s}\left|-\frac1{r}\right|\psi_\mathrm{3s}\right>\\
&=&-\frac1{6561\times 3\pi}\int \frac1{r}\left(27-18r+2r^2\right)^2e^{-\frac23r}\,dV\\
&=&-\frac4{6561\times 3}\int \left(729r-972r^2+432r^3-72r^4+4r^5\right)e^{-\frac{2r}3}\,dr\\
&=&\frac4{6561\times 3}\left[729\Gamma(2)\left(\frac32\right)^2-972\Gamma(3)\left(\frac32\right)^3+432\Gamma(4)\left(\frac32\right)^4-72\Gamma(5)\left(\frac32\right)^5+4\Gamma(6)\left(\frac32\right)^6\right]\\
&=&\frac4{3^9}\left(\frac{3^8}4-3^8+2\times 3^8-2\times 3^8+\frac52\times 3^7\right)\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{3px}
\begin{eqnarray}
\psi_\mathrm{3p_x}&=&\frac4{\sqrt{6561\times 4\pi}}r\left(6-r\right)e^{-\frac{r}3}\sin \theta\cos \varphi\\
&=&\frac4{\sqrt{6561\times 6}}r\left(6-r\right)e^{-\frac{r}3}\frac1{\sqrt 2}(-Y_{1,+1}+Y_{1,-1})
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3p_x}|\psi_\mathrm{3p_x}\right>&=&\frac{16}{6561\times 6}\int r^2\left(6-r\right)^2e^{-\frac23r}\frac{|-Y_{1,+1}+Y_{1,-1}|^2}2\,dV\\
&=&\frac{16}{6561\times 6}\int (r^6-12r^5+36r^4)e^{-\frac23r}\,dr\int\!\!\!\int\frac{|-Y_{1,+1}+Y_{1,-1}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 6}\left[\Gamma(7)\left(\frac32\right)^7-12\Gamma(6)\left(\frac32\right)^6+36\Gamma(5)\left(\frac32\right)^5\right]\\
&=&\frac{16}{6561\times 6}\left(\frac58\times 3^9-\frac52\times 3^8+3^8\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p_x}\left|-\frac12\nabla^2\right|\psi_\mathrm{3p_x}\right>\\
&=&-\frac4{6561\times 6}\int r\left(6-r\right)e^{-\frac{r}3}(-Y_{1,+1}+Y_{1,-1})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r\left(6-r\right)e^{-\frac{r}3}(-Y_{1,+1}+Y_{1,-1})\,dV\\
&=&-\frac8{6561\times 6}\int (6r-r^2)\left(-\frac19r^4+\frac83r^3-12r^2\right)e^{-\frac23r}\,dr\int\!\!\!\int\frac{|-Y_{1,+1}+Y_{1,-1}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&-\frac8{6561\times 6}\int\left(\frac19r^6-\frac{10}3r^5+28r^4-72r^3\right)e^{-\frac23r}\,dr\\
&=&-\frac8{6561\times 6}\left[\frac19\Gamma(7)\left(\frac32\right)^7-\frac{10}3\Gamma(6)\left(\frac32\right)^6+28\Gamma(5)\left(\frac32\right)^5-72\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac8{6561\times 6}\left(\frac58\times 3^7-\frac{25}4\times 3^6+7\times 3^6-3^7\right)\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p_x}\left|-\frac1{r}\right|\psi_\mathrm{3p_x}\right>\\
&=&-\frac{16}{6561\times 6}\int r^3(6-r)^2e^{-\frac23r}\,dr\int\!\!\!\int \frac{|-Y_{1,+1}+Y_{1,-1}|^2}2\,d\theta\,d\varphi\\
&=&-\frac{16}{6561\times 6}\int (r^5-12r^4+36r^3)e^{-\frac23r}\,dr\\
&=&-\frac{16}{6561\times 6}\left[\Gamma(6)\left(\frac32\right)^6-12\Gamma(5)\left(\frac32\right)^5+36\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac{16}{6561\times 6}\left(\frac58\times 3^7-3^7+\frac12\times 3^7 \right)\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{3py}
\begin{eqnarray}
\psi_\mathrm{3p_y}&=&\frac4{\sqrt{6561\times 4\pi}}r\left(6-r\right)e^{-\frac{r}3}\sin \theta\sin \varphi\\
&=&\frac4{\sqrt{6561\times 6}}r\left(6-r\right)e^{-\frac{r}3}\frac{i}{\sqrt 2}(Y_{1,+1}+Y_{1,-1})
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3p_y}|\psi_\mathrm{3p_y}\right>&=&\frac{16}{6561\times 6}\int r^2\left(6-r\right)^2e^{-\frac23r}\frac{|Y_{1,+1}+Y_{1,-1}|^2}2\,dV\\
&=&\frac{16}{6561\times 6}\int (r^6-12r^5+36r^4)e^{-\frac23r}\,dr\int\!\!\!\int\frac{|Y_{1,+1}+Y_{1,-1}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 6}\left[\Gamma(7)\left(\frac32\right)^7-12\Gamma(6)\left(\frac32\right)^6+36\Gamma(5)\left(\frac32\right)^5\right]\\
&=&\frac{16}{6561\times 6}\left(\frac58\times 3^9-\frac52\times 3^8+3^8\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p_y}\left|-\frac12\nabla^2\right|\psi_\mathrm{3p_y}\right>\\
&=&-\frac4{6561\times 6}\int r\left(6-r\right)e^{-\frac{r}3}(Y_{1,+1}+Y_{1,-1})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r\left(6-r\right)e^{-\frac{r}3}(Y_{1,+1}+Y_{1,-1})\,dV\\
&=&-\frac8{6561\times 6}\int (6r-r^2)\left(-\frac19r^4+\frac83r^3-12r^2\right)e^{-\frac23r}\,dr\int\!\!\!\int\frac{|Y_{1,+1}+Y_{1,-1}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&-\frac8{6561\times 6}\int\left(\frac19r^6-\frac{10}3r^5+28r^4-72r^3\right)e^{-\frac23r}\,dr\\
&=&-\frac8{6561\times 6}\left[\frac19\Gamma(7)\left(\frac32\right)^7-\frac{10}3\Gamma(6)\left(\frac32\right)^6+28\Gamma(5)\left(\frac32\right)^5-72\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac8{6561\times 6}\left(\frac58\times 3^7-\frac{25}4\times 3^6+7\times 3^6-3^7\right)\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p_y}\left|-\frac1{r}\right|\psi_\mathrm{3p_y}\right>\\
&=&-\frac{16}{6561\times 6}\int r^3(6-r)^2e^{-\frac23r}\,dr\int\!\!\!\int \frac{|Y_{1,+1}+Y_{1,-1}|^2}2\,d\theta\,d\varphi\\
&=&-\frac{16}{6561\times 6}\int (r^5-12r^4+36r^3)e^{-\frac23r}\,dr\\
&=&-\frac{16}{6561\times 6}\left[\Gamma(6)\left(\frac32\right)^6-12\Gamma(5)\left(\frac32\right)^5+36\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac{16}{6561\times 6}\left(\frac58\times 3^7-3^7+\frac12\times 3^7 \right)\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{3pz}
\begin{eqnarray}
\psi_\mathrm{3p_z}&=&\frac4{\sqrt{6561\times 8\pi}}r\left(6-r\right)e^{-\frac{r}3}\cos \theta\\
&=&\frac4{\sqrt{6561\times 6}}r\left(6-r\right)e^{-\frac{r}3}Y_{1,0}
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3p_z}|\psi_\mathrm{3p_z}\right>&=&\frac{16}{6561\times 6}\int r^2\left(6-r\right)^2e^{-\frac23r}|Y_{1,0}|^2\,dV\\
&=&\frac{16}{6561\times 6}\int (r^6-12r^5+36r^4)e^{-\frac23r}\,dr\int\!\!\!\int|Y_{1,0}|^2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 6}\left[\Gamma(7)\left(\frac32\right)^7-12\Gamma(6)\left(\frac32\right)^6+36\Gamma(5)\left(\frac32\right)^5\right]\\
&=&\frac{16}{6561\times 6}\left(\frac58\times 3^9-\frac52\times 3^8+3^8\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p_z}\left|-\frac12\nabla^2\right|\psi_\mathrm{3p_z}\right>\\
&=&-\frac4{6561\times 6}\int r\left(6-r\right)e^{-\frac{r}3}Y_{1,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r\left(6-r\right)e^{-\frac{r}3}Y_{1,0}\,dV\\
&=&-\frac8{6561\times 6}\int (6r-r^2)\left(-\frac19r^4+\frac83r^3-12r^2\right)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,0}|^2\sin \theta \,d\theta\,d\varphi\\
&=&-\frac8{6561\times 6}\int\left(\frac19r^6-\frac{10}3r^5+28r^4-72r^3\right)e^{-\frac23r}\,dr\\
&=&-\frac8{6561\times 6}\left[\frac19\Gamma(7)\left(\frac32\right)^7-\frac{10}3\Gamma(6)\left(\frac32\right)^6+28\Gamma(5)\left(\frac32\right)^5-72\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac8{6561\times 6}\left(\frac58\times 3^7-\frac{25}4\times 3^6+7\times 3^6-3^7\right)\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3p_z}\left|-\frac1{r}\right|\psi_\mathrm{3p_z}\right>\\
&=&-\frac{16}{6561\times 6}\int r^3(6-r)^2e^{-\frac23r}\,dr\int\!\!\!\int |Y_{1,0}|^2\,d\theta\,d\varphi\\
&=&-\frac{16}{6561\times 6}\int (r^5-12r^4+36r^3)e^{-\frac23r}\,dr\\
&=&-\frac{16}{6561\times 6}\left[\Gamma(6)\left(\frac32\right)^6-12\Gamma(5)\left(\frac32\right)^5+36\Gamma(4)\left(\frac32\right)^4\right]\\
&=&-\frac{16}{6561\times 6}\left(\frac58\times 3^7-3^7+\frac12\times 3^7 \right)\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{3dxy}
\begin{eqnarray}
\psi_\mathrm{3d_{xy}}&=&\frac4{\sqrt{6561\times 8\pi}}r^2e^{-\frac{r}3}\sin^2 \theta\sin \varphi\cos \varphi\\
&=&\frac{4}{\sqrt{6561\times 30}}r^2e^{-\frac{r}3}\frac{i}{\sqrt 2}(Y_{2,+2}-Y_{2,-2})
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3d_{xy}}|\psi_\mathrm{3d_{xy}}\right>&=&\frac{16}{6561\times 30}\int r^4e^{-\frac23r}\frac{|Y_{2,+2}-Y_{2,-2}|^2}2\,dV\\
&=&\frac{16}{6561\times 30}\int r^6e^{-\frac23r}\,dr\int\!\!\!\int\frac{|Y_{2,+2}-Y_{2,-2}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 30}\Gamma(7)\left(\frac32\right)^7\\
&=&\frac{16}{6561\times 30}\left(\frac58\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{xy}}\left|-\frac12\nabla^2\right|\psi_\mathrm{3d_{xy}}\right>\\
&=&-\frac4{6561\times 30}\int r^2e^{-\frac{r}3}(Y_{2,+2}-Y_{2,-2})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r^2e^{-\frac{r}3}(Y_{2,+2}-Y_{2,-2})\,dV\\
&=&-\frac8{6561\times 30}\int \left(\frac19r^6-2r^5\right)e^{-\frac23r}\,dr\int\!\!\!\int \frac{|Y_{2,+2}-Y_{2,-2}|^2}2\,d\theta d\varphi\\
&=&-\frac8{6561\times 30}\left[\frac19\Gamma(7)\left(\frac32\right)^7-2\Gamma(6)\left(\frac32\right)^6\right]\\
&=&-\frac8{6561\times 30}\left(\frac58\times 3^7-\frac54\times 3^7\right)\\
&=&\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{xy}}\left|-\frac1{r}\right|\psi_\mathrm{3d_{xy}}\right>\\
&=&-\frac{16}{6561\times 30}\int r^3e^{-\frac23r}\frac{|Y_{2,+2}-Y_{2,-2}|^2}2\,dV\\
&=&-\frac{16}{6561\times 30}\int r^5e^{-\frac23r}\,dr\int\!\!\!\int \frac{|Y_{2,+2}-Y_{2,-2}|^2}2\sin \theta\,d\theta d\varphi\\
&=&-\frac{16}{6561\times 30}\Gamma(6)\left(\frac32\right)^6\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{3dyz}
\begin{eqnarray}
\psi_\mathrm{3d_{yz}}&=&\frac4{\sqrt{6561\times 8\pi}}r^2e^{-\frac{r}3}\sin \theta\cos \theta\sin \varphi\\
&=&\frac{4}{\sqrt{6561\times 30}}r^2e^{-\frac{r}3}\frac{i}{\sqrt 2}(-Y_{2,+1}-Y_{2,-1})
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3d_{yz}}|\psi_\mathrm{3d_{yz}}\right>&=&\frac{16}{6561\times 30}\int r^4e^{-\frac23r}\frac{|-Y_{2,+1}-Y_{2,-1}|^2}2\,dV\\
&=&\frac{16}{6561\times 30}\int r^6e^{-\frac23r}\,dr\int\!\!\!\int\frac{|-Y_{2,+1}-Y_{2,-1}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 30}\Gamma(7)\left(\frac32\right)^7\\
&=&\frac{16}{6561\times 30}\left(\frac58\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{yz}}\left|-\frac12\nabla^2\right|\psi_\mathrm{3d_{yz}}\right>\\
&=&-\frac4{6561\times 30}\int r^2e^{-\frac{r}3}(-Y_{2,+1}-Y_{2,-1})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r^2e^{-\frac{r}3}(-Y_{2,+1}-Y_{2,-1})\,dV\\
&=&-\frac8{6561\times 30}\int \left(\frac19r^6-2r^5\right)e^{-\frac23r}\,dr\int\!\!\!\int \frac{|-Y_{2,+1}-Y_{2,-2}|^2}2\,d\theta d\varphi\\
&=&-\frac8{6561\times 30}\left[\frac19\Gamma(7)\left(\frac32\right)^7-2\Gamma(6)\left(\frac32\right)^6\right]\\
&=&-\frac8{6561\times 30}\left(\frac58\times 3^7-\frac54\times 3^7\right)\\
&=&\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{yz}}\left|-\frac1{r}\right|\psi_\mathrm{3d_{yz}}\right>\\
&=&-\frac{16}{6561\times 30}\int r^3e^{-\frac23r}\frac{|-Y_{2,+1}-Y_{2,-1}|^2}2\,dV\\
&=&-\frac{16}{6561\times 30}\int r^5e^{-\frac23r}\,dr\int\!\!\!\int \frac{|-Y_{2,+1}-Y_{2,-1}|^2}2\sin \theta\,d\theta d\varphi\\
&=&-\frac{16}{6561\times 30}\Gamma(6)\left(\frac32\right)^6\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{3dzx}
\begin{eqnarray}
\psi_\mathrm{3d_{zx}}&=&\frac4{\sqrt{6561\times 8\pi}}r^2e^{-\frac{r}3}\sin \theta\cos \theta\cos \varphi\\
&=&\frac{4}{\sqrt{6561\times 30}}r^2e^{-\frac{r}3}\frac1{\sqrt 2}(-Y_{2,+1}+Y_{2,-1})
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3d_{zx}}|\psi_\mathrm{3d_{zx}}\right>&=&\frac{16}{6561\times 30}\int r^4e^{-\frac23r}\frac{|-Y_{2,+1}+Y_{2,-1}|^2}2\,dV\\
&=&\frac{16}{6561\times 30}\int r^6e^{-\frac23r}\,dr\int\!\!\!\int\frac{|-Y_{2,+1}+Y_{2,-1}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 30}\Gamma(7)\left(\frac32\right)^7\\
&=&\frac{16}{6561\times 30}\left(\frac58\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{zx}}\left|-\frac12\nabla^2\right|\psi_\mathrm{3d_{zx}}\right>\\
&=&-\frac4{6561\times 30}\int r^2e^{-\frac{r}3}(-Y_{2,+1}+Y_{2,-1})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r^2e^{-\frac{r}3}(-Y_{2,+1}+Y_{2,-1})\,dV\\
&=&-\frac8{6561\times 30}\int \left(\frac19r^6-2r^5\right)e^{-\frac23r}\,dr\int\!\!\!\int \frac{|-Y_{2,+1}+Y_{2,-2}|^2}2\,d\theta d\varphi\\
&=&-\frac8{6561\times 30}\left[\frac19\Gamma(7)\left(\frac32\right)^7-2\Gamma(6)\left(\frac32\right)^6\right]\\
&=&-\frac8{6561\times 30}\left(\frac58\times 3^7-\frac54\times 3^7\right)\\
&=&\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{zx}}\left|-\frac1{r}\right|\psi_\mathrm{3d_{zx}}\right>\\
&=&-\frac{16}{6561\times 30}\int r^3e^{-\frac23r}\frac{|-Y_{2,+1}+Y_{2,-1}|^2}2\,dV\\
&=&-\frac{16}{6561\times 30}\int r^5e^{-\frac23r}\,dr\int\!\!\!\int \frac{|-Y_{2,+1}+Y_{2,-1}|^2}2\sin \theta\,d\theta d\varphi\\
&=&-\frac{16}{6561\times 30}\Gamma(6)\left(\frac32\right)^6\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage



\subsection{3dx${}^2$-y${}^2$}
\begin{eqnarray}
\psi_\mathrm{3d_{x^2-y^2}}&=&\frac4{\sqrt{6561\times 8\pi}}r^2e^{-\frac{r}3}\sin^2 \theta(\cos^2 \varphi-\sin^2\varphi)\\
&=&\frac{4}{\sqrt{6561\times 30}}r^2e^{-\frac{r}3}\frac1{\sqrt 2}(Y_{2,+2}+Y_{2,-2})
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3d_{x^2-y^2}}|\psi_\mathrm{3d_{x^2-y^2}}\right>&=&\frac{16}{6561\times 30}\int r^4e^{-\frac23r}\frac{|Y_{2,+2}+Y_{2,-2}|^2}2\,dV\\
&=&\frac{16}{6561\times 30}\int r^6e^{-\frac23r}\,dr\int\!\!\!\int\frac{|Y_{2,+2}+Y_{2,-2}|^2}2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 30}\Gamma(7)\left(\frac32\right)^7\\
&=&\frac{16}{6561\times 30}\left(\frac58\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{xy}}\left|-\frac12\nabla^2\right|\psi_\mathrm{3d_{xy}}\right>\\
&=&-\frac4{6561\times 30}\int r^2e^{-\frac{r}3}(Y_{2,+2}+Y_{2,-2})^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r^2e^{-\frac{r}3}(Y_{2,+2}+Y_{2,-2})\,dV\\
&=&-\frac8{6561\times 30}\int \left(\frac19r^6-2r^5\right)e^{-\frac23r}\,dr\int\!\!\!\int \frac{|Y_{2,+2}+Y_{2,-2}|^2}2\,d\theta d\varphi\\
&=&-\frac8{6561\times 30}\left[\frac19\Gamma(7)\left(\frac32\right)^7-2\Gamma(6)\left(\frac32\right)^6\right]\\
&=&-\frac8{6561\times 30}\left(\frac58\times 3^7-\frac54\times 3^7\right)\\
&=&\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{xy}}\left|-\frac1{r}\right|\psi_\mathrm{3d_{xy}}\right>\\
&=&-\frac{16}{6561\times 30}\int r^3e^{-\frac23r}\frac{|Y_{2,+2}+Y_{2,-2}|^2}2\,dV\\
&=&-\frac{16}{6561\times 30}\int r^5e^{-\frac23r}\,dr\int\!\!\!\int \frac{|Y_{2,+2}+Y_{2,-2}|^2}2\sin \theta\,d\theta d\varphi\\
&=&-\frac{16}{6561\times 30}\Gamma(6)\left(\frac32\right)^6\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\newpage

\subsection{3d3x${}^2$-r${}^2$}
\begin{eqnarray}
\psi_\mathrm{3d_{3z^2-r^2}}&=&\frac1{\sqrt{6561\times 6\pi}}r^2e^{-\frac{r}3}(3\cos^2 \theta-1)\\
&=&\frac{4}{\sqrt{6561\times 30}}r^2e^{-\frac{r}3}Y_{2,0}
\end{eqnarray}

\subsubsection{normalization}
\begin{eqnarray*}
\left<\psi_\mathrm{3d_{3z^2-r^2}}|\psi_\mathrm{3d_{3z^2-r^2}}\right>&=&\frac{16}{6561\times 30}\int r^4e^{-\frac23r}|Y_{2,0}|^2\,dV\\
&=&\frac{16}{6561\times 30}\int r^6e^{-\frac23r}\,dr\int\!\!\!\int|Y_{2,0}|^2\sin \theta \,d\theta\,d\varphi\\
&=&\frac{16}{6561\times 30}\Gamma(7)\left(\frac32\right)^7\\
&=&\frac{16}{6561\times 30}\left(\frac58\times 3^9\right)\\
&=&1
\end{eqnarray*}

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{3z^2-r^2}}\left|-\frac12\nabla^2\right|\psi_\mathrm{3d_{3z^2-r^2}}\right>\\
&=&-\frac4{6561\times 30}\int r^2e^{-\frac{r}3}Y_{2,0}^\ast\left(\frac{\partial^2}{\partial r^2}+\frac2{r}\frac{\partial}{\partial r}+\frac1{r^2}\Lambda\right)r^2e^{-\frac{r}3}Y_{2,0}\,dV\\
&=&-\frac8{6561\times 30}\int \left(\frac19r^6-2r^5\right)e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,0}|^2\,d\theta d\varphi\\
&=&-\frac8{6561\times 30}\left[\frac19\Gamma(7)\left(\frac32\right)^7-2\Gamma(6)\left(\frac32\right)^6\right]\\
&=&-\frac8{6561\times 30}\left(\frac58\times 3^7-\frac54\times 3^7\right)\\
&=&\frac1{18}
\end{eqnarray*}

\newpage

\subsubsection{$\left$}

\begin{eqnarray*}
\left&\equiv&\left<\psi_\mathrm{3d_{3z^2-r^2}}\left|-\frac1{r}\right|\psi_\mathrm{3d_{3z^2-r^2}}\right>\\
&=&-\frac{16}{6561\times 30}\int r^3e^{-\frac23r}|Y_{2,0}|^2\,dV\\
&=&-\frac{16}{6561\times 30}\int r^5e^{-\frac23r}\,dr\int\!\!\!\int |Y_{2,0}|^2\sin \theta\,d\theta d\varphi\\
&=&-\frac{16}{6561\times 30}\Gamma(6)\left(\frac32\right)^6\\
&=&-\frac19
\end{eqnarray*}

\subsubsection{$\left$}
\begin{eqnarray*}
\left&=&\left+\left\\
&=&\frac1{18}-\frac19\\
&=&-\frac1{18}
\end{eqnarray*}

\end{document}

National Center For University Entrance Examinations

Link

• National Center For University Entrance Examinations
  
  http://www.dnc.ac.jp/ [another window]
  



• Kei-net
  
  http://hiw.oo.kawai-juku.ac.jp/nyushi/ [another window]
  

---

COMMENT

• 数学IA
  
  とても簡単かつ典型的な問題のみ。困るところはなかったように思える。
  
• 数学IIB
  
  とても典型的。しかしながら積分の計算は時間がかかる。時間配分をうまく考えて、十分解ききってほしい問題である。