bug in multi-dimensional embedding solved

crp2.m

@@ -81,6 +81,9 @@ function xout=crp2(varargin)
 % $Revision$
 %
 % $Log$
+% Revision 5.6  2006/02/06 13:46:17  marwan
+% plugin for order patterns recurrence plots supported
+%
 % Revision 5.5  2005/11/23 07:30:30  marwan
 % modified interdependent algorithm
 % bug in showing RP fixed
@@ -225,7 +228,7 @@ if isnumeric(varargin{1}) 		% read commandline input
   end
 
     Nx=size(x,1); Ny=size(y,1);
-    NX=Nx-t*(m-1);NY=Ny-t*(m-1);  
+    NX=Nx-t*(m0-1);NY=Ny-t*(m0-1);  
     x0=zeros(Nx,m);y0=zeros(Ny,m);
     x0(1:size(x,1),1:size(x,2))=x; 
     y0(1:size(y,1),1:size(y,2))=y; 
@@ -237,12 +240,12 @@ if isnumeric(varargin{1}) 		% read commandline input
 
   if ~isempty(find(isnan(x))), for k=1:size(x,2),  x(find(isnan(x(:,k))),:)=[]; end, end
   if ~isempty(find(isnan(y))), for k=1:size(y,2),  y(find(isnan(y(:,k))),:)=[]; end, end
-  if size(x,1) < t*(m-1)+1 | size(y,1) < t*(m-1)+1
+  if size(x,1) < t*(m0-1)+1 | size(y,1) < t*(m0-1)+1
      error(['Too less data',10,...
             'Either too much NaN or the number of columns in the vectors do not match.'])
   end
-  if method==7 & m > 1, 
-      m=1; 
+  if method==7 & m0 > 1, 
+      m0=1; 
       disp('Warning: For order matrix a dimension of one is used.')
   end
   ds=eye(m);