#include #include /**************** p.9のエイトケン加速適用の結果 *****************/ #define n 9 //textのnに対応(S^(0)_n) #define m (int)((n-2)/2) double S_n(int nn){ //級数の和S_n int k; double y=0.; //小さな値から順に加算して丸めの誤差を防ぐ for(k=nn; k>=0; k--) y += pow(-1,k)/(2*(double)k+1); return y; } /***** エイトケン加速の関数 *****/ double S_ni(int i, int j, double s[m+1][n+1]){ double sum; sum = (-s[i-1][j]*s[i-1][j]+s[i-1][j+1]*s[i-1][j-1]) /(-2*s[i-1][j]+s[i-1][j+1]+s[i-1][j-1]); return sum; } int main(void){ int i, j; //S^(i)_j(textに対応) double s[m+1][n+1]; for(j=0; j