diff --git a/twinsurr.m b/twinsurr.m
index 7b52ca53c6cb865649ec2c363bbece6010bee8f0..9f35fd4135bb58347777c70d05e33139c93fcd7c 100644
--- a/twinsurr.m
+++ b/twinsurr.m
@@ -33,6 +33,9 @@ function y = twinsurr(varargin)
 % $Revision$
 %
 % $Log$
+% Revision 5.1  2008/07/01 13:09:27  marwan
+% initial import
+%
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% programme properties
 
@@ -40,6 +43,7 @@ global errcode props
 
 init_properties
 nsur_init = 100;
+sil = 0;
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% check the input
 
@@ -73,21 +77,26 @@ i_double=find(cellfun('isclass',varargin,'double'));
 i_char=find(cellfun('isclass',varargin,'char'));
 check_meth={'ma','eu','mi','nr','rr','fa','in','om','op','di'}; 	% maxnorm, euclidean, nrmnorm,  fan, distance
 check_norm={'non','nor'};                        % nonormalize, normalize
+check_sil={'ve','si'};                         % verbose, silent
 
 if isnumeric(varargin{1}) 		% read commandline input
    % check the text input parameters for method, gui 
     temp_meth=0;
     temp_norm=0;
+    temp_sil=0;
     if ~isempty(i_char)
          for i=1:length(i_char), 
             varargin{i_char(i)}(4)='0';
             temp_norm=temp_norm+strcmpi(varargin{i_char(i)}(1:3),check_norm'); 
             temp_meth=temp_meth+strcmpi(varargin{i_char(i)}(1:2),check_meth'); 
+            temp_sil=temp_sil+strcmpi(varargin{i_char(i)}(1:2),check_sil'); 
          end
          method_n=min(find(temp_meth));
          nonorm=min(find(temp_norm))-1;
+         sil=min(find(temp_sil))-1;
          for i=1:length(i_char); temp2(i,:)=varargin{i_char(i)}(1:3); end
 
+         if isempty(sil), sil=0; end
          if isempty(nonorm), nonorm=1; end
          if nonorm>1, nonorm=1; end
          if isempty(method_n), method_n=1; end
@@ -117,9 +126,9 @@ else
 end
 
 
-N = length(x); 
 if size(x,1) < size(x,2), x = x'; end
 m = size(x,2);
+N = length(x); 
 
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% make surrogates
@@ -128,12 +137,12 @@ X = crp(x,m,t,e,method,norm_str,'sil');
 NX = length(X);
 
 %% find twins
-h = waitbar(0,'Searching Twins');
-for i = 1:NX, waitbar(i/NX)
+if ~sil, h = waitbar(0,'Searching Twins'); end
+for i = 1:NX, if ~sil, waitbar(i/NX); end
     A = repmat(X(:,i),1,NX);
     S{i} = find(all(X == A));
 end
-delete(h)
+if ~sil, delete(h); end
 
 for k = 1:nsur
     %% chose randomly first point