view release on metacpan or  search on metacpan

cpt/differential_comoving_volume.pl  view on Meta::CPAN

my $lv_b = log10( $v_b );
my $lv_c = log10( $v_c );
my $lv_d = log10( $v_d );
my $lv_e = log10( $v_e );

# plot the graph
my $win = PDL::Graphics::PGPLOT::Window->new();

$win->env( log10(0.03), log10(100), log10(0.001), log10(20), 0, 30 );
$win->label_axes( "redshift z",
		  "comoving volume derivative H\\d0\\u\\u3\\d dV/(dzd\\gW)", 
		  "Fig 6 of Carroll, Press & Turner ARAA 1992 vol 30, 499-542" );

$win->hold;

$win->line( $lz, $lv_a, { LINESTYLE => 'solid' } );
$win->line( $lz, $lv_b, { LINESTYLE => 'dashed' } );
$win->line( $lz, $lv_c, { LINESTYLE => 'dot-dash' } );
$win->line( $lz, $lv_d, { LINESTYLE => 'dotted' } );
$win->line( $lz, $lv_e, { LINESTYLE => 'dash-dot-dot' } );

hogg/comoving_volume.pl  view on Meta::CPAN

# - ie the differential comoving volume element
my $d_a = $a->dcomov_vol( $z );
my $d_b = $b->dcomov_vol( $z );
my $d_c = $c->dcomov_vol( $z );

# plot the graph
my $win = PDL::Graphics::PGPLOT::Window->new();

$win->env( 0, 5, 0, 1.1 );
$win->label_axes( "redshift z",
		  "comoving volume element [1/D\\dH\\u]\\u3\\d dV\\dc\\u/dz/d\\gW",
		  "Fig 5 of Hogg, D.W. astro-ph/9905116" );

$win->hold;

$win->line( $z, $d_a, { LINESTYLE => 'solid' } );
$win->line( $z, $d_b, { LINESTYLE => 'dotted' } );
$win->line( $z, $d_c, { LINESTYLE => 'dashed' } );

$win->release;
$win->close;