#include "simplex.h"
#include "rxn.h"

extern void lsoda();
extern void* derivs();

double f(x, d)
double    *x;
double   *d;
{
double tmp, obj;
int i, m;
/* t = 0.0; */

for (i=0; i<M;i++) {m = plist[i]; R[m].k = x[i];}

tout = d[0];
if (tout < t) { 
	istate=1 ; t = 0.0; 
	for(i=1; i <= NS; i++)  Y[i] = DY[i];
	}
  lsoda(derivs,NS,Y,&t,tout,itol,rtol,atol,itask,&istate,iopt,jt,
	iwork1, iwork2, iwork5, iwork6, iwork7, iwork8, iwork9,
	rwork1,rwork5,rwork6,rwork7);

if (istate < 0) bail("in f() lsoda failed");
obj = 0.0;
for (i=1; i <= NSF; i++){
  m = splst[i]; tmp = Y[m] - d[i];
  obj += tmp*tmp;
  /* obj = Y[m]; */
}
return (obj);
}