Page 301 - MATLAB an introduction with applications

286 ——— MATLAB: An Introduction with Applications

Solution:

>> % Penalty function method
>> clear,clf
>> f=’f321p’;
>> x0=[0.4 0.5];
>> TolX=1e–5;TolFun=1e–9;alpha0=1;
>> TolX=1e–5;
>> MaxIter=100;
>> [x0_Nelder,f0_Nelder]=opt_Nelder(f,x0,TolX,TolFun,MaxIter) % Nelder
method
x0_Nelder =
1.4423 0.6540
f0_Nelder =
0.3176
>> [fc_Nelder,f0_Nelder,c0_Nelder]=f321p(x0_Nelder) % Its results
fc_Nelder =
0.3176
f0_Nelder =
0.3176
c0_Nelder =
–1.4423
–0.6540
–0.0002
–1.4613
–0.0002

>> [x0_s,f0_s]=fminsearch(f,x0) %MATLAB built-in fminsearch()
x0_s =
1.4421 0.6540
f0_s =
0.3178
>> [fc_s,f0_s,c0_s]=f321p(x0_s) % its results
fc_s =
0.3178
f0_s =
0.3178
c0_s =
–1.4421
–0.6540
–0.0009
–1.4618
–0.0001