参考: パソコンで見る天体の動き, pp.131-142
既知のいくつかのポイントを通過する未知の曲線が偶数次多項式であると假定してその既知のポイント群の外側に存在する未知のポイントを推定する。聯立方程式を解いて求める。
この数列を既知のポイントとする。
この 5 元聯立方程式を解く。
x:={((1)/(2)),1,2,3,4} © 既知のポイントの数列を定義する。 y:={((277)/(64)),−8,−167,−472,277} © 既知のポイントの数列を定義する。 p(n):=polyEval({a4,a3,a2,a1,a0},x[n]^(2)) © 偶数次多項式を定義する。 funcs:=expand(seqn(y[n]=p(n),5)) © 聯立方程式を定義する。expand() は見栄えをよくするためだけであってここでは意味はない。 linSolve(funcs,a4,a3,a2,a1,a0) © 聯立方程式を解く。 polyEval(ans,x^(2)) © 元の y に戻るかどうか確かめる。
解は () である。
すなわち補外値は である。
――――――――――――――――――――――――――――――――――――――――――――
すなおにこうしたほうがわかりやすい。
linSolve(system(y[n]=a4*x[n]^(8)+a3*x[n]^(6)+a2*x[n]^(4)+a1*x[n]^(2)+a0|n=1,y[n]=a4*x[n]^(8)+a3*x[n]^(6)+a2*x[n]^(4)+a1*x[n]^(2)+a0|n=2,y[n]=a4*x[n]^(8)+a3*x[n]^(6)+a2*x[n]^(4)+a1*x[n]^(2)+a0|n=3,y[n]=a4*x[n]^(8)+a3*x[n]^(6)+a2*x[n]^(4)+a1*x[n]^(2)+a0|n=4,y[n]=a4*x[n]^(8)+a3*x[n]^(6)+a2*x[n]^(4)+a1*x[n]^(2)+a0|n=5),a4,a3,a2,a1,a0)