% Himmelblau's Function # 3
% Range of initial points: -2 <= x <= 2
% Global minima: x=0.409951714917356
% f(x)=-2.267543938143762
% Coded by: Ali R. Alroomi | Last Update: 28 Jan. 2015 | www.al-roomi.org
 
clear
clc
warning off
 
xmin=-1.4;
xmax=2;
R=100000; % steps resolution
x=xmin:(xmax-xmin)/R:xmax;
 
for i=1:length(x)
    f(i)=((1-x(i)).^2)*((x(i)+1).^4)*((x(i)-2).^3)*x(i);
end
 
plot(x,f,'r','LineWidth',2);grid;set(gca,'FontSize',12);axis([-1.4 2 -4 8]);
xlabel('x','FontName','Times','FontSize',20,'FontAngle','italic');
ylabel('f(x)','FontName','Times','FontSize',20,'FontAngle','italic');