void test() {
int i,k;
double x,sum[2001],temp;
for (i=0; i<2001;i++)
{
sum[i] = 0;
temp = 0;
for ( k=1; k<=1100; k++)
{
temp += 1/(k*(k+1)*(k+2)*(k+i*0.001));
}
sum[i] = (1-i*0.001) * (2-i*0.001) * temp + 0.25*(1-i*0.001)+1;
}
x = 0.0;
// for (i=0; i<2001; i++)
// {
// printf("%6.3f %13.12f\n", x + (double)i * 0.001,sum[i]);
// }
printf("%6.3f %13.12f\n ...\n", 0.000, sum[0]);
printf("%6.3f %13.12f\n ...\n", 1.000, sum[1000]);
printf("%6.3f %13.12f\n ...\n", 2.000, sum[2000]);
return 0;
}
int main() {
int i,k;
double x,sum[2001],SUM;
for (i=0; i<2001;i++)
{
sum[i] = 0;
for ( k=1; k<=1100; k++)
{
sum[i] += 1/(k*(k+i*0.001));
}
}
x = 0.0;
// for (i=0; i<3; i++)
// {
//printf("%6.3f %13.12f\n", x + (double)i * 0.001,sum[i]);
// }
printf("%6.3f %13.12f\n ...\n", 0.000, sum[0]);
printf("%6.3f %13.12f\n ...\n", 1.000, sum[1000]);
printf("%6.3f %13.12f\n ...\n", 2.000, sum[2000]);
}
for (i = 0; i < 300; i++)
{
index = i * 10, bias = 0;
Sum = 0;
for (x = 0.1; x <= 0.9; x += 0.1)
{
bias++;
Sum = 0;
for (k = 1; k <= 1500; k++)
{
Sum += 1.0 / (k*(k + i + x)*(k + i)*(k + i + 1));
}
sum[index + bias] = (Sum*(1 - x) + (sum[index] - sum[index + 10])) * (-x) + sum[index];
}
}`