X-Git-Url: http://gitweb.michael.orlitzky.com/?p=octave.git;a=blobdiff_plain;f=tests%2Ffixed_point_method_tests.m;fp=tests%2Ffixed_point_method_tests.m;h=df3b05e721053e8de3129de68d0fc9cb737e7b52;hp=0000000000000000000000000000000000000000;hb=ef3085ef3c40eab925219bad1c6dbb2c6e106727;hpb=33e4dd7a64f35fb22774d59af11f5670ab566fbc

diff --git a/tests/fixed_point_method_tests.m b/tests/fixed_point_method_tests.m
new file mode 100644
index 0000000..df3b05e
--- /dev/null
+++ b/tests/fixed_point_method_tests.m
@@ -0,0 +1,19 @@
+g = @(x) 1 + atan(x);
+expected_fp = 2.1323;
+tol = 1 / 10^10;
+x0 = 2.4;
+unit_test_equals("Homework #2 problem #5 fixed point is correct", ...
+		 expected_fp, ...
+		 fixed_point_method(g, tol, x0));
+
+
+h = 0.5;
+g1 = @(u) 1 + h*exp(-u(1)^2)/(1+u(2)^2);
+g2 = @(u) 0.5 + h*atan(u(1)^2 + u(2)^2);
+my_g = @(u) [g1(u), g2(u)];
+tol = 1 / 10^9;
+u0 = [1,1];
+expected_fp = [1.0729, 1.0821];
+unit_test_equals("Homework #3 problem #3i fixed point is correct", ...
+		 expected_fp, ...
+		 fixed_point_method(my_g, tol, u0));