xsobseq.c

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

/* Driver for routine sobseq */

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

main()
{
	int i,n1=(-1),n2=3;
	float *x;

	x=vector(1,n2);
	sobseq(&n1,x);
	for (i=1;i<=32;i++) {
		sobseq(&n2,x);
		printf(" %10.5f %10.5f %10.5f %5d\n",x[1],x[2],x[3],i);
	}
	free_vector(x,1,n2);
	return 0;
}