relative_error.m

best routing protocol

function rerr = relative_error(estimate, real)
      rerr=abs(estimate-real)./(min(estimate,real)+0.1); %0.1 to avoid divided by zero
      mask  = (real>0);
      rerr= rerr.*mask;
      rerr = rerr(:);