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