xzbrac.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 30 行

/* Driver for routine zbrac */

#include <stdio.h>
#include "nr.h"

static float fx(x)
float x;
{
	return bessj0(x);
}

main()
{
	int succes,i;
	float x1,x2;

	printf("%21s %23s\n","bracketing values:","function values:");
	printf("%6s %10s %21s %12s\n","x1","x2","bessj0(x1)","bessj0(x2)");
	for (i=1;i<=10;i++) {
		x1=i;
		x2=x1+1.0;
		succes=zbrac(fx,&x1,&x2);
		if (succes) {
			printf("%7.2f %10.2f %6s %12.6f %12.6f \n",
				x1,x2," ",fx(x1),fx(x2));
		}
	}
	return 0;
}