f1dim.cpp

实现共轭梯度法的实例

#include "stdafx.h"
extern int ncom;	/* defined in LINMIN */
extern double *pcom,*xicom,(*nrfunc)(double[]);

double f1dim(double x)
{
	int j;
	double f,*xt;

	xt=new double[ncom+1];
	for (j=1;j<=ncom;j++) xt[j]=pcom[j]+x*xicom[j];
	f=(*nrfunc)(xt);
	delete xt;
	return f;
}