*! version 1.0.0 Fri Apr 23 13:19:04 PDT 1993
program define eqmat /* m1 m2 */
	version 3.1
	local m1 "`1'"
	local m2 "`2'"
	tempname dif sc
	mat `dif' = `m1' - `m2' 
	mat `dif' = `dif' ' * `dif'
	scalar `sc' = trace(`dif')
	capture assert `sc'<1e-14
	if _rc { 
		di "assertion is false; squared difference = " `sc'
		assert 1==2
	}
end