diff --git a/analysis.m b/analysis.m
index 27bfe2a..458e416 100644
--- a/analysis.m
+++ b/analysis.m
@@ -76,10 +76,35 @@ for i=1:length(param.dates)
     report.data(num).EP_Etemp(i)=mymean(Etemp,1e10);
 end
 
+report.data(num).RPA_HHeden_lim=mylim(report.data(num).RPA_HHeden,param);
+report.data(num).RPA_O2den_lim=mylim(report.data(num).RPA_O2den,param);
+report.data(num).EP_Eden_lim=mylim(report.data(num).EP_Eden,param);
+report.data(num).EP_Etemp_lim=mylim(report.data(num).EP_Etemp,param);
+
+
 param.report=report;
 save(fname_rep,'report');
 end
 
+function out=mylim(y,param)
+Q1=quantile(y,0.25);
+Q2=quantile(y,0.5);
+Q3=quantile(y,0.75);
+L1=Q2-param.k*(Q3-Q1);
+L2=Q2+param.k*(Q3-Q1);
+
+out=[];
+for i=1:length(y)
+    if (y(i)<=L1)||(y(i)>=L2)&&(isnan(y(i)))        
+        out(i)=(y(i)-Q2)*100/(L1-Q2);% отклонение в процентах
+        %out(i)=1;
+    else
+        out(i)=0;
+    end
+end
+end
+
+
 function num=give_number(report,f,R)
 flag=0; % Проверяет, есть ли в отчете данные со спутника f
 [~,b]=size(report.data);