function compm(m1,m2)
%

plot(m1(60:84,60:84),m2(60:84,60:84),'+g', ...
     m1(60:84,10:40),m2(60:84,10:40),'.r', ...
     m1(10:40,60:84),m2(10:40,60:84),'+b', ...
     m1(10:40,10:40),m2(10:40,10:40),'.c');