X-Git-Url: http://gitweb.michael.orlitzky.com/?p=octave.git;a=blobdiff_plain;f=optimization%2Ftest_functions%2Fhimmelblau.m;fp=optimization%2Ftest_functions%2Fhimmelblau.m;h=c973a866ceb71428965abed7f9ba79764c6a6942;hp=0000000000000000000000000000000000000000;hb=4cf95f960c61165130795e7b5197cfe6e2829692;hpb=336f1181f7f065a07ee58bd772a875f7f4b39247

diff --git a/optimization/test_functions/himmelblau.m b/optimization/test_functions/himmelblau.m
new file mode 100644
index 0000000..c973a86
--- /dev/null
+++ b/optimization/test_functions/himmelblau.m
@@ -0,0 +1,17 @@
+function f = himmelblau(x1, x2)
+  ##
+  ## The eponymous function defined by Himmelblau in his Applied
+  ## Nonlinear Programming, 1972.
+  ##
+  ## It has four identical local minima,
+  ##
+  ##   * f(3,2)                                    == 0
+  ##
+  ##   * f(-2.805118086952745,  3.131312518250573) == 0
+  ##
+  ##   * f(-3.779310253377747, -3.283185991286170) == 0
+  ##
+  ##   * f(3.584428340330492, -1.848126526964404)  == 0
+  ##
+  f = (x1^2 + x2 - 11)^2 + (x1 + x2^2 - 7)^2;
+end