function R=isamicable(M,N)
if norm(M.*N) > 0
 R = 0;
elseif M*N' == N*M'
 R = 1;
elseif M*N' == -N*M'
 R = -1;
else
 R = NaN;
end
