function da=cosmosf2(t,a)
global Omega

f=Omega./a+1-Omega;
%if f<0
%    da=-sqrt(f);
%else
%    da=sqrt(f);
%end
da=-sqrt(f);