function mplot(A,ind)
% pairwise scatterplots of the selected columns of A
clf;
[p n]=size(A);
n=length(ind); 
a=min(min(A));b=max(max(A));
for i=1:n
    subplot(n,1,i)
    plot(A(:,ind(i)))
    hold on
    axis([0 p a b]);
    ylabel(['A',int2str(ind(i))]);
    plot([0,p+1],[0,0],'r--')
    hold off;
end;
%