最小自乗法

右辺に出てくる合計を
Σ_i y_i 　 を t[0]
Σ_i y_i x_i¹ を t[1]
Σ_i y_i x_i² を t[2]
　：
Σ_i y_i x_i^N-1 を t[N-1]

のようにあらかじめ求めておいて

A_{1 N+1} = t[0]
A_{2 N+1} = t[1]
A_{3 N+1} = t[2]

　：
A_{N N+1} = t[N-1] にしたいとき

for( i=1; i<=N; i++){ for( j=1; j<=N; j++){ 　a[i][j] = s[i+j-2]; }　　　 　a[i][N+1] = t[i-1]; } 戻る