--- /dev/null
+## Test the optimal point.
+
+unit_test_equals("powell1([0;0;0;0]) == 0", ...
+ 0, ...
+ powell1([0;0;0;0]));
+
+## It should fail with the wrong number of coordinates.
+f = powell1([1;2;3]);
+unit_test_equals("powell1 fails with too few coordinates", ...
+ true, ...
+ isna(f));
+
+f = powell1([1;2;3;4;5]);
+unit_test_equals("powell1 fails with too many coordinates", ...
+ true, ...
+ isna(f));