function [out]=relation(e,ec,u,rule)
n1=length(rule);
n2=length(e);
n3=length(ec);
n4=length(u);
l=n2*n3;
out=zeros(l,n4);
for k=1:n1
    v1=e(rule(k,1),:);
    v2=ec(rule(k,2),:);
    v3=u(rule(k,3),:);
    d=v1mutiv2(v1,v2);
    v=mtov(d);
    d=v1mutiv2(v,v3);
    out=max(out,d);
end