wolffd@0: function p = approxeq_pot(A, B, tol)
wolffd@0: 
wolffd@0: if nargin < 3, tol = 1e-3; end
wolffd@0: 
wolffd@0: p = approxeq(A.p, B.p, tol) & approxeq(A.u, B.u, tol);