xpade.c

一些常用数值算法的 C程序源代码 对初级学数值计算很有帮助

#define N 8

void padeps(float *a, int n, int m, float *s);

main()
{
    float a0[N+1],a[N], s;

    int m,i;
      printf("           M         S           Pi\n");
      a0[1] = 4.0;
      a[1] = 4.0;
      for(i=2; i<=N; i++) {
        a0[i] = -a0[i-1] * (2.0*i-3) / (2.0*i-1.0);
        a[i] = a[i-1] + a0[i];
      }
      for(m=N/2-1; m<=N/2+1; m++) {
        padeps(a,N,m,&s);
        printf("%d,   %f,    %f\n",  m, s, 4.0 * atan(1.0) );
      }
}