function SMALL_plot(SMALL)
  % Ivan Damnjanovic 2009
  % Function gets as input SMALL structure and plots  the solution and
  % reconstructed signal

  figure;
 
  m=size(SMALL.solver,2);
  n=size(SMALL.solver(1).reconstructed,2)+1;
  for i =1:m
  
  subplot(m,n, (i-1)*n+1); plot(1:length(SMALL.solver(i).solution), SMALL.solver(i).solution, 'b')
        title([ SMALL.solver(i).name,'(', SMALL.solver(i).param,')'])
          xlabel('Coefficient') 
     
  % Plot reconstructed signal against original

  
  for j=2:n
      
      subplot(m,n,(i-1)*n+j); plot(1:length(SMALL.solver(i).reconstructed(:,j-1)), SMALL.solver(i).reconstructed(:,j-1)  ,'b.-',  1:length(SMALL.Problem.signal(:,j-1)), SMALL.Problem.signal(:,j-1),'r--')
          %legend(SMALL.solver(i).name,'Original signal',0);
              title('Reconstructed and original signals');
  end
end
