Add the fixed point method and some tests.

[octave.git] / run-tests.m

diff --git a/run-tests.m b/run-tests.m
index 1682d18d20acbdf7938ba5d81c68780bd2df3351..aa285c32e38fb3f70747a75402e57ef6a8515ed2 100755 (executable)
--- a/run-tests.m
+++ b/run-tests.m
@@ -51,3 +51,24 @@ expected_A = [1, 0, 0, 0, 0; ...
 unit_test_equals("Homework #1 problem #1 Poisson matrix is correct", ...
                 true, ...
                 expected_A == poisson_matrix(4, 0, 1));
+
+
+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));