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

#define NTRIAL 5  // 大きくするな！
#define TOLX 1.0e-6
#define N 4
#define TOLF 1.0e-6

int main()
{	
	int i,j,k,kk;
	double xx,*x,*fvec,**fjac;

	fjac=matrix(N,N);
	fvec=vector(N);
	x=vector(N);
	for (kk=1;kk<=2;kk++) {
		for (k=1;k<=3;k++) {
			xx=0.2001*k*(2*kk-3);
			printf("Starting vector number %2d\n",k);
			for (i=0;i<4;i++) {
				x[i]=xx+0.2*i;
				printf("%7s%1d%s %5.2f\n","x[",i,"] = ",x[i]);
			}
			printf("\n");
			for (j=1;j<=NTRIAL;j++) {
				mnewt(1,x,N,TOLX,TOLF);
				usrfun(x,N,fvec,fjac);
				printf("%5s %13s %13s\n","i","x[i]","f");
				for (i=0;i<N;i++)
					printf("%5d %14.6f %15.6f\n",i,x[i],fvec[i]);
				printf("\npress RETURN to continue...\n");
	//			(void) getchar();
			}
		}
	}
	free_vector(x);
	free_vector(fvec);
	free_matrix(fjac);

	return 0;
}