Найти значение по индексу без потери местонахождения в векторе

А4=GUnic(:,1);% Месяцы
 k4 = find(A4==4);% берем 4 если он есть 
 AN4=GUnic(:,2);% значение в соответствующих месяцах
Nes4=AN4(k4);% получили значения всех значений в 4 месяце
 
А5=GUnic(:,1);% Месяцы
 k5 = find(A5==5);% берем 5 если он есть 
 AN5=GUnic(:,2);% значение в соответствующих месяцах
Nes5=AN5(k5);% получили значения всех значений в 5 месяце

%GUnic файл из приложения, если добавить колонку с годами впереди, будет работать и по годам 
А4=GUnic(:,1);% Месяцы
 k4 = find(A4==4);% берем 4 если он есть 
 AN4=GUnic(:,2);% значение в соответствующих месяцах
Nes4=AN4(k4);% получили значения всех значений в 4 месяце
 
А5=GUnic(:,1);% Месяцы
 k5 = find(A5==5);% берем 5 если он есть 
 AN5=GUnic(:,2);% значение в соответствующих месяцах
Nes5=AN5(k5);% получили значения всех значений в 5 месяце
% Далее все месяцы получим

txtfiles = dir( '*.txt' );
files={txtfiles(~[txtfiles.isdir]).name};
ss=1:size(files,2);
Nss=ss(1,end);
Nss1=Nss/2;
N=1;
for Nn=1:N;
for fnc = 1:ss(1,Nn);
filenameP1=[num2str(txtfiles(fnc).name) ];
ncidP1=importdata(files{fnc});
%end
A=ncidP1.data;
%A=importdata('res/d0424469.txt');
data1=A;
F=data1(:,[1,2,3,6,7,8,10]);
%r=find(all(F<9999,2) & data1(:,4)>=0);
r=find(all(F<9999,2) & data1(:,4)>=0 & data1(:,2)>=4); % с 4 мес
F1=F(r,:);
F1(:,8)=F1(:,4).*(F1(:,4)-F1(:,7)).*(F1(:,5)<2.5);
Y=F1(:,1);
miny=min(Y);
M=F1(:,2)+12*(Y-miny);
T=F1(:,8);
Summa=splitapply(@sum,T,findgroups(M));
f=find(diff(M)~=0)+1;
K=T(f);
Summa(2:end)=Summa(2:end)+Summa(1:end-1).*(K~=0);
UM=unique(M);
C=fix(UM/12);
Month=mod(UM,12);
Year=miny+C;
mz=find(Month==0);
Year(mz)=Year(mz)-1;
Month(mz)=12;
ff=sort([1;f;find(T==0);length(T)+1]);
dff=diff(ff);
A=mat2cell(T,dff,1);
CSum=cell2mat(cellfun(@cumsum,A,'UniformOutput',false));
F2=[F1 CSum];
F3=F2(CSum>1000,:);
X=(F3(:,1));
Ygod=F3(:,1);
miny1=min(Ygod);maxy2=max(Ygod);
G=(maxy2-miny1)+1;
S2=histcounts(X,G);
%dlmwrite(' filenameP1', S2,'\t');
dlmwrite([num2str(filenameP1) 'G.txt'],S2);
end
end
 
 
Az=Ygod;
d=diff(Az)~=0;
cs=cumsum([0 d']);
g=groupcounts(cs');% вектор из 59 значений, за каждый год 
god1=miny;
% для объяснения алгоритма
   g1=g(1);% первый год получили количество значений 
   g2=g(2);% второй год общее количество значений
   g3=g(3);
   g4=g(4);
 
% нарезка интервалов месяцев
   XX1=(F3(1:g1,2));
   XX2=(F3(g1+1:g1+g2,2));
% нарезка интервалов лет
 
GOD1=(F3(1:g1,1));
   GOD2=(F3(g1+1:g1+g2,1));
% получили годы все
 
 
 AzXX1=XX1; % за каждый год получить сколько есть значений в  каждом месяце и уникальные месяцы, но если один месяц выпадает, то значение уже не совпадает с уникальным месяцем
 dXX1=diff(AzXX1)~=0;
  csXX1=cumsum([0 dXX1']);
   gXX1=groupcounts(csXX1');
    Cuniq1 = unique(XX1);
     C1 = cat(1,Cuniq1,gXX1);
     
     
     AzXX2=XX2;
dXX2=diff(AzXX2)~=0;
 csXX2=cumsum([0 dXX2']);
  gXX2=groupcounts(csXX2');
   Cuniq2 = unique(XX2);
    C2 = cat(1,Cuniq2,gXX2);
 CunikALL = cat(1,Cuniq1    ,Cuniq2 ,Cuniq3 ,Cuniq4 ,Cuniq5 ,Cuniq6 ,Cuniq7 ,Cuniq8 ,Cuniq9 ,Cuniq10    ,Cuniq11    ,Cuniq12    ,Cuniq13    ,Cuniq14    ,Cuniq15    ,Cuniq16    ,Cuniq17    ,Cuniq18    ,Cuniq19    ,Cuniq20    ,Cuniq21    ,Cuniq22    ,Cuniq23    ,Cuniq24    ,Cuniq25    ,Cuniq26    ,Cuniq27    ,Cuniq28    ,Cuniq29    ,Cuniq30    ,Cuniq31    ,Cuniq32    ,Cuniq33    ,Cuniq34    ,Cuniq35    ,Cuniq36    ,Cuniq37    ,Cuniq38    ,Cuniq39    ,Cuniq40    ,Cuniq41    ,Cuniq42    ,Cuniq43    ,Cuniq44    ,Cuniq45    ,Cuniq46    ,Cuniq47    ,Cuniq48    ,Cuniq49    ,Cuniq50    ,Cuniq51    ,Cuniq52    ,Cuniq53    ,Cuniq54    ,Cuniq55    ,Cuniq56    ,Cuniq57    ,Cuniq58);
 gXXALL = cat(1,gXX1    ,gXX2   ,gXX3   ,gXX4   ,gXX5   ,gXX6   ,gXX7   ,gXX8   ,gXX9   ,gXX10  ,gXX11  ,gXX12  ,gXX13  ,gXX14  ,gXX15  ,gXX16  ,gXX17  ,gXX18  ,gXX19  ,gXX20  ,gXX21  ,gXX22  ,gXX23  ,gXX24  ,gXX25  ,gXX26  ,gXX27  ,gXX28  ,gXX29  ,gXX30  ,gXX31  ,gXX32  ,gXX33  ,gXX34  ,gXX35  ,gXX36  ,gXX37  ,gXX38  ,gXX39  ,gXX40  ,gXX41  ,gXX42  ,gXX43  ,gXX44  ,gXX45  ,gXX46  ,gXX47  ,gXX48  ,gXX49  ,gXX50  ,gXX51  ,gXX52  ,gXX53  ,gXX54  ,gXX55  ,gXX56  ,gXX57  ,gXX58);
GUnic=cat(1,CunikALL,gXXALL);
k4 = find(A111==4);
AN=GUnic(:,2);
% также и с позицией по году
GOD1=(F3(1:g1,1));
   GOD2=(F3(g1+1:g1+g2,1));
AzGOD1=GOD1;
     dGOD1=diff(AzGOD1)~=0;
     csGOD1=cumsum([0 dGOD1']);
     gGOD1=groupcounts(csGOD1');
     CuniqG = unique(GOD1);
     CGOD1 = cat(2,CuniqG,gGOD1);
  AzGOD2=GOD2;
  dGOD2=diff(AzGOD2)~=0;
  csGOD2=cumsum([0 dGOD2']);
  gGOD2=groupcounts(csGOD2');
  CuniqG2 = unique(GOD2);
  C2 = cat(2,CuniqG2,gGOD2);
 
  CunikALLGOD = cat(1,CuniqG    ,CuniqG2    ,CuniqG3    ,CuniqG4    ,CuniqG5    ,CuniqG6    ,CuniqG7    ,CuniqG8    ,CuniqG9    ,CuniqG10   ,CuniqG11   ,CuniqG12   ,CuniqG13   ,CuniqG14   ,CuniqG15   ,CuniqG16   ,CuniqG17   ,CuniqG18   ,CuniqG19   ,CuniqG20   ,CuniqG21   ,CuniqG22   ,CuniqG23   ,CuniqG24   ,CuniqG25   ,CuniqG26   ,CuniqG27   ,CuniqG28   ,CuniqG29   ,CuniqG30   ,CuniqG31   ,CuniqG32   ,CuniqG33   ,CuniqG34   ,CuniqG35   ,CuniqG36   ,CuniqG37   ,CuniqG38   ,CuniqG39   ,CuniqG40   ,CuniqG41   ,CuniqG42   ,CuniqG43   ,CuniqG44   ,CuniqG45   ,CuniqG46   ,CuniqG47   ,CuniqG48   ,CuniqG49   ,CuniqG50   ,CuniqG51   ,CuniqG52   ,CuniqG53   ,CuniqG54   ,CuniqG55   ,CuniqG56   ,CuniqG57   ,CuniqG58);

A=[4 1
5 4
6 8
7 12
8 8
9 5
10 12
11 21
12 2
5 12
6 4
9 3
10 15
11 3
4 16
5 8
7 13
8 7
9 11
10 9
11 5];
a=A(:,1);
V=[0;cumsum(diff(a)<0)];
Z=(min(a):max(a))';
m=length(Z);
n=V(end)+1;
B=repmat(Z,n,1);
[~,k]=ismember(a,B);
N=repmat(NaN,length(B),1];
N(k+m*V)=A(:,2);
M=[B,N];