function plotmst(mst,tg,m,nhd,lntp)
if nargin <3,m=1;end
if nargin <4,nhd=1;end
if nargin <5,lntp='-';end

global ubnd nf;

figure(m),
pn = 0;
    if pn == -1
       subplot(211), plot(tg, mst(1,:),lntp);
       if nhd ==1, hold on; end
       subplot(212), plot( tg,ubnd * exp(-exp(-mst(3,:)/10)),lntp);
       if nhd ==1, hold on; end
    elseif pn == 0
       subplot(311),
       for j = 1:nf
           plot(tg, mst(1+(j-1)*3,:),lntp);
           if nhd ==1, hold on; end
       end
       ylabel('amplitude')
       subplot(312), 
       for j = 1:nf
           plot( tg, mst(2+(j-1)*3,:),lntp);     
           if nhd ==1, hold on; end   
       end
       ylabel('phase')
       subplot(313), 
       for j = 1:nf
           plot( tg, mst(3+(j-1)*3,:),lntp);     
           if nhd ==1, hold on; end   
       end
       ylabel('frequency')       
    elseif pn == 1
       subplot(211),   plot(tg, sqrt(mst(1,:).^2 + mst(2,:).^2) ,lntp);
       if nhd ==1, hold on; end
       subplot(212), plot( tg,ubnd * exp(-exp(-mst(3,:)/10)),lntp);
       if nhd ==1, hold on; end
    elseif pn == 2
        subplot(211),plot(tg, mst(1,:))
        subplot(212), plot(tg, ubnd * exp(-exp(-mst(2,:)/10)));
    elseif pn == 3        
       subplot(211), plot(tg, sqrt(mst(1,:).^2 + mst(2,:).^2) ,lntp);
       if nhd ==1, hold on; end
       subplot(212), plot( tg,mst(3,:),lntp);
       if nhd ==1, hold on; end
        
    end