// McKinnon function requiring restart to reach the minimum
//

// tau = 1; theta = 15; phi =  10;
// tau = 2; theta =  6; phi =  60;
   tau = 3; theta =  6; phi = 400;

lambdap = (1+sqrt(33))/8; lambdam = (1-sqrt(33))/8;

b2 = c1*(1+c1) + ((c0 <= 0) ? theta*phi*pow(fabs(c0),tau) : theta*pow(c0,tau));

{c0:c1} = search( b2, 1, 1, 1, 1, lambdap, lambdam, 0, 0);

print all; debug on; eval; print all;