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

#include "nrutil.h"

#define VarNum		10	// 変数の数     ndata

int main()
{
	double *x,*y,*sig;  // 説明変数x, 目的変数y, データ（目的変数）とモデル標準偏差sig
	double siga,sigb;   // a,bの標準偏差
	double a,b;			// 求めたいの線形関数のパラメータ
	double chi2, q;		// chi2:χ^2(カイ)２乗評価関数，q:当てはまりの良さの確率
	int mwt = 0;
	int i;

// 動的メモリ確保
	x = vector(VarNum);
	y = vector(VarNum);
	sig = vector(VarNum);

// 説明変数を適当に入れる．
	for (i=0;i<VarNum;i++){
		x[i] = i+1;
	}
// 目的変数を適当に入れる
	y[0] = 1;
	y[1] = 3;
	y[2] = 3;
	y[3] = 6;
	y[4] = 4;
	y[5] = 5;
	y[6] = 4;
	y[7] = 7;
	y[8] = 6;
	y[9] = 7;
// 求めたい線形関数は，a+bx
	fit(x, y, VarNum, sig, mwt, &a, &b, &siga, &sigb, &chi2, &q);

	printf("%lf %lf %lf %lf %lf %lf\n",a,b,siga,sigb,chi2,q);

	free_vector(x);
	free_vector(y);
	free_vector(sig);



	return 0;
}