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'));