Matlab - como llamar archivos de entrada

function [covmat, eofvec, eofval, eofvar, power, phase, powert, phaset] = eof(x, iop)

%______________________________________________________________________________

%

%[COVMAT, EOFVEC, EOFVAL, EOFVAR, MODES] = EOF(X, 0) real case

%

%[COVMAT,EOFVEC,EOFVAL,EOFVAR,POWER,PHASE,POWERT,PHASET] = EOF(X, [1 or 2]) complex case

%______________________________________________________________________________

%

%                               EOF

%            Calculates Empirical Orthogonal Functions

%

%

% X        : Time Series Matrix

% IOP      : Real(iop = 0) or Complex (iop = 1, 2) EOF's

%

% COVMAT   : Real or Complex Covariance Matrix.

% EOFVEC   : Real or Complex Eigenvectors Matrix (N by N).

% EOFVAL   : Real or Complex Eigenvalues Vector (N).

% EOFVAR   : Real or Complex Relative Variance Vector (N).

% POWER    : Spatial Power Functions Matrix (N by N).

% PHASE    : Spatial Phase Functions Matrix (N by N).

% POWERT   : Temporal Power Functions Matrix (M by N).

% PHASET   : Temporal Phase Functions Matrix (M by N).

%

%

% JJSP 1995

%

% help when called without input arguments

if nargin == 0

 help eof

 return

end

tfac = 0.50 / 86400.00; rfac = 180.0 / pi;

covmat = []; eofvec = []; eofval = []; eofvar = []; power = [];

phase = []; powert = []; phaset = [];

x = detrnd(x);

[n, m] = size(x);

ind_nan = find(isnan(x) == 1);

if ~isempty(ind_nan)

 x(ind_nan) = zeros(size(ind_nan));

end

%

% IOP = 0 Real Covariance Matrix and EOF's

%

if iop == 0

 covmat = (1 / n) * x' * x;

 [eofvec, eofval, eofvar] = reig(covmat);

 powert = []; phase = []; phaset;

 power = x * eofvec

size(power)

%

% IOP = 1 Complex Covariance Matrix and CEOF's

%

elseif iop == 1

 x = hilbert(x);

 covmat = (1 / n) * x' * x;

 [eofvec, eofval, eofvar, power, phase] = cpxeig(covmat);

 y = x * eofvec;

 powert = abs(y);

 phaset = angle(y);

 phaset = unwrap(phaset) * rfac;

 clear y

%

% IOP = 2 Complex Covariance Matrix and CEOF's

%

elseif iop ==2

 x = x(:, 1: 2: (n-1)) + i .* x(:, 2: 2: n);

 [n, m] = size(x);

 covmat = (1 / n) * x' * x;

 [eofvec, eofval, eofvar, power, phase] = cpxeig(covmat);

 y = x * eofvec;

 powert = abs(y);

 phaset = angle(y);

 phaset = unwrap(phaset) * rfac;

 clear y

 [xmax, ind] = max(power);

 [n, m] = size(power);

 for k = 1: m

  power(:, k) = power(:, k) / xmax(k);

  z(:, k) = phase(:, k) - phase(ind(k), k);

 end

 phase = z;

end

>> x=load('2013Dic1.mat'); iop=0;  [covmat, eofvec, eofval, eofvar, power, phase, powert, phaset] = eof(x, iop)

Error using cast

Unsupported data type for conversion: 'struct'.

Error in detrend>detrendInternal (line 138)

        W = cast(W,'like',x);

Error in detrend (line 89)

    y = detrendInternal(x,polyDeg,bp,s,continuity,N);

Error in eof (line 30)

x = detrend(x);