#include "nr.h" #define TOL 2.0e-4 int ncom; double *pcom,*xicom,(*nrfunc)(double []); void linmin(double p[], double xi[], int n, double *fret, double (*func)(double [])) { // 直線に沿った最小化プログラム int j; double xx,xmin,fx,fb,fa,bx,ax; ncom=n; pcom=vector(n); xicom=vector(n); nrfunc=func; // function へのポインタを代入 for (j=0;j