function da=univf(t,a)
global Omegar Omegam OmegaL Omega0 alpha

f=Omegar./a.^2+Omegam./a+OmegaL*a.^2+1-Omega0;

da=alpha*sqrt(f);