#include <stdio.h>
#include <math.h>
#include <stdlib.h>

#include "nr.h"

#define DIM  2

double function(double x[])
{
	return fabs((x[0]-2.0)*(x[0]-2.0)*(x[0]-2.0)) + (x[1]-10.0)*(x[1]-10.0);
}

int main()
{
	int i, j, ndim, iter;
	double *p, **xi, fmin;
	double ftol = 1.0e-8;

	ndim = DIM;

	xi = matrix(ndim, ndim); // 2次元のメモリ確保,方向行列
	p = vector(ndim);		// 1次元のメモリ確保

	for (i=0;i<DIM;i++){
		p[i] = 1.0;
		for (j=0;j<DIM;j++){
			xi[i][j] = (double)rand() / (double)RAND_MAX;
		}
	}

	powell(p, xi, ndim, ftol, &iter, &fmin, function); // 入力:p,xi,ftol, 出力:p,xi,iter,fmin

	for (i=0;i<DIM;i++){
		printf("極値を満たす変数:p[%d]    %lf\n",i,p[i]);
		for (j=0;j<DIM;j++){
			printf("方向ベクトル:xi[%d][%d]   %lf  ",i,j,xi[i][j]);
		}
		printf("\n\n");
	}
	printf("反復回数:%d  極小値:%lf\n",iter,fmin);

	free_matrix(xi); //
	free_vector(p);

	return 0;
}