int i;
NP = atoi(argv[1]);
+ epsilon = atof(argv[2]);
gsl_rng *rng = gsl_rng_alloc(gsl_rng_ranlxd2);
/* NINPUT + 1 doubles: startpoint, epsilon for residual */
for (i=0; i<NINPUT; i++)
INPUT[i] = scan1double();
- epsilon = scan1double();
gsl_rng_set(rng,0);