for mq1=mq
  subplot(length(mq),1,find(mq==mq1))
  axis([-1 1 -0.5 1])
end

shg;
clear afterframe