xchint.c

Numerical Recipes Software 提供的算法子程序集

/* Driver for routine chint */

#include <stdio.h>
#include <math.h>
#define NRANSI
#include "nr.h"

#define NVAL 40
#define PIO2 1.5707963

float func(float x)
{
	return x*x*(x*x-2.0)*sin(x);
}

float fint(float x)
{
	return 4.0*x*(x*x-7.0)*sin(x)-(x*x*(x*x-14.0)+28.0)*cos(x);
}

int main(void)
{
	int i,mval;
	float a=(-PIO2),b=PIO2,x;
	float c[NVAL],cint[NVAL];

	chebft(a,b,c,NVAL,func);
	/* test integral */
	for (;;) {
		printf("\nHow many terms in Chebyshev evaluation?\n");
		printf("Enter n between 6 and %2d. (n=0 to end).\n",NVAL);
		scanf("%d",&mval);
		if ((mval <= 0) || (mval > NVAL)) break;
		chint(a,b,c,cint,mval);
		printf("\n%9s %14s %16s\n","x","actual","Cheby. integ.");
		for (i = -8;i<=8;i++) {
			x=i*PIO2/10.0;
			printf("%12.6f %12.6f %12.6f\n",
				x,fint(x)-fint(-PIO2),chebev(a,b,cint,mval,x));
		}
	}
	return 0;
}
#undef NRANSI
/* (C) Copr. 1986-92 Numerical Recipes Software . */