程序：定義函數(shù)


%黃金分割法
function f=gnew(f,a,b,epsilon);%定義函數(shù)和精度
la=0.618;
x1=b-la*(b-a);f1=f(x1);%左試點
x2=a+la*(b-a);f2=f(x2);%右試點
for k=1:(log(epsilon/(b-a))/log(0.618)+2)
    fprintf(1,'   迭代次數(shù) k= %3.0f\n',k)
    if f1<=f2
        b=x2;x2=x1;f2=f1;
        x1=b-la*(b-a);f1=f(x1);
    else
        a=x1;x1=x2;f1=f2;
        x2=a+la*(b-a);f2=f(x2);
    end
    x=0.5*(b+a);
    fprintf(1,'   迭代區(qū)間-左端 a=% 3.4f\n',a)
    fprintf(1,'   試點1坐標值    x1=% 3.4f\n',x1)
    fprintf(1,'   函數(shù)值         f1 =% 3.4f\n',f1)
    fprintf(1,'   迭代區(qū)間-右端 b=% 3.4f\n',b)
    fprintf(1,'   試點2坐標值    x2=% 3.4f\n',x2)
    fprintf(1,'   函數(shù)值         f2=% 3.4f\n',f2)
   fprintf(1,'   區(qū)間中點        x=% 3.4f\n',x)
   disp ' '
end



調(diào)用函數(shù)，輸入條件，運行結(jié)果；>> format long;gnew(inline('x^2+2*x'),-3,5,0.001),format short;
   迭代次數(shù) k=   1
   迭代區(qū)間-左端 a=-3.0000
   試點1坐標值    x1=-1.1114
   函數(shù)值         f1 =-0.9876
   迭代區(qū)間-右端 b= 1.9440
   試點2坐標值    x2= 0.0560
   函數(shù)值         f2= 0.1151
   區(qū)間中點        x=-0.5280
 
   迭代次數(shù) k=   2
   迭代區(qū)間-左端 a=-3.0000
   試點1坐標值    x1=-1.8326
   函數(shù)值         f1 =-0.3068
   迭代區(qū)間-右端 b= 0.0560
   試點2坐標值    x2=-1.1114
   函數(shù)值         f2=-0.9876
   區(qū)間中點        x=-1.4720
 
   迭代次數(shù) k=   3
   迭代區(qū)間-左端 a=-1.8326
   試點1坐標值    x1=-1.1114
   函數(shù)值         f1 =-0.9876
   迭代區(qū)間-右端 b= 0.0560
   試點2坐標值    x2=-0.6654
   函數(shù)值         f2=-0.8881
   區(qū)間中點        x=-0.8883
 
   迭代次數(shù) k=   4
   迭代區(qū)間-左端 a=-1.8326
   試點1坐標值    x1=-1.3868
   函數(shù)值         f1 =-0.8504
   迭代區(qū)間-右端 b=-0.6654
   試點2坐標值    x2=-1.1114
   函數(shù)值         f2=-0.9876
   區(qū)間中點        x=-1.2490
 
   迭代次數(shù) k=   5
   迭代區(qū)間-左端 a=-1.3868
   試點1坐標值    x1=-1.1114
   函數(shù)值         f1 =-0.9876
   迭代區(qū)間-右端 b=-0.6654
   試點2坐標值    x2=-0.9410
   函數(shù)值         f2=-0.9965
   區(qū)間中點        x=-1.0261
 
   迭代次數(shù) k=   6
   迭代區(qū)間-左端 a=-1.1114
   試點1坐標值    x1=-0.9410
   函數(shù)值         f1 =-0.9965
   迭代區(qū)間-右端 b=-0.6654
   試點2坐標值    x2=-0.8358
   函數(shù)值         f2=-0.9730
   區(qū)間中點        x=-0.8884
 
   迭代次數(shù) k=   7
   迭代區(qū)間-左端 a=-1.1114
   試點1坐標值    x1=-1.0061
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.8358
   試點2坐標值    x2=-0.9410
   函數(shù)值         f2=-0.9965
   區(qū)間中點        x=-0.9736
 
   迭代次數(shù) k=   8
   迭代區(qū)間-左端 a=-1.1114
   試點1坐標值    x1=-1.0463
   函數(shù)值         f1 =-0.9979
   迭代區(qū)間-右端 b=-0.9410
   試點2坐標值    x2=-1.0061
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0262
 
   迭代次數(shù) k=   9
   迭代區(qū)間-左端 a=-1.0463
   試點1坐標值    x1=-1.0061
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9410
   試點2坐標值    x2=-0.9812
   函數(shù)值         f2=-0.9996
   區(qū)間中點        x=-0.9936
 
   迭代次數(shù) k=  10
   迭代區(qū)間-左端 a=-1.0463
   試點1坐標值    x1=-1.0214
   函數(shù)值         f1 =-0.9995
   迭代區(qū)間-右端 b=-0.9812
   試點2坐標值    x2=-1.0061
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0138
 
   迭代次數(shù) k=  11
   迭代區(qū)間-左端 a=-1.0214
   試點1坐標值    x1=-1.0061
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9812
   試點2坐標值    x2=-0.9966
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0013
 
   迭代次數(shù) k=  12
   迭代區(qū)間-左端 a=-1.0061
   試點1坐標值    x1=-0.9966
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9812
   試點2坐標值    x2=-0.9907
   函數(shù)值         f2=-0.9999
   區(qū)間中點        x=-0.9937
 
   迭代次數(shù) k=  13
   迭代區(qū)間-左端 a=-1.0061
   試點1坐標值    x1=-1.0002
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9907
   試點2坐標值    x2=-0.9966
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-0.9984
 
   迭代次數(shù) k=  14
   迭代區(qū)間-左端 a=-1.0061
   試點1坐標值    x1=-1.0025
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9966
   試點2坐標值    x2=-1.0002
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0013
 
   迭代次數(shù) k=  15
   迭代區(qū)間-左端 a=-1.0025
   試點1坐標值    x1=-1.0002
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9966
   試點2坐標值    x2=-0.9988
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-0.9995
 
   迭代次數(shù) k=  16
   迭代區(qū)間-左端 a=-1.0025
   試點1坐標值    x1=-1.0011
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9988
   試點2坐標值    x2=-1.0002
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0007
 
   迭代次數(shù) k=  17
   迭代區(qū)間-左端 a=-1.0011
   試點1坐標值    x1=-1.0002
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9988
   試點2坐標值    x2=-0.9997
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0000
 
   迭代次數(shù) k=  18
   迭代區(qū)間-左端 a=-1.0011
   試點1坐標值    x1=-1.0006
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9997
   試點2坐標值    x2=-1.0002
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0004
 
   迭代次數(shù) k=  19
   迭代區(qū)間-左端 a=-1.0006
   試點1坐標值    x1=-1.0002
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9997
   試點2坐標值    x2=-1.0000
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0001
 
   迭代次數(shù) k=  20
   迭代區(qū)間-左端 a=-1.0002
   試點1坐標值    x1=-1.0000
   函數(shù)值         f1 =-1.0000
   迭代區(qū)間-右端 b=-0.9997
   試點2坐標值    x2=-0.9999
   函數(shù)值         f2=-1.0000
   區(qū)間中點        x=-1.0000
 

ans =

     Inline function:
     ans(x) = x^2+2*x