function test_hnf (A)
  quad (@nested1, 0, 1)
  function z = nested1 (x)
    z = A * x;
  end
  quad2d (@nested2, 0, 1, 0, 1)
  function z = nested2 (x,y)
    z = A * x .* y;
  end
end
