clear all; x = randn(4,7); y = randn(3,7); D1 = Dx(x); D2 = Dx(y); E = (sqrt(D1) + sqrt(D2))/2; isedm(E.*E) F = (log(1+D1) + log(1+D2))/2; isedm(exp(F)-1)