function [ppomer] = CO2aCO(t,X_CO2,X_CO,tau,out,k,meze)
% [t_CO,CO,l]=vyhazet_radky(t,X_CO);
% [t_CO2,CO2,l]=vyhazet_radky(t,X_CO2);
[tspol,CO2spol,COspol,l]=vyhazet_radky2(t,X_CO2,X_CO);
n=length(tspol);
for i=1:n
    pomer(i)=CO2spol(i)/COspol(i); 
end;
pomer=pomer';
if out==1;
    [Q,I]=outliers([tspol pomer],1.5);
    t=Q(:,1);
    X=Q(:,2);
else
    t=tspol;
    X=pomer;
end;
    
if k==1;
    b_lin=robustfit(t,X);
    y=b_lin(1)+b_lin(2)*tau;
end;    
if k==2;
    b_moc=robustfit(log(t),log(X));
    y=exp(b_moc(1))*(tau^b_moc(2));
end; 
if k==3;
    f1=fit(t,X,'poly1');
    b5=coeffvalues(f1);
    y=b5(2)+b5(1)*tau;
end;
if k==4;
    f2=fit(t,X,'exp1');
    b2=coeffvalues(f2);
    y=b2(1)*exp(b2(2)*tau);
end;
if k==5;
    r=ksr(t,X);
    f3=fit(t,r.f','poly2');
    b3=coeffvalues(f3);
    y=b3(1)*tau^2+b3(2)*tau+b3(3);
end;
if k==6;
    klp=smooth(X,3);
    f4=fit(t,klp,'poly2');
    b4=coeffvalues(f4);
    y=b4(1)*tau^2+b4(2)*tau+b4(3);
end;
ppomer=y;

if ppomer >= meze(2)
verdikt='Papírová izolace je naprosto v pořádku';
end;
if ppomer >=meze(1) && ppomer <meze(2)
verdikt='Papírová izolace je zřejmě v pořádku';
end;
if ppomer >=meze(3) && ppomer <meze(1)
verdikt='Papírová izolace je zřejmě v pořádku';
end;



end

