#!/usr/bin/perl
#
# compute the asymptotic behavior of a beta-model
#    cluster, showing the observed counts as a 
#    function of distance at large distances
#
# MWR 4/12/2011

use POSIX;
$debug = 1;

$core_radius = 1.0;



$start_r = 2*$core_radius;
$end_r = 20*$core_radius;
$d_r = 1.0*$core_radius;

# first, we do the calc for beta = 0.8
$K = 0.5E23;
$beta = 0.8;
$power = -3.0*$beta + 0.5;
for ($r = $start_r; $r < $end_r; $r += $d_r) {

  $y = $K*(1.0 + ($r*$r))**($power);
  printf " %9.4f %9.4e \n", $r, $y;

}

printf "\n\n\n";

# now, we do the calc for beta = 1.2
$K = 0.3E22;
$beta = 1.2;
$power = -3.0*$beta + 0.5;
for ($r = $start_r; $r < $end_r; $r += $d_r) {

  $y = $K*(1.0 + ($r*$r))**($power);
  printf " %9.4f %9.4e \n", $r, $y;

}

printf "\n\n\n";










exit 0;