Advertisement
Adytzu04

TO-l2

Mar 8th, 2013
160
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.65 KB | None | 0 0
  1.  clc;
  2.  clear;
  3.  xm=-10;
  4.  xM=10;
  5. % xm=x1;
  6.  %xM=x2-d/2;
  7.  s=(sqrt(5)-1)/2;
  8.  x1=xm+(1-s)*(xM-xm);
  9.  x2=xm+s*(xM-xm);
  10.  E=0.000001;
  11.  f1=f(x1);
  12.  f2=f(x2);
  13.  while abs(xM-xm)>=E
  14.  if(f1<f2)
  15.      xM=x2;
  16.      x2=x1;
  17.      x1=xm+(1-s)*(xM-xm);
  18.      f1=f(x1);
  19.      f2=f(x2);
  20.  else
  21.      if(f1>f2)
  22.          xm=x1;
  23.          x1=x2;
  24.          x2=xm+s*(xM-xm);
  25.          f1=f(x1);
  26.          f2=f(x2);
  27.      else
  28.          if(f1==f1)
  29.              xm=x1;
  30.              xM=x2;
  31.              x1=xm+(1-s)*(xM-xm);
  32.                 x2=xm+s*(xM-xm);
  33.                 f1=f(x1);
  34.                 f2=f(x2);
  35.          end
  36.      end
  37.  end
  38.  end
  39.    
  40.  xmin=(xm+xM)/2;
  41.  xmin
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement