4-х мерная интерполяция

clc; clear all; clf;
 
rm = importdata('Эксперимент.txt');
Mrm = rm.data (:,1);
alpharm = rm.data (:,2);
deltarm = rm.data (:,3);
phirm = rm.data (:,4);
 
X=importdata('Массивы.txt');
Mt = X(:,1); % t-значит табличное
alphat = X(:,2);
deltat = X(:,3);
phit = X(:,4);
Cxt = X(:,5);
 
[Mt] = unique(Mt);
[alphat] = unique(alphat);
Ml = length(Mt);
alphal = length(alphat);
 
s=0;
n=0;
a =NaN(1);
M =NaN(1);
 
for i=1:size (X,1)/247
   x{i}=X(247*(i-1)+1:247*i,:);
   
   Mt = x{i}(:,1); % t-значит табличное
   alphat = x{i}(:,2);
   deltat = x{i}(:,3);
   phit = x{i}(:,4);
   Cxt = x{i}(:,5);
   
   [Mt] = unique(Mt);
   [alphat] = unique(alphat);
   [deltat] = unique(deltat);
   [phit] = unique(phit);
   [Cxt] = reshape(Cxt,length(unique(deltat)),[]);
  
   if ([a]~=[alphat])&(s <= alphal)
    s= s+1
    [a]=[alphat]
   end;
   
   if  s > alphal
    s=1
   end;
    
   if ([M]~=[Mt])&(n <= Ml) 
    n= n+1
    [M]=[Mt]
   end;
   
   if  n > Ml 
    n=1     
   end;
   
   Cx(:, :, s, n) = [Cxt];
   
end
 
Cxrm = interpn(deltat,phit,alphat,Mt,Cx, deltarm, phirm,alpharm,Mrm, 'linear');

[dt,pt,at,mt] = ndgrid(unique(deltat),unique(phit),...
    unique(alphat),unique(Mt));
 
mCxt = x (:,5);
Cxt=reshape(mCxt,19,13,5,8);
 
for i=1:length(Mrm)
    if ~isnan(phirm(i))
        Cxrm(i,:) = interpn(dt,      pt,      at,    mt, Cxt, ...
           deltarm(i), phirm(i), alpharm(i), Mrm(i),'linear');
        disp([Mrm(i) alpharm(i) deltarm(i) phirm(i) Cxrm(i)])
    end
end
 
fid = fopen('res.txt', 'w');
if fid == -1
    error('File is not opened');
end
fprintf(fid, 'Mrm      alpharm deltarm phirm Cxrm\r\n');
fprintf(fid, '%.5f %.5f %.5f %.5f %.7f\r\n', [Mrm alpharm deltarm phirm Cxrm]');
fclose(fid);