#include "simplex.h"

#define sqr(x) ((x) * (x))

sum_residual(x)
double    *x;
{
     register int   i;

     x[N-1] = 0.0;

     for (i=0 ; i<np ; i++)
/*          x[N-1] += sqr(f(x, data[i]) - data[i][1]); */
       x[N-1] +=  f(x, data[i]); 
/* mod to force all params to be > 0 */
     for (i=0 ; i< M ; i++)
	{
	if (x[i] < 0.0) x[N-1]= 9.999e99;
	}
}