dhrt0.c

数值分析算法描述与习题解答

  #include "math.h"
  #include "stdio.h"
  #include "dhrt.c"
  main()
  { int i,n;
    int m=6;
    double x[6],dhrtf();
    n=dhrt(-2.0,5.0,0.2,0.000001,x,m,dhrtf);
    printf("M=%d\n",n);
    for (i=0; i<=n-1; i++)
      printf("x(%d)=%13.7e\n",i,x[i]);
    printf("\n");
  }

  double dhrtf(x)
  double x;
  { double z;
    z=(((((x-5.0)*x+3.0)*x+1.0)*x-7.0)*x+7.0)*x-20.0;
    return(z);
  }