月xを入力したら、その月の日数を出力する関数f(x)の作り方

x=1だったら31
x=2だったら28
x=3だったら31
x=4だったら30
x=5だったら31
x=6だったら30
x=7だったら31
x=8だったら31
x=9だったら30
x=10だったら31
x=11だったら30
x=12だったら31

を出力する多項式を生成するには、自由度が12個必要なので、
xのゼロ乗から11乗までの、11次多項式が必要になります。


という話を耳にしまして

どうやったらこの式が導出できるのでしょうか。


まず、以下のような式を仮定します。

このxに1～12を入れた際に、先ほどのような日数が出るように
係数Anを定めていくわけです。


などといったように、1～12月までの12本の式を、12元連立方程式にします。
これを行列で表すと、以下のようになるので

これをXA=Bと定義して
両辺左側から、Xの逆行列を掛け算すると、係数A0～A11までが求まります。

おそらく、クラメルの方法を用いてもよいでしょう。
(行ベクトルバージョンのクラメルも試してみたい)


昨日の行列式の意味はこれです。


ツェラーとは何の関係もないことがわかりました。

ついでに、gifアニメも載せておきます。(分母と分子が逆になってました)