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

diff --git a/tests/himmelblau_tests.m b/tests/himmelblau_tests.m
new file mode 100644
index 0000000..75a145f
--- /dev/null
+++ b/tests/himmelblau_tests.m
@@ -0,0 +1,26 @@
+## Test the optimal points.
+
+unit_test_equals("himmelblau(3,2) == 0", ...
+		 0, ...
+		 himmelblau(3,2));
+
+x1 = -2.805118086952745;
+x2 = 3.131312518250573;
+msg = sprintf("himmelblau(%f, %f) == 0", x1, x2);
+unit_test_equals(msg, ...
+		 true, ...
+		 norm(himmelblau(x1, x2)) < 1e-12);
+
+x1 = -3.779310253377747;
+x2 = -3.283185991286170;
+msg = sprintf("himmelblau(%f, %f) == 0", x1, x2);
+unit_test_equals(msg, ...
+		 true, ...
+		 norm(himmelblau(x1, x2)) < 1e-12);
+
+x1 = 3.584428340330492;
+x2 = -1.848126526964404;
+msg = sprintf("himmelblau(%f, %f) == 0", x1, x2);
+unit_test_equals(msg, ...
+		 true, ...
+		 norm(himmelblau(x1, x2)) < 1e-12);