A = rand(10);
A_hat = perturb(A, 1e-12);
dA = A_hat - A;

unit_test_equals("dA has the correct infinity norm", ...
		 1e-12, ...
		 norm(dA, 'inf'));
