楕円軌道と内心の軌跡

楕円の軌跡は焦点からの線分の長さの和が一定という条件で描くことができる。楕円の長半径を$a$，短半径を$b$，焦点の座標を${\rm O}_1=(-c,0)$，${\rm O}_2=(c,0)$，楕円上の点Pの座標を$(x,y)$とする。例えば，$\ell_1={\rm O}_1{\rm P}$，$\ell_2={\rm O}_2{\rm P}$として，$\ell_1+\ell_2=2 a$と一定となる。

このとき，三角形${\rm O}_1{\rm O}_2 P$の内心（内接円の中心）Qの軌跡はどんな図形を描くだろうか。twitterでアニメーションをみかけたが，楕円に見えたので確かめてみよう。点Ｐは次の楕円の方程式の上を動く。
\begin{equation}
\dfrac{x^2}{a^2} + \dfrac{y^2}{b^2} = 1 \quad a^2=b^2+c^2
\end{equation}
このとき，$\ell_1, \ell_2$を求めてみる。
\begin{equation}
\begin{aligned}
\ell_1 &= \sqrt{ (x+c)^2 + y^2 } = \sqrt{(x+c)^2 + b^2 (1 - (x/a)^2 ) }\\
&= \sqrt{(x+c)^2 + (1-c^2/a^2) (a^2 - x^2) } = a + \dfrac{c x}{a}\\
\ell_2 &= \sqrt{ (x-c)^2 + y^2 } = \sqrt{(x-c)^2 + b^2 (1 - (x/a)^2 ) } \\
&= \sqrt{(x-c)^2 + (1-c^2/a^2) (a^2 - x^2) } = a - \dfrac{c x}{a}
\end{aligned}
\end{equation}
次に，内心Qの座標を，$(p,q)$とする。$q$は内接円の半径と等しい。三角形の内接円の半径$r$は，三角形の面積$S$と$2S=(\ell_1+\ell_2+2c) r$の関係がある。ヘロンの公式より，
\begin{equation}
\begin{aligned}
s &= (\ell_1+\ell_2+2c)/2 = a+c \\
S &=\sqrt{s(s-2c)(s-\ell_1)(s-\ell_2)}=\sqrt{(a+c)(a-c)(c - c x/a)(c + c x/a)}\\
\therefore q &= r= \dfrac{S}{a+c}=c\sqrt{\frac{a-c}{a+c}(1-x^2/a^2)}= \dfrac{c y}{b}\sqrt{\frac{a-c}{a+c}}\equiv  \dfrac{c\ y\ \varepsilon}{b}
\end{aligned}
\end{equation}
また，角${\rm PO_1 O_2}=\phi$，角${\rm PO_2 O_1}=\theta$とすると，余弦定理から，
\begin{equation}
\begin{aligned}
\cos\phi = \dfrac{\ell_1^2+(2c)^2-\ell_2^2}{4\ell_1 c}=\dfrac{x+c}{\ell_1}\\
\cos\theta = \dfrac{\ell_2^2+(2c)^2-\ell_1^2}{4\ell_1 c}=\dfrac{c-x}{\ell_2}
\end{aligned}
\end{equation}
内心の性質から，角${\rm QO_1 O_2}=\phi/2$，角${\rm QO_2 O_1}=\theta/2$であり，半角の公式から，
\begin{equation}
\begin{aligned}
\tan{\frac{\phi}{2}} &=\sqrt{\dfrac{1-\cos\phi}{1+\cos\phi}} =\sqrt{\dfrac{\ell_1-(x+c)}{\ell_1+(x+c)} }\\
&=\sqrt{\dfrac{a+cx/a-(x+c)}{a+cx/a+(x+c)}} = \varepsilon \sqrt{\dfrac{a-x}{a+x}}\\
\tan{\frac{\theta}{2}} &=\sqrt{\dfrac{1-\cos\theta}{1+\cos\theta}}=\sqrt{\dfrac{\ell_2-(c-x)}{\ell_2+(c-x)} }\\
&=\sqrt{\dfrac{a-cx/a-(c-x)}{a-cx/a+(c-x)}} = \varepsilon \sqrt{\dfrac{a+x}{a-x}}\\
\end{aligned}
\end{equation}
${\rm O_1}$から角度$\phi/2$で望む内心Qのy座標が，${\rm O_2}$から角度$\theta/2$で望むものと等しいことから，
\begin{equation}
\begin{aligned}
(p+c) \tan \dfrac{\phi}{2} &= (c-p) \tan \dfrac{\theta}{2} \\
p \bigl( \tan \dfrac{\phi}{2} +  \tan \dfrac{\theta}{2} \bigr) &= c \bigl(  \tan \dfrac{\theta}{2} - \tan \dfrac{\phi}{2} \bigr) \\
p \varepsilon \Bigl( \sqrt{\dfrac{a-x}{a+x}} + \sqrt{\dfrac{a+x}{a-x}} \Bigr)
& = c \varepsilon  \Bigl( \sqrt{\dfrac{a+x}{a-x}} - \sqrt{\frac{a-x}{a+x}} \Bigr) \\
2 a p \varepsilon &= 2 x c \varepsilon\\
\therefore p &= \dfrac{c}{a} x
\end{aligned}
\end{equation}
これから，内心Ｑの満足する軌跡の方程式は，Pが描く楕円の軌跡の式を用いて以下のように求まった。
\begin{equation}
\dfrac{p^2}{c^2} + \dfrac{q^2}{(c \varepsilon)^2} = 1
\end{equation}

図　楕円の軌跡（青）と内心の軌跡（赤）

公転速度と公転周期

万有引力定数は，G=6.7 × 10^-11 m ^3 kg^-1 s^-2 である。高等学校の物理の教科書にのっているように，質量 M kg の天体の周りを質量 m kg の天体が速さ v m/s の速度で半径 R m の等速円運動するとき，$ \dfrac{G M m }{R^2} = \dfrac{m v^2}{R}$ から，$ v= \sqrt{\dfrac{GM}{R}} $，周期は$ T = \dfrac{2 \pi R }{v}$であった。これで，ブラックホール（BH）のまわりの"惑星"の公転速度と公転周期が求まる。だからどうしたといわれても。

地球−月（M = 6.0 × 10^24 kg, R = 3.8 × 10^8 m）→（v= 1 km/s, T= 0.76 y）
太陽−地球（M = 2.0 × 10^30 kg, R = 1.5 × 10^11 m）→（v= 30 km/s, T= 1.0 y）
銀河中心-太陽（M = 2.2 × 10^41 kg，R= 2.6 × 10^20 m）→（v= 240 km/s, T= 2.0 × 10^8 y）
BH−"惑星"（M = 8.2 × 10^36 kg，R = 1.0 ×10^17 m）→（v=74 km/s, T=2.7 × 10^5 y）

［１］国立天文台，最新の観測による銀河中心〜太陽系の距離や回転速度を発表
［２］超大質量ブラックホール（Wikipedia）

ブラックホールの周りの惑星

物理ではなかなかびっくりする話がない。いや，あるのだが，驚嘆するためにはそれなりの基礎知識が必要なのでたいへんだ。数学でもそうかもしれない。その意味では望月新一さんのIUTはとても大きなトピックだった。一方，天文学では結構な頻度で，驚きのニュースが飛び込んでくるような気がする。最近のそれは，ブラックホールの周りの惑星だ。

ブラックホールの周りには降着円盤があることは，福江先生の得意分野でもあるので知っていた。そこに原始太陽系の形成モデルのシミュレーションを適用すると，銀河中心にあるような太陽の1000万倍の質量を持つ巨大ブラックホールのまわりのミクロな塵から，数億年をかけて，ブラックホール中心から10光年ほどのところに，地球質量の10倍程度の惑星が１万個以上形成されるというものだ。となりの惑星との距離は0.1光年のオーダーかな。生命が誕生するような熱源は確保できるのだろうか。ＳＦ作家の夢が広がるだろうか。

［１］K. Wada, Y. Tsukamoto and E. Kokubo,  Planet Formation around Supermassive Black Holes in the Active Galactic Nuclei（2019.11.26）

児童生徒１人１台ＰＣ

11月13日の経済財政諮問会議で，経済対策（未来投資）として義務教育の児童生徒１人１台のPCを配備するという考えが示された。経済再生担当大臣の西村康稔は，1995年から1997年にかけて通産省からの出向で石川県の商工課長を務めていた。そのころ父が亡くなり，葬儀後の挨拶で県庁もまわった際に名刺だけを置いてきた。西村は2003年に衆議院議員になったが，その前後からしばらく後援会からの案内が届いていたことがあった。選挙区は離れているのであまり効果はないと思うのだが，とりあえず関西圏なので。

さて，児童生徒１人１台PCに話を戻す。これが，小学校５年生から中学３年生までとなると，小学校 5, 6 年が 213万人程度，中学校1, 2, 3 年が 322万人なので，合計540万人。現在あるPCの普及率は 5.4 人に1台なので，人数として現有のもので必要数の 19% 程度がカバーされているとするなら，540万 × 81%で 440万台新たに必要になる（現有機器がデスクトップならばもっと必要かもしれない）。１台10万円として，4400億円だ（いつまでに実行するのだろう？）。その後，更新を続けるとすれば，毎年数百億円が必要になる。で，これを学校で管理できるのだろうか。活用できるのだろうか。どこかでBOYDにスライドできるのだろうか。謎は深まるが，あまりまともな制度設計がされることは期待できないように思う。

P. S. 読売新聞によると，2022年までに小５〜中３，2024年までに小１〜小４という説もあるようだ。

P. P. S. 12/3 PISA2018結果発表を前に，日経新聞によると2023年が目標年度であり，5000億円を児童生徒１人１台のPCまたはタブレットにつっこむそうだ。初年度は1500億円とか。南無阿弥陀仏。

山羊問題

山羊問題（Goat Problem）は次のような問題である。中心Ａ，半径$r=1$の草の生えた円形の土地Ｓがある。その周の１点Ｏから長さ$a$のヒモにつながれた山羊を放し飼いにすると，Ｓのうち，半径$a$のOを中心とし半径$a$の円内の草が食べられてしまう。その面積がＳの半分になるようなヒモの長さ$a$はいくらか？図形はOAを結ぶ線に対称なので，半分だけ考えてみよう。

図　山羊問題

山羊が食べた草地の面積は，扇型O-BP+扇型A-OP-三角形AOPであり，これが$\pi r^2/4$になればよい。角POA=$\phi$とすると，角PAO=$\pi-2\phi$である。これを式で表すと，

\begin{equation}
\begin{aligned}
\dfrac{\pi r^2}{4} &= \dfrac{r^2}{2}(\pi - 2 \phi) + \dfrac{a^2}{2} \phi - \dfrac{r^2}{2} \sin 2 \phi \\
a &= 2 r \cos \phi
\end{aligned}
\end{equation}
したがって，
\begin{equation}
4 \phi \cos^2 \phi = \sin 2 \phi +2 \phi -\dfrac{\pi}{2}
\end{equation}
これをMathematicaで解くと，$a/r = 1.15873$となった。
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
In[1]:= Clear[x]; sol1 = FindRoot[4 Cos[x]^2  x == Sin[2 x] + 2 x - Pi/2, {x, 1}]; x = x /. sol1
Out[1]= 0.952848
In[2]:= a = 2 Cos[x]
Out[2]= 1.15873
In[3]:= Clear[x]; sol2 = FindRoot[Sqrt[1 - (x - 1)^2] == Sqrt[a^2 - x^2], {x, 1}]; b = x /. sol2
Out[3]= 0.671326
In[4]:= c = NIntegrate[Sqrt[1 - (x - 1)^2], {x, 0, b}] + NIntegrate[Sqrt[a^2 - x^2], {x, b, a}]
Out[4]= 0.785398

In[5]:= g1 = Plot[{Sqrt[a^2 - x^2], Sqrt[1 - (x - 1)^2],
   Sqrt[a^2 - b^2]/b x, -Sqrt[a^2 - b^2]/(a - b) (x - a), -Sqrt[a^2 - b^2]/(1 - b) (x - 1)}, {x, 0, 2},
  AspectRatio -> Automatic, PlotStyle -> {, , Dashed, Dashed, Dotted}, PlotRange -> {0, 1.2}]
In[6]:= g2 = Graphics[{Text[O, {0.05, 0.03}], Text[A, {0.95, 0.03}],
   Text[B, {1.20, 0.03}], Text[P, {0.69, 1.0}]}];
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

パラメータ励振（５）

（パラメータ励振（４）からの続き）

\begin{equation}
\begin{aligned}
\ddot{x}+\dfrac{g-\ddot{\ell}}{\ell}x &= 0\\
\ell &= \ell_0(1-\epsilon \cos(\omega t + \delta))\\
g/\ell_0 &= \omega_0^2\\
\end{aligned}
\end{equation}
このブランコのモデルで，エネルギーが励振項から供給される様子を調べてみる。
微分方程式を変形すると次のようになる。
\begin{equation}
\ddot{x}+\omega_0^2 x = - \epsilon(\omega_0^2-\omega^2) \cos (\omega t + \delta ) x
\end{equation}
$\dot{x}$を微分方程式の両辺にかけると，
\begin{equation}
\dot{x} \ddot{x}+\omega_0^2 \dot{x} x = - \epsilon(\omega_0^2-\omega^2) \cos (\omega t + \delta ) \dot{x} x
\end{equation}
したがって，
\begin{equation}
\dfrac{d}{dt} (\dot{x}^2+\omega_0^2  x^2 )  = - \epsilon(\omega_0^2-\omega^2) \cos (\omega t + \delta ) \dot{x} x \equiv {\rm R}
\end{equation}
つまり，左辺は運動エネルギー＋位置エネルギーの微分となるため，右辺がこれに対するエネルギーソースの働きをしている。そこで，右辺Rの符号を調べる。

ここで，$\omega=2\omega_0, \lambda = \dfrac{-3\epsilon \omega_0^2}{4\pi} < 0$とする。また，$\delta =(0, \pi/2, \pi, 3\pi/2)$に対して，$\varphi = (\pi/4, \pi/2, \pi/4, 0)$に注意する。
\begin{equation}
\begin{aligned}
x &= C_1 e^{\lambda t} \cos(\omega_0 t + \varphi) + C_2 e^{-\lambda t} \sin (\omega_0 t + \varphi)\\
\dot{x} &= \lambda \{ C_1 e^{\lambda t} \cos(\omega_0 t + \varphi) - C_2 e^{-\lambda t} \sin (\omega_0 t + \varphi) \} \\
&+\omega_0 \{ -C_1 e^{\lambda t} \sin(\omega_0 t + \varphi) + C_2 e^{-\lambda t} \cos (\omega_0 t + \varphi) \}
\end{aligned}
\end{equation}
初期条件として，$x(0)=c >0, \dot{x}(0)=0$とする。
\begin{equation}
\begin{aligned}
c &= C_1  \cos \varphi + C_2 \sin \varphi\\
0 &= \lambda \{ C_1 \cos \varphi - C_2 \sin \varphi \} \\
&+\omega_0 \{ -C_1 \sin \varphi + C_2 \cos \varphi \}
\end{aligned}
\end{equation}
これを解いて，
\begin{equation}
\begin{aligned}
C_1 &= c \dfrac{\lambda \sin \varphi - \omega_0 \cos \varphi}{\lambda \sin 2\varphi -\omega_0}\\
C_2 &= c \dfrac{\lambda \cos \varphi - \omega_0 \sin \varphi}{\lambda \sin 2\varphi -\omega_0}
\end{aligned}
\end{equation}
近似解は，$x(t)=c cos(\omega_0 t)$であるから，$\dot{x} \ x = -\dfrac{c^2}{2} \sin(2\omega_0 t) $ である。したがって，$\delta=\pi/2の場合$
\begin{equation}
{\rm R}= \epsilon(\omega_0^2-\omega^2) \cos (2\omega_0 t + \delta ) \dfrac{c^2}{2} \sin(2 \omega_0 t) = \dfrac{3 \omega_0^2 c^2}{2} \sin^2 (2\omega_0 t) > 0
\end{equation}
となって，エネルギーが増加することがわかる。また，$\delta = 3\pi/2$では符号が逆転する。

パラメータ励振（４）

以下のブランコのモデルに対する解析的な近似解を考える。
\begin{equation}
\begin{aligned}
\ddot{x}+\dfrac{g-\ddot{\ell}}{\ell}x &= 0\\
\ell &= \ell_0(1-\epsilon \cos(\omega t + \delta))\\
g/\ell_0 &= \omega_0^2\\
\end{aligned}
\end{equation}
もとの微分方程式を $x(t), y(t)$ の１階連立微分方程式の形に表す。
\begin{equation}
\begin{aligned}
\dot{x} &= y\\
\dot{y} &= -\{\omega_0^2+\epsilon(\omega_0^2-\omega^2)\ \cos(\omega t + \delta)\}\ x
\end{aligned}
\end{equation}
ここで，$x=a(t) \cos (\omega_0 t + \phi(t)), y= -a \omega_0 \sin (\omega_0 t + \phi(t))$とおいて，
上の連立微分方程式を$a(t), \phi(t) $の連立微分方程式に書き直す。
\begin{equation}
\begin{aligned}
\dot{a} \cos (\omega_0 t + \phi) -a \dot{\phi} \sin (\omega_0 t + \phi) &= 0\\
\dot{a} \sin (\omega_0 t + \phi) +a \dot{\phi} \cos (\omega_0 t + \phi) &= \dfrac{ \epsilon(\omega_0^2-\omega^2)}{\omega_0}\ \cos(\omega t + \delta) a \cos(\omega_0 t + \phi )\\
\end{aligned}
\end{equation}
整理すると次のような２式となる。
\begin{equation}
\begin{aligned}
\dot{a} &= \dfrac{ \epsilon(\omega_0^2-\omega^2)}{2\omega_0}\
\cos(\omega t + \delta) \sin(2 \omega_0 t + 2 \phi ) \ a \\
\dot{\phi} &= \dfrac{ \epsilon(\omega_0^2-\omega^2)}{2\omega_0}\
\cos(\omega t + \delta) \{ 1 + \cos(2 \omega_0 t + 2 \phi ) \} \\
\end{aligned}
\end{equation}
ここで，$\omega = 2 \omega_0$とし，$a(t),\phi(t)$の時間変化が緩いとして上式の右辺を$t=0$から周期$T=2 \pi/\omega$まで時間で積分した量を周期で割った量で置き換える。
\begin{equation}
\begin{aligned}
\dot{a}  &= \dfrac{\omega}{2\pi} \int_0^{2\pi/\omega} \dfrac{\epsilon (\omega_0^2-\omega^2)}{2 \omega_0} \cos(\omega t + \delta) \sin(2 \omega_0 t + 2 \phi ) \ a dt\\
\dot{\phi} &=  \dfrac{\omega}{2\pi} \int_0^{2\pi/\omega} \dfrac{\epsilon (\omega_0^2-\omega^2)}{2 \omega_0} \cos(\omega t + \delta) \{ 1 + \cos(2 \omega_0 t + 2 \phi ) \} dt
\end{aligned}
\end{equation}
このようにして平均化された$\dot{a},\dot{\phi}$に対して次式が成り立つ。
\begin{equation}
\begin{aligned}
\lambda &= \dfrac{\omega}{2\pi}\dfrac{\epsilon (\omega_0^2-\omega^2)}{4 \omega_0}\\
\dot{a} &= - \lambda \sin (\delta - 2 \phi) \ a \\
\dot{\phi} &= \lambda \cos (\delta - 2 \phi) \\
\end{aligned}
\end{equation}
さらに，$u=a \cos \phi, v= a \sin \phi$と置き上式を代入して加法定理を用い，さらに整理すると，
\begin{equation}
\begin{aligned}
\dot{u} &= \dot{a} \cos \phi - a \dot{\phi} \sin \phi = -\lambda a \sin (\delta - \phi)
= -\lambda (u \sin \delta - v \cos \delta)\\
\dot{v} &= \dot{a} \sin \phi + a \dot{\phi} \cos \phi = \lambda a \cos (\delta - \phi)
= \lambda (u \cos \delta + v \sin \delta)
\end{aligned}
\end{equation}
結局，
\begin{equation}
\dfrac{d}{dt}\begin{pmatrix} u \\ v \end{pmatrix}
= \lambda \begin{pmatrix} - \sin \delta & \cos \delta \\ \cos \delta & \sin \delta \end{pmatrix} \begin{pmatrix} u \\ v \end{pmatrix}
\end{equation}
さらに，$(u,v)=(A,B)e^{p \lambda t}$と置くと，
\begin{equation}
\begin{pmatrix} p+ \sin \delta & -\cos \delta \\ -\cos \delta & p-\sin \delta \end{pmatrix} \begin{pmatrix} A \\ B \end{pmatrix}=0
\end{equation}
自明でない解を持つ条件から，$p^2-\sin^2 \delta -\cos^2 \delta =0$より，$p=\pm 1$
$p=1$の場合
\begin{equation}
(1+\sin \delta)A -\cos\delta B=0 \quad \therefore (A,B)=(\sqrt{\dfrac{1-\sin \delta}{2}}, \sqrt{\dfrac{1+\sin \delta}{2}})\\
a=e^{\lambda t}, \quad \cos\phi = \sqrt{\dfrac{1-\sin \delta}{2}}, \quad \sin \phi = \sqrt{\dfrac{1+\sin \delta}{2}}
\end{equation}
$p=-1$の場合
\begin{equation}
(-1+\sin \delta)A -\cos\delta B=0 \quad \therefore (A,B)=(\sqrt{\dfrac{1+\sin \delta}{2}},
-\sqrt{\dfrac{1-\sin \delta}{2}})\\
a=e^{-\lambda t}, \quad \cos\phi = \sqrt{\dfrac{1+\sin \delta}{2}}, \quad \sin \phi =
 -\sqrt{\dfrac{1-\sin \delta}{2}}
\end{equation}
したがって，$x=a \cos (\omega_0 t + \phi)$の一般解は，２つのモードの重ね合わせとして表現される。
\begin{equation}
x(t)=\ell(t) \theta(t) =  C_1 e^{\lambda t} \cos (\omega_0 t+ \varphi)
+ C_2 e^{-\lambda t} \sin (\omega_0 t+ \varphi)\\
\varphi = \tan^{-1}\sqrt{\dfrac{1+\sin \delta}{1-\sin \delta}}
\end{equation}
$\omega = 2\omega_0$としたので，$\lambda = -\dfrac{3 \epsilon \omega_0^2}{4 \pi} < 0$である。$\delta=\pi/2$の場合，指数関数的に増大する項は，$C_2 e^{-\lambda t} \cos(\omega_0 t)$という形になる。

参考文献
［１］対話・非線形振動（楽しい物理ノート KENZOU）
［２］Parametirc Oscillator（Wikipedia en:）

パラメータ励振（３）

ブランコのモデルでは，重心の位置が振動の両端で最も高くなるのが良いのだと思い込んでいた。そのため，前回のモデルでは，$\epsilon=a/\ell_0 > 0$として，$\ell=\ell_0 (1-\epsilon \cos \omega t)$と考えた。つまり，$t=0$と$t=2\pi/\omega$で，$\ell=\ell_0 (1-\epsilon)$と重心が高くなり，$t=\pi/\omega$で，$\ell=\ell_0 (1+\epsilon)$と重心が低くなるわけだ。

しかし，どうやらそうではなかった。$\ell=\ell_0 (1-\epsilon \cos (\omega t + \pi/2))$のときに，最も励振が大きくなるのだ。前半の1/4周期に重心の落ちる速度が正，後半の1/4周期に重心の落ちる速度が負になるような運動の場合で，こちらの方が実際のブランコの漕ぎ方に直感的に一致しているような気がする。

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
d = Pi/2; e = 0.05; w0 = Pi; w = 2 w0;
Table[sol[i] = NDSolve[{x''[t] +
       2 e w Sin[w t + d*i] /(1 - e Cos[w  t + d*i]) x'[t] +
       w0^2/(1 - e Cos[w t + d*i]) x[t] == 0, x[0] == 0.1,
     x'[0] == 0}, x, {t, 0, 150}], {i, 0, 3}];
f[i_, t_] := x[t] /. sol[i][[1]]
Plot[Evaluate@Table[f[i, t], {i, 0, 3}], {t, 0, 10},
PlotRange -> {-0.4, 0.4}, PlotStyle -> {Red, Gray, Blue, Black}]

#

# 振幅をPi/4に強調した重心の軌跡

g1=ParametricPlot[ Evaluate@Table[{(1 - e Cos[w  t + d0*i]) Cos[w0 t], -(1 - e Cos[w  t + d0*i]) Sin[w0 t]}, {i, 0, 3}], {t, 0, 2 Pi/w}, PlotStyle -> {Red, Gray, Blue, Black},
  PlotRange -> {{-0.8, 0.8}, {-1.2, 0}}]
g2 = Plot[-Abs[x], {x, -1/Sqrt[2], 1/Sqrt[2]},   PlotRange -> {{-0.8, 0.8}, {-1.2, 0}}]
Show[g1,g2]

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

 図　パラメタ励振（支点と重心の距離は $\ell_0(1-\epsilon \cos (2 \omega_0 t + d)$ ）

図　パラメタ励振（位相d = 0赤，Pi/2灰，Pi青，3Pi/2黒）

パラメータ励振（２）

ブランコのパラメータ励振単振子モデルをMathematicaで解いてみる。
$\epsilon = a / \ell_0$は，重心の上下振幅$a$ともとの振り子の長さ$\ell_0$の比であり，$\omega_0= \sqrt{g/\ell_0}$は重心が動かない場合の振り子の固有角振動数だ。
\begin{equation*}
\ddot{\phi} + \epsilon \sin \omega t \ \dot{\phi}  + \omega_0^2 (1+\epsilon \cos\omega t) \phi = 0
\end{equation*}
重心の上下振動の振動数を振り子の振動数の２倍にしたときに励振が起こる。
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
e = 0.05; w0 = 3; w = 2 w0;
sol = NDSolve[{x''[t] + 2 e w Sin[w t] x'[t] 
    + w0^2 (1 + e Cos[w t]) x[t] == 0,
    x[0] == 0.1, x'[0] == 0}, x, {t, 0, 30}];
f[t_] := x[t] /. sol[[1]]
Plot[f[t], {t, 0, 20}]
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

図　パラメータ励振の数値計算例

パラメータ励振（１）

物理Ⅲの授業は来週中間テストでいよいよ後半戦に向う。テーマは振動と波動で，小形正男さんの裳華房の振動・波動をテキストにしている。第５章の減衰振動と強制振動のところだ。小学校５年生の振り子の単元でもブランコを説明に遣うが，ここでも強制振動の説明にブランコをイメージしてもらう。ところが他人に押してもらわない自分で漕ぐブランコは，パラメータ励振で説明しなければならない。

ブランコのモデルとしてひもの長さが与えられた時間の関数として変化する単振り子を考える。運動を記述する変数として，鉛直方向からの振り子のなす角度を $\phi(t)$ とする。支点のまわりの角運動量は，$L=m\ell ^2 \dot{\phi}$ であり，支点のまわりの重力のモーメントは $N=-mg\sin \phi$  から，運動方程式$\dfrac{d}{dt}L=N$ は，$m \dfrac{d}{dt}(\ell^2 \phi) = -m g \sin \phi$ となり，微小振動を仮定して整理すると，$\ddot{\phi} + 2\dfrac{\dot{\ell}}{\ell}\dot{\phi} + \dfrac{g}{m} \phi = 0$ を得る。

例えば，$\ell(t) = \ell_0 - a \cos \omega t$ に対して，この方程式が簡単に解ければ，初等的な力学の教科書に例題としてどんどん載せればいいようなものだけれど，そうは問屋が卸さないのだった。ちなみにこの場合の近似運動方程式は，$\epsilon = a / \ell_0$，$\omega_0
= \sqrt{g/\ell}$として，次のようになる。
\begin{equation*}
\ddot{\phi} + \epsilon \sin \omega t \ \dot{\phi}  + \omega_0^2 (1+\epsilon \cos\omega t) \phi = 0
\end{equation*}

眉村卓

11月3日に眉村卓が亡くなった。中学生のころからＳＦを意識して読み始めるようになった。その当時の日本のＳＦ作家では，小松左京，筒井康隆，星新一に次いで読んでいたのが眉村卓かもしれない。阪大経済学部で柔道部というところにもちょっと引っ掛かったかもしれない。眉村の作品で書棚に並んでいるものを調べてみた。燃える傾斜（ハヤカワ文庫），準Ｂ級市民（ハヤカワＳＦシリーズ），万国博がやってくる（ハヤカワＳＦシリーズ），幻影の構成（世界ＳＦ文学全集），EXPO87（日ハヤカワ文庫），滅びざるもの（徳間文庫），消滅の光輪（ハヤカワ文庫），司政官全短編（創元文庫）などであり，1960年から70年代にかけての作品に集中している。

眉村卓といえばインサイダー文学論だ。コリン・ウィルソンのアウトサイダーは大学に入ってから非常に興味深く読んだ。ただ，この実存主義的なアウトサイダーとインサイダー文学論のインサイダーはベクトルの向きが正反対というわけではないかもしれない。それでもＳＦ読者としては，アウトサイダーという立場がなぜかしっくり来ることが多い。そう，自分が普通の人類と違っていたらどうだろうというあの感覚だ。

インサイダー文学論は，ＳＦの対象として，そのようなアウトサイダー的傾向に親和的な層だけではなく，社会組織の中で歯車として生きているインサイダー層にフォーカスするような作品が必要であるという眉村自身の創作の立場を説明するものだ。それは初期作品にも見られたし，中期から後期にかけての司政官シリーズで端的に表現されていたのだと思う。

［１］私の失敗「眉村卓さん」（産経新聞文化部）
［２］眉村卓『司政官　全短編』あとがき（2008年1月）
［３］眉村卓『消滅の光輪』あとがき（2008年7月）
［４］眉村卓「霧を行く」（橄欖追放，東郷雄二）

スピン

日本の主な国語辞典のスピンの項目には，まだスピン（パブリック・リレーションズ）としての語意が載っていない。英語版のWikipediaでは，Spin（Propaganda）となっていて，日本でよく用いられるスピン報道とは若干ニュアンスが違うような気もする。

マスコミは基本的に読者の目を引いてメディアへのアクセスを増やし，広告の獲得や購買につながるニュースを仕立てている。その事情を熟知していて，かつ，自由に情報発信を制御できる政府は，あたかも自然現象のようにニュースを作ることができるかもしれない。

［１］「桜を見る会」と芸能報道から考える「結果スピン」の効能（荻上チキ）

公共財としての教育ビッグデータ（２）

例えば，JAPAN e-Portfolio は関連大学が構成員となる一般社団法人教育情報管理機構が運用主体となっている。その会員についてのページを見ると，正会員は大学である。国立大学が６法人，公立大学が３法人，私立大学が１８法人となっている。賛助会員には教育産業の企業が名を連ねている。ベネッセが含まれる特別賛助会員が４社，指定賛助会員が４社，賛助会員が１社である。年会費300万円特別賛助会員の説明は以下の通り。

当該会員が運営する学習支援システム事業，ポートフォリオ事業，SNS事業，データベース事業等これらに類する事業において取得したデータ又は本機構が運営する高大接続ポータルサイト「JAPAN e-Portfolio」(以下、「JeP」という。)が所有するデータ等，JePに連携し，蓄積した情報を活用して事業を行う会員が支払う会費
＊『「JAPAN e-Portfolio」(以下、「JeP」という。)が所有するデータ等」』とは「JAPAN e-Portfolio」の入力項目などシステム本体の情報を指します。生徒が「JAPAN e-Portforio」に入力し，蓄積された個人情報は，一切共有されることはありません。

ベネッセの個人情報漏れが起こったときには，これでもうおしまいだなと思ったのだが，あっという間に，復活を果たしただけでなく，ありとあらゆる教育公共事業案件に首を突っ込み，この優位性を背景に各学校への攻勢を強めている。これまでは，全国学力テストだったが，これに加えて，大学入試共通テストや，高校生のポートフォリオなど膨大なデータを手中に収めようとしている。このままでいけば文部科学省は民営化されて，ベネッセに取って代わられる日も近い。

公共財としての教育ビッグデータ（１）

日経朝刊で，公正取引委員会の杉本和行委員長が，「医療や金融分野などのビッグデータは多くの事業者が利用可能な「公共財」になりうるとの考え方を示した」。ここには教育分野が含まれていないが，まさに公共財として位置づけて，広く公開と流通を図ることが必要な部分は含まれているはずだ。もちろん，個人情報の機微に関る部分も多いため，その取り扱いにも注意は必要だが，それは医療や金融でも同じことだ。

数年前から清和会・政府の進めている教育政策の特徴は，一つには右翼的イデオロギー貫徹による教育支配であり，もう一つが，公共財としての教育を売り払う利権支配である。そして，その矛盾が噴出したのが，森友・加計・大学入試問題だ。巷では，モリ・カケ・サクラとして安倍の不正をあげつらっているが，本当にこわいのは，もてあそばれている日本の教育システムの破壊ではないのか。

岡山という地域をベースにしたベネッセと加計学園グループの深い結びつきについては，いろいろとネット上に情報が流れているが，そこに様々な公益法人や文部科学省からの天下りなどが絡み合い，複雑な病巣を構成している。そこに巣食う企業人や大学人や官僚崩れが，トップダウンの準公的な諮問会議や私的なワーキンググループを隠れみのに互いの利権を貪りながら，十分に吟味されておらず，当事者や専門家を排除した政策を立案している。

［１］教育情報管理機構（JAPAN e-Portfolioの実施大学の団体＋特定賛助会員など）
［２］進学基準機構（佐藤禎一が理事長のベネッセ系一般財団法人）
［３］学力評価研究機構（ベネッセグループの株式会社）
［４］福武教育文化振興財団（ベネッセグループの公益財団法人）
［５］高校生のための学びの基礎診断（文部科学省）

米山保三郎

日経朝刊の連載小説が伊集院静の道草先生だ。最初はフォローしていなかったが，漱石の話だったことに気づき時々読んでいる。夏目金之助（1867-1916）の大学予備門（第一高等中学校）での友人として金沢出身の米山保三郎（1869-1897）がでてきた。金沢といえば宇ノ気町出身の西田幾多郎（1870-1945）も同時代の空気を吸っている。

米山保三郎については，Yahoo！知恵袋に北國・富山新聞の情報として，生没年と生誕地が野町２丁目であること，父が加賀藩の算用侍であったことなどが記されていた。東京帝大の哲学科で学び，かなりの秀才だったようだ。「吾輩は猫である」では米山をモデルとした天然居士＝曾呂崎が登場している。米山は実際に天然居士という号を持っている。米山の助言によって金之助＝漱石が建築ではなく文学を目指すことになった。

まだ子供のとき，財産がなかったので，一人で食わなければならないという事は知っていました。忙がしくなく時間づくめでなくて飯が食えるという事について非常に考えました。しかし立派な技術を持ってさえいれば，変人でも頑固でも人が頼むだろうと思いました。佐々木東洋という医者があります。この医者が大へんな変人で，患者をまるで玩具か人形のように扱う，愛嬌のない人です。それではやらないかといえば不思議なほどはやって，門前市をなす有様です。あんな無愛想な人があれだけはやるのはやはり技術があるからだと思いました。それだから建築家になったら，私も門前市をなすだろうと思いました。丁度ちょうどそれは高等学校時分の事で，親友に米山保三郎という人があって，この人は夭折しましたが，この人が私に説諭しました。セント・ポールズのような家は我国にははやらない。下らない家を建てるより文学者になれといいました。当人が文学者になれといったのはよほどの自信があったからでしょう。私はそれで建築家になる事をふっつり思い止まりました。私の考は金をとって，門前市をなして，頑固で，変人で，というのでしたけれども，米山は私よりは大変えらいような気がした。二人くらべると私が如何にも小ぽけなように思われたので，今までの考をやめてしまったのです。そして文学者になりました。（夏目漱石「無題」，大正三年一月十七日東京高等工業学校において）

［１］落第（夏目漱石，青空文庫）
［２］無題（夏目漱石，青空文庫）
［３］処女作追懐談（夏目漱石，青空文庫）
［４］明治二十四，五年頃の東京文科大学選科（西田幾多郎，青空文庫）
［５］漱石と自分（狩野亨吉，青空文庫）
［６］吾輩は猫である（夏目漱石，青空文庫）
［７］漱石の悼む大怪物（かわうそ亭）
［８］漱石とメリメ : 『吾輩は猫である』における『カルメン』の水浴（高木雅恵）
［９］漱石覚え書補篇（柴田宵曲）

（写真：夏目漱石と米山保三郎，熊本市広報記事「漱石とくまもと」より引用）

東大寺七重塔

奈良東大寺の創建当時，東西に七重塔があったようだ。最近の東塔の発掘調査でも東塔を囲む回廊の東門の跡が見つかっている。塔の高さは70mとも100mともいわれており，東大寺はその再建を目標としているとのこと。ほんとですか。

日本で一番高い木造の塔は，54.8mの東寺の五重の塔だ。世界で一番高い木造建築物は，58.5mのカナダのブリティッシュコロンビア大学の学生寮だ。70mとか100mは可能なのだろうか？建築基準法上は燃えなければ可能みたい。

その東大寺の七重塔は，1970年の大阪万博の際に古河グループの古河パビリオンとして復元されたことがある。ただしこれは木造ではなくプレハブ構造で１年で作ったようだ。その際に作られた七重塔相輪のレプリカは，東大寺大仏殿の横に設置されている。

金沢から大阪まで，雷鳥の臨時特急に乗って万博に行ったのは高校２年の夏休みだった。が，古河パビリオンにはあまり印象がない。岡本太郎のデザインで，小松左京のコンセプトによる内部展示を楽しみにしていただけにちょっと期待外れだった70mの太陽の塔には登ったけれど。

哲学 v.s. 物理学

勁草書房の〈現在〉という謎：時間の空間化批判（森田邦久編著）における哲学者と物理学者の論争は，場外にまで発展した。

第４章　客観的現在と心身相関の同時性［青山拓央］
　コメント：哲学者に考えてもらいたいこと［谷村省吾］
　リプライ：まず問いの共有を［青山拓央］

第５章　時間に「始まり」はあるか――哲学的探究［森田邦久］
　コメント：物理学者が哲学者の時間論を読むとこうなる［谷村省吾］
　リプライ：哲学者も物理学を無視しない――形而上学と物理学の関係性［森田邦久］

これを受けて，名古屋大学の谷村省吾は，『〈現在〉という謎』をめぐる議論 のページで，『一物理学者が観た哲学』　2019年11月6日公開　PDFファイル　総ページ数 110
（書籍『〈現在〉という謎―時間の空間化批判』への補足ノート 増補版）を書いた。

かなり刺激的に哲学者（青山拓央）に喧嘩を売っている。谷村の論点はほぼ理解できる。しかし，谷村と同様に自分にも哲学者の方の言い分はよくわからない。これまでに何冊も科学哲学の本を読んでいるが，哲学者の基本用語がほとんど〈それ〉として理解できないままだ。どうしたものだろうか。

［１］物理学者と哲学者は「時間」を語ってどうすれ違うのか。（重ね描き日記）
［２］哲学者は物理学者の本気の拳をどう受け止めるか…谷村省吾「一物理学者が観た哲学」を読んで（丸山隆一 note）

大学改革の迷走：佐藤郁哉

ちくま新書1451の「大学改革の迷走」を読了した。478ページでかなり分厚い新書本だが２日で読むことができた。ほとんど，自分が体験してきたことと並行しているので，容易に頭に入ったからだ。著者は東北大学の大学院で心理学を専攻した後にシカゴ大学で社会学の Ph. D. を取り，一橋大学の教授から現在は同志社大学の商学部の教授である。社会調査方法論や組織社会学が専門であり，大学改革の真っただ中でのご自分の経験をも踏まえつつ，社会学者らしい視点が随所に見られる。

各章のタイトルではなくサブタイトルとキーワードを並べたほうがその内容が分かる。
第１章　—和製シラバスの呪縛（シラバス，画一化）
第２章　—和製マネジメント・サイクルの幻想（ＰＤＣＡ，経営ごっこ）
第３章　—残念な破滅的誤解から創造的誤解へ（選択と集中，ＫＰＩ）
第４章　—実質化と形骸化のミスマネジメント・サイクル（慢性改革病，改革ごっこ）
第５章　—大学院拡充政策の破綻と「無責任の体系」（集団無責任体制）
第６章　—改革劇のドラマツルギー（作劇術）を越えて（道徳劇）
第７章　—「大人の事情」を越えて（エビデンス，ＥＢＰＭ，ＧＩＧＯ）

自分も大学の中で，この壮大なインチキを支えて自他を半ば騙しながら改革ごっこの推進側でもがいていた。したがって，あたかも第三者のように大学改革の迷走の犯人・悪者探しをすることはできないはずなのだが，あらゆる場面で道徳劇への参加の誘惑は強い。

５万人の学生調査で問題を抱えたまま，60万人を対象とする調査を繰り返すという話があったが，41万人の全国学習調査のことだろうか。あるいは，大学３年生60万人を対象とした授業評価調査のことだろうか。

桜を見る会

桜を見る会は，例年４月に新宿御苑で行われる内閣総理大臣主催の公的行事である。1700万円の予算を大幅に超過する支出が続いたことが問題視されたかと思ったら，逆に実態を追認して来年度予算を３倍の5700万円にする始末。何年か前に，安倍を批判していた爆笑問題がよばれて安倍と並んでいたときにいやーな感じはしたよね。11月8日の参議院予算委員会での日本共産党の田村智子議員の質疑が注目されている。最初，大手マスコミはほとんど沈黙したままだったが，少しだけ流れが変わったか。これも新しい燃料がなければすぐに沙汰やみになるのだろう。

P. S. 2019.11.13 橋下−NHK−菅 ラインで，燃料を除去（来年度中止）して火消しに走った模様。安倍晋三後援会のホテルニューオータニ前夜祭＋富士急行バスツアーの問題も政治資金報告書の問題も何も解決していないはずだけど・・・

P. P. S. 2019.11.15 NHK社会部も結構攻勢をかけていたが，安倍が説明するといっているので（実際，夕刻に記者のぶらさがりで実施），つっぱねるための言い訳の作文ができたのだろう。追求側は証拠が足りず今回も詰めが甘かったということか・・・（今場所の遠藤の相撲みたい）

［１］「桜を見る会」を全部見る！テレ東NEWS：テレビ東京

英語民間試験の国会質疑から

京都工芸繊維大学の羽藤（はとう）由美先生は，京都工芸繊維大学の「大学入試への英語スピーキングテスト導入プロジェクト」のリーダーである。京都工芸繊維大学では，入試改革の一環として，英語スピーキングテスト導入に向けた学際的な研究グループを立ち上げており，国内でも最も先進的な取り組みをしている外国語教育の専門家である。先日の衆議院の文部科学委員会の参考人意見陳述に感銘を受けたので，ここに文字起した。

羽藤由美　参考人（京都工芸繊維大学教授）
意見陳述　11/5　衆議院・文部科学委員会
高大接続改革［英語民間試験］9:08-10:32

英語四技能の育成は重要です。しかし教育のインフラともいえる共通テストを民間に委ねるということは，国にとって重大な判断です。それによって得られるものと失なうものの大きさの比較検討が一切なされていません。
財や名を成した素人が，どこか高いところに集まって，個人的な経験や感想を言いあい，その中で決めた現実味のない教育政策が，推進に無批判に協力するごく少数の研究者や教員を利用する形で，そのまま現場に降りてきます。この現状こそどうぞ改善して下さい。
この国には，英語教育，言語テスト，テスト理論など能力の高い研究者がたくさんいます。教育現場にも地味に研鑽を積み，着実な成果をあげている先生方がいらっしゃいます。どうかその人たちの専門知を結集して，入試に頼らない教育の在り方も含めて，実現可能な最適解を探す努力をして下さい。
今回の検討会議が，そういう会議となるようなご配慮をお願いいたします。

［１］東京大学入学者選抜方法検討ワーキング・グループ答申（2018.07.12）

仮名手本忠臣蔵（八段目より十一段目まで）

先の金曜日，国立文楽劇場の仮名手本忠臣蔵（八段目から十一段目）をみてきた。大ホール前の台風19号などの水害募金チャリティで，吉田玉男さんが珍しく姫の人形を持って募金者との記念写真に応じていたので，記念撮影してもらった。

八段目は，道行旅路の嫁入りだ。どうも道行きは苦手なので途中半分は寝ていたかもしれない。加古川本蔵の妻の戸無瀬と娘の小浪が許嫁の大星力弥に会うため京都の山科まで向う冬の場面。浄瑠璃の出発点は薩埵峠（現在の由比町と静岡市と境で富士山がよく見える）。なにやらなまめかしいセリフがつづくなあと思っていたら，どんどん進んで行って，背景も琵琶湖にかわり，そろそろ山科につくところで終了。道行のトップを務めていた竹本津駒太夫は来年からは六代目竹本錣太夫になる。竹澤宗助の三味線の裏拍子がうまいなあと思った。人形は戸無瀬が吉田和生，小浪が吉田一輔。

九段目は，雪転しの段。雪達磨になぞらえながら祇園一力茶屋から帰ってきた大星由良之助が息子の力弥と禅問答をする。そして，いよいよ山科閑居の段。これは松竹座の歌舞伎も含めて何回かみている。戸名瀬と小浪が山科の大星家に到着し，大星由良之助の妻のお石が二人に冷たく，ピリピリとしたやり取りが延々と続くのであった。竹本千歳太夫はあいかわらずうなっていたけれど，豊澤富助はセーブしている感じ。確かにご無用の出所がふめいだった。豊竹藤太夫は，結構重い球を投げていて良かったのではないか。三味線の鶴澤藤蔵がうなりすぎるのをなんとかしてほしいという家人の説に同感である。歌舞伎で見たときに，雪が積もった竹で戸を開ける仕掛けの場面が印象深かったので，今か今かと待っていたら，この段の最後にやっとでてきた。しかも今回はちょっといまいちの戸の外れ方で，期待の方が外れてしまった。加古川本蔵は重要な役なのに，大星由良之助の仕掛けの披露に紛れていつの間にか死んでしまいかわいそうである。もう少しスポットを当てればよかったのに＞並木宗輔？。で，はっきり天河屋に続くことを宣言してこの段は幕を閉じる。

十段目の天河屋の段が，今回のスペシャル趣向になっている。住太夫の相三味線だった野澤錦糸師匠が復曲して，102年振りに原作通りの口・奥の形で上演された。前半は以前も見ているけれど，天河屋義平も怒るよね。竹本小住太夫は迫力はないかもしれないけれど，登場人物の語り分けがうまいし，心配なしで聞くことができる。今回の豊竹靖太夫はまあボチボチか。後半のストーリーがなかなか大変で，子どもを持つお母さん達にはキツイのではないか。桐竹勘介の足の黒子だけ，ちょっと高さが他とずれていたので改善してほしい。吉田簑紫郎は丁稚伊吾の役であり，人形が人形を遣うという珍しいものだった。

十一段目は花水橋引揚げより光明寺焼香の段。普通はどちらか一方だけでおしまいなのだが，今回は両者が続く。女性の人形遣いが登場したのかと思ってよく見ると，千崎弥五郎の吉田玉延君だった模様。若い人が増えるのはいいですね。それにしても太夫の層が薄いのが気になる。豊竹咲太夫師匠が７月に人間国宝の認定を受けたが，最近あまりメジャーな役をしないのは，若手育成のためだとおっしゃっている。これはこれでありかもしれない。なお，騒動のそもそもの原因である桃井若狭之助が最後に騎馬で登場して，エエカッコしているのはどうなのよという文をどこかで見たが，その通りだと思う。

興福院特別公開

奈良市法蓮町の興福院（こんぶいん）は奈良時代に創建された浄土宗知恩院派の鄙びた尼寺だ。奈良高校の裏手辺りの佐保山の丘陵地に立地し，本堂からは奈良の町並みが一望できる。普段は一般公開されていないのだが，今年20年ぶりの特別公開が行われ，予約無しで拝観することができた。秋晴れの日に参詣客がそぞろ訪れていたが，さすがにここまでは外国人観光客が到達していない模様。

国重要文化財の客殿には，襖絵や袱紗などが展示されていた。雨ざらしになりそうな渡り廊下の階段を登って，足を踏み抜きそうな橋部を慎重に踏み渡っていくと本殿に至る。本殿にあるご本尊は，国重要文化財である天平時代の阿弥陀三尊像だ。木心乾漆阿弥陀如来の印相は珍しい中品下生のもので，両脇侍は観音菩薩と勢至菩薩。庭も美しく手入れされており，くちなしの若芽が養生中だった。

  

写真：興福院の大門と本堂（2019.11.08撮影）

TEAC LP-R550

TEACのターンテーブル/カセット付きCDレコーダー LP-R550が壊れてしまった。先日，新大阪にあるTEAC修理センターに持ち込んだところ，修理不能のため，有料で新品のLP-R550USBに交換され送られてきた。家にあるレコードを聴くためのミニマムのセットである。

大学時代の下宿ではSONYのラジカセを愛用し，レコードを聴くのは自宅に帰省したときに限られていた。あ，そうではない。妹が隣棟の２Ｆに下宿していて，SONYのオーディオセットを買っており，レコードはそこで聴いていた。あれは，たぶん1978年に彼女が結婚した際に持って行ったはずだ。

1981年に自分が結婚したときには，テクニクスのシンプルなコンボを買った。父が亡くなる前にこれを金沢に持って行ったのは，レコードを聴いてもらうためだっただろうか。その後，廃棄して後継はSONYにしたのだが，こちらのほうは記憶にすらない。このころから購入したSONY製品の劣化が進んで，テレビやビデオなども軒並みにすぐ調子が悪くなった。

CDの時代になってレコードは聴かなくなったと思ったら，iPod全盛期を迎え，家族全員が持っていた。そのうちCDも買わなくなり，レンタルCDをダビングする時代になった。iPodはiPhoneに移行し，ついにstreamingの時代になってしまった。CDの全盛期は1985年から2005年までしか続かなかった。で，再びアナログレコードの復権である。かな？

 　　　

写真：テクニクスの1980年頃のモデル（Technicsより），TEAC LP-R550USB（TEACより）

子ども提灯行列

中原中也の「金沢の思い出」には次の文があった。

「今以てそれは不思議といへば，一度町内の子供が全部揃つて，忠臣蔵の真似をして練り歩いたことがある。長い列であつた。その先頭ではホラ貝を吹いてゐた。子供達ばかりでやつたことゝしては統制があり過ぎた。あれは親達も手を貸したのであらうか？　それともあゝした習慣が金沢にはあつたのであらうか？」

これ，百万石まつりのときに行われる小学生の子ども提灯太鼓行列のことかと思った。しかし，金沢百万石まつりが始まったのは，1952年の商工まつりから。あるいはルーツを辿っても，尾山神社での封国祭に合わせて，1923年（大正12年）から1945年（昭和20年）まで金沢市祭として行われてきた奉祝行事ということなので，中原中也が金沢にいた1912-1914年とはあわない。

自分たちが提灯行列に参加したのは小学校５年のときだったのではないかと思う。泉野小学校の校庭に集まり，兼六園の付近まで全員で移動する。そこには市内の小学生が終結しており，夕方になると各校下にむけて，赤い提灯を掲げながら行進した。寺町通りをある程度進むと自宅が近くなるので，流れ解散となったような気もするし，もう暗くなった夜道を小学生がばらばらと帰って行くのも危ないと思うがどうだったのか。

行列の行進の際は，金沢市民の歌（金沢市歌のほうではない）や石川県民の歌など（小学校でよく練習している）をみんなで声を張り上げてうたい気勢をあげるのであった。たぶん，何かのDNAで中也の思い出と結ばれているような気がした。

中原中也と金沢

NHKＥテレの「にほんごであそぼ」は良い。今は亡きうなりやベベン（国本武春）や文楽の竹本織太夫（豊竹咲甫太夫）や鶴澤清介，野村萬斎など豪華出演陣とよく練られた構成が相俟ってたまに見るととても良い。

昨日は，中原中也（1907-1937）の「サーカス（山羊の歌）」の一部を織太夫・清介コンビが朗読していた。朗読というのかな「ゆあーん　ゆよーん　ゆやゆよん」だけなので。

この詩「サーカス」は中也の金沢時代（1912-1914）の印象がヒントになっているといわれている。青空文庫には「金沢の思ひ出」という小文があり，中也が父に連れられて軽業の興行を見に行った神明館横の空地で，映画館神明館の弁士の息子とやりとりする場面が記されている。この神明館は犀川を寺町台地側に渡ったすぐにある今の神明宮にあった。

この神明宮で私の父の芳一が小林小児科の院長と若いときに喧嘩したと聞いたことがある。その院長の息子の小林博人君は皮膚科医になったが，小学校から高等学校までの同級生だ。野田中の２年のときは同じクラスだったので，特に仲よくしていた。マジカルミステリーツアー，ハローグッドバイ，レディマドンナの時代のビートルズに感化されたのは彼からだった。寺町で小林小児科のあとに皮膚科医として開業していたが，後に病院勤務になっていた。

中原中也の金沢でのゆかりの場所について，金沢の思い出の文を辿りながらたいへん詳しく調べているのがこちらのブログ記事「中原中也の金沢の思い出」である。若干付け加えるとするとタカジアスターゼの弟やその三男の話がある。タカジアスターゼとは高峰譲吉のことである。「高峰譲吉博士③ルーツ，高峰譲吉博士④金沢の所縁の地を歩く」によれば「高峰家の菩提寺は，金沢市寺町5丁目の臨済宗国泰寺で，高峰家のお墓には、譲吉博士の父、母、弟3人、妹3人がまつられています」，また「譲吉博士は、大正2年（1913）5月2日，亡き両親の法要のためこの国泰寺を訪れ，金沢第一中学で講演もしています。当時，寺町，大桜向いに弟の家が有ったと聞きます」とのことで，中原中也はちょうどこの時に，高峰譲吉のその弟の家の隣に住んでいたと思われる。

南円堂と北円堂

家人に勧められて，興福寺国宝特別公開【南円堂，北円堂】に行った。秋晴れの爽やかな朝，すでに30人ほどが興福寺に設置された特別展の ticket counter に並んでいた。重要文化財南円堂では，東向きの国宝木造不空羂索観音菩薩座像（康慶作）のまわりに国宝木造四天王立像（康慶作）が多聞天（北東）持国天（南東）増長天（南西）広目天（北西）と囲んでおり，本尊の左手と右手の脇に国宝木造法相六祖座像（康慶作）が玄賓，行賀，常騰と玄昉，善集，神叡が据えてあった。南円堂よりひとまわり小さな国宝北円堂には国宝木造弥勒如来座像（運慶作）の背後に国宝木造無著菩薩立像と世親菩薩立像（運慶作）が，脇の左右に法苑林菩薩と大妙相菩薩が，そして国宝木心乾漆造四天王立像が囲んでいる。

運慶の無著菩薩は前回見たときは凄い印象があったが今回は弥勒如来ともども，康慶の迫力に負けてしまった。南円堂は何周もして四天王像の特徴を観察した。両足の上げ方や火焔光背の炎の向きの対称性，両手に持つ武具と宝物，衣装や兜，体型などすべてが緻密に計算されて形造られて本尊の周りを囲んでいた。一般に，四天王は持国天（東）増長天（南）広目天（西）多聞天（北）とそれぞれの守護範囲が決まっているのだが，本尊の向きとこの四方を護持する関係で実際の方位は異なる場合もあるようだ（見張りのお坊さんに尋ねてみた応答からの推測）。

その後，再建された中金堂にも回ってみたが，こちらは着飾った新婚外国人カップルがプロカメラマンを従えて中金堂を背景に記念スチール写真の撮影会をやっているような状況で，これはこれで大変雰囲気にマッチしていたが，中に収められた江戸時代の釈迦如来座像などは今一つだった。国宝木造四天王立像や重要文化財薬王・薬上菩薩立像なども再建された中金堂の空間の雰囲気になじめず，残念な状況だった。

　　　　　　

写真：南円堂と北円堂（撮影　2019.11.05）

反対称テンソルと対称和

ある量の対称和が必要な場合がある。これを，レヴィ=チヴィタ記号$\varepsilon_{ijk}$で表現する例をいくつか考えてみた。$\varepsilon_{ijk}^2 = (1-\delta_{ij}) (1-\delta_{jk}) (1-\delta_{ki})$が成り立つような気がする。
\begin{equation}
\begin{aligned}
A_1 B_2 C_3 + A_2 B_3 C_1 + A_3 B_1 C_2 &= \dfrac{1}{2} \sum_{ijk} (\varepsilon_{ijk} + \varepsilon_{ijk}^2) A_i  B_j C_k\\
A_1 B_2 B_3 + A_2 B_3 B_1 + A_3 B_1 B_2 &= \dfrac{1}{2} \sum_{ijk} \varepsilon_{ijk}^2 A_i  B_j B_k\\
A_1 + A_2 + A_3 &= \dfrac{1}{2} \sum_{ijk}\varepsilon_{ijk}^2 A_i
\end{aligned}
\end{equation}

民間試験導入としての記述問題

大学入学共通テストへの英語の民間試験導入は，幅広い層の反対論をベースに萩生田のオウンゴールがアシストとなって，11月1日の共通ＩＤ提出開始のタイミングで延期となった。『11月１日に文部科学大臣から、令和３年度大学入学者選抜から導入予定であった英語民間試験活用のための「大学入試英語成績提供システム」の導入を見送るとの発表がありました』

しかし，文部科学省は，もう一つのベネッセ案件である国語・数学における記述試験はこのまま進めると主張している。その問題点については，「大学入試改革｢炎上｣の裏に潜むもう1つの火種英語のみならず数学と国語でも民間試験導入」で指摘されている。結局すべてが，公共財を私企業へ売り渡たす利権モデルということか。大阪維新と完全に相似形の安倍政権。

［１］高校も大学も頭を抱える「センター試験改革」あまりにも多すぎる問題点（2019.1.19）
［２］高等教育局大学振興課説明資料（2019.2.12）
［３］英語だけじゃない…大学入試改革の「国語記述式問題導入」の害悪（2019.11.01）
［４］残る選択肢は「センター試験」続行？英語の次は国語……大学入試改革で起こるドミノ倒し（2019.11.06）

Makie.jl

Juliaのパッケージ Grassmann.jl がリニューアルされたというニュースが伝わってきたので，早速調べてみたら， Makie.jl を用いたベクトル場の流線表示の図形が載っていた。Makieは日本語の蒔絵に由来してネーミングされた，GPUを用いる高機能なjulia用グラフィックスのパッケージのようだ。

さっそく，Makie.jlをPkg.add("Makie") してみたがなかなかうまくいかない。そもそも例題が実行できないのだ。あの物性理論の永井祐紀さんが，「Juliaで綺麗なプロットを作る：Makie.jlのインストールと使い方」という記事を2018年12月に書いていたので，早速試してみた。

まず，AbstractPlotting.jlとMakie.jl と GLMakie.jl をインストールせよとある。そうなのか。
サンプルはすべて，FileIO.jl を使って，save("filename.png", scene)としている。そうなのか。ちなみに，ファイルに保存せずに直接画面に出力しようとすると，No renderer could be found for output. It has the following MIME types: というエラーが出てしまう。

とりあえず，指示に従うと，サンプルファイルはほぼ再現できたが，残念ながら，minimum(　) があるものはすべて，収束しないエラーで挫折してしまった。あと，streamplotの中のvectorfieldやlinesの中のpointなどがないといわれる。何が足りないのか？

成功した例は次のようなもの。
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
using AbstractPlotting

 x = range(0, stop = 2pi, length = 80)
 f1(x) = sin.(x)
 f2(x) = exp.(-x) .* cos.(2pi*x)
 y1 = f1(x)
 y2 = f2(x)

 scene = lines(x, y1, color = :blue)
 scatter!(scene, x, y1, color = :red, markersize = 0.1)

 lines!(scene, x, y2, color = :black)
 scatter!(scene, x, y2, color = :green, marker = :utriangle, markersize = 0.1)
 save("graph.png",scene)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

図　Makieを用いた関数のプロットの例

ブルーシート

ブルーシートの記事があった。日本のトップメーカーで，岡山県にある萩原工業が，化学繊維のフラットヤーン技術を確立して，1965年にブルーシートの原型の万能シートを売り出す。最初はオレンジ色だったのが，カドミウムが含まれるという誤ったウワサを払拭するため，1970年代にブルーになった。国産シートの９割のシェアを持つが，最近は中国産が主流となっているなど，興味深い話だった。耐久性等，品質の観点からは輸入品には負けないらしい。

米国ではカラーコードによって強度と厚みが共通化されており
　青＝ライト・デューティ＝0.13-0.15 mm
　黄色/オレンジ＝ミディアム・デューティ0.18-0.20 mm
　緑＝ミディアム・デューティ＝0.23-0.25 mm
　シルバー＝ヘビーデューティ＝0.28-0.30 mm
　茶＝スーパーヘビーデューティ＝0.41mm

ということで，アシスト瓦のシルバーはこれに相当するということなのだろうか。