gsl_vector_set_all(step_size, 0.01);
gsl_multimin_fminimizer *minimiser =
- gsl_multimin_fminimizer_alloc(gsl_multimin_fminimizer_nmsimplex2rand, NX);
+ gsl_multimin_fminimizer_alloc(gsl_multimin_fminimizer_nmsimplex2, NX);
for (;;) {
/* NINPUT + 1 doubles: startpoint, epsilon for residual */