function [pomer] = Pomer_CO2aCO(t,X_CO2,X_CO,k,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);
if out==1;
    [Q,I]=outliers([t1 X1],1.5);
    t=Q(:,1);
    X=Q(:,2);
else
    t=t1;
    X=X1;
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);



end

