function eta4e = MESTIMATE(c4n,n4e,n4ed,ed4e,e4ed,x)
m = size(n4e,1); n = size(e4ed,1);
eta4e = zeros(m,1); eta4ed = zeros(n,1);
for j = 1 : m
    area = abs(det([1 1 1;c4n(n4e(j,:),:)'])/2);
    eta4e(j) = (area).^2;
    t = c4n(n4ed(ed4e(j,:),2),:)-c4n(n4ed(ed4e(j,:),1),:);
    n = [c4n(n4e(j,:),1),c4n(n4e(j,:),2)];
    signum = (-1).^(e4ed(ed4e(j,:),1)~=j);
    sigma = (signum.*x(ed4e(j,:)))'...
        *(sqrt(sum(t.^2,2)).*(sum(n,1)/3-n))/(2*area);
    eta4ed(ed4e(j,:)) = eta4ed(ed4e(j,:)) + signum.*(t*sigma');
end
eta4ed(e4ed(:,2)==0) = 0;
eta4e = sqrt(eta4e + 0.5*sum(eta4ed(ed4e).^2,2));
end