#include <stdio.h>
#include <math.h>
/* golden.c : 1999 3 27 by T.Oguni */
#define N 30
double f(double x){ return (x-1.0)*(x-1.0)*(x-1.0)*(x-1.0); }
/* function */
golden(a,b,tol,result,itr)
  double a,b,tol,result[N];
  int itr[2];
{
  double beta, gamma, x0,x1,x2,x3, d, f0,f1,f2,f3;
  int it;
  beta=0.5*(sqrt(5.0)-1.0);  gamma=1.0-beta;
  x0=a; x3=b;
  if (a>b) { x0=b; x3=a; }
  d=x3-x0;
  x1=x0+gamma*d; x2=x3-gamma*d;
  f0=f(x0); f1=f(x1); f2=f(x2); f3=f(x3);
  it=0;
  while (it<itr[0]){
    it=it+1; d=beta*d;
    if (f2<f1){
      x0=x1; x1=x2; x2=x3-gamma*d;
      f0=f1; f1=f2; f2=f(x2);}
    else {
      x3=x2; x2=x1; x1=x0+gamma*d;
      f3=f2; f2=f1; f1=f(x1); }
    result[it]=0.5*(x1+x2);
    if (fabs(x2-x1)<tol) break;
  }
  itr[1]=it;
  return;
}
/* main */
main(){
  double a=-2.0, b=5.0, tol=0.001, result[N];
  int itr[2]={N,1}, it;
  golden(a,b,tol,result,itr);
  it=itr[1];
  printf("%d %f\n", it, result[it]);
}