% fibon, fibonaccimaximering av unimodal funktion
  clear
  a=0; b=144;  n2=89;  F2=Ff(n2);  n1=a+b-n2; F1=Ff(n1);
  disp('     a    n1    n2    b ')
  disp([a  n1  n2  b]) 
  while b-a>1
    if F1>F2, b=n2; n2=n1; F2=F1; n1=a+b-n2; F1=Ff(n1); 
       else   a=n1; n1=n2; F1=F2; n2=a+b-n1; F2=Ff(n2);
    end
    disp([a  n1  n2  b])    
  end
  maxF=[F1 F2]


