1001.cpp

来自「数值分析的c语言几种算法实现」· C++ 代码 · 共 64 行

#include <stdio.h>



void Series_Sum( double sum[] );



int main()

{

    int i;

    double x, sum[3001];



    Series_Sum( sum );



    x = 0.0;

    for (i=0; i<11; i++)

        printf("%6.2lf %16.12lf\n", x + (double)i * 0.10, sum[i]);



    return 0;

}

void Series_Sum( double sum[] )
{
    int i,j;
    double dX = 0;
    double Sum = 0.0;
    double S = 0.0;
 for(i=1400000; i>30000; i--)
 {
  S += 1/((double)(i)*(double)(i)*(double)(i+1));
  }
 for(i=0; i<11; i++,dX+=0.1)
 {
  Sum = S;
  for(j=30000; j>=1; j--)
   {
   Sum +=  1/(((double)(j)*(double)(j+1))*double(j+dX));
  }
  Sum *= (1-dX);
  sum[i]=1+Sum;

 }
 //printf("%f\n",dX);
  for(i=11; i<3001; i++,dX+=0.1)
 {
    sum[i]=((double)(1)/(double)(dX)+sum[i-10]*(double)(dX-1))/(double)(dX);

 }
}