;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;; PROCEDURE OVST      ;;;
;;; By Darrin Casebeer  ;;;
;;;    Kevin Tubbs      ;;;
;;;    Summer 98        ;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;

PRO ovst,OM,OL,t
  COMMON COSMO
  Ho=60.0
  C=3.1e19/(3.16e7*60)
  deltaz=0.01

  znow=DBLARR(LIMITS)
  z=DBLARR(LIMITS)
  R=DBLARR(LIMITS)
  dR=DBLARR(LIMITS)
  dt=DBLARR(LIMITS)
  mu=DBLARR(LIMITS)


  holder=0.01


  FOR I=0,LIMITS-1 DO BEGIN
    znow=0.00
    z(I)=holder
    holder=holder+deltaz
  ENDFOR


  t=QSIMP('dt',0.00,z)
  t=C*t
  R=1/(1+z)
  FOR I=1,LIMITS-2 DO BEGIN
    dR[I]=R[I+1]-R[I-1]
    dt[I]=t[I+1]-t[I-1]
  ENDFOR
  FOR I=1,LIMITS-2 DO BEGIN
    OM[I]=R[I]^(-1.0)*(dR[I]/dt[I])^(-2.0)
    OL[I]=R[I]^2*(dR[I]/dt[I])^(-2.0)
  ENDFOR
  OM[0]=OM[1]
  OM[LIMITS-1]=OM[LIMITS-2]
  C=Omegamass/OM[0]
  OM=C*OM
  OL[0]=OL[1]
  OL[LIMITS-1]=OL[LIMITS-2]
  C=OmegaLambda/OL[0]
  OL=C*OL
END