#!/usr/bin/perl
#
# Compute values of the angular diameter of an object
#   for a range of distances, and for a range of 
#   Omega_M values.  We assume that Omega_lambda = 0.
#
# MWR 5/9/2011

$debug = 1;

$H = 1.0;
$c = 1.0;
$S = 1.0;

$start_z = 0.01;
$end_z = 100.0;
$mult_z = 1.05;

@omega_list = (0.1, 0.5, 1.0, 2.0, 10.0);

for ($z = $start_z; $z <= $end_z; $z *= $mult_z) {

  printf " %6.4e ", $z;
  foreach $omega (@omega_list) {


     $top = (($S*$H)/(2*$c)) * ((1 + $z)*(1 + $z)) * ($omega*$omega);
     $bot = $omega*$z - (2.0 - $omega)*( sqrt(1 + $omega*$z) - 1.0);
     $theta = $top / $bot;
     printf " %9.4e ", $theta;

  }
  printf "\n";

}







exit 0;
