NAME

Math::Random::Cauchy - Random numbers following a Cauchy PDF

SYNOPSIS

  use Math::Random::Cauchy;
  my $cauchy = Math::Random::Cauchy->new(
    fwhm  => 1,  # the width (full width, half maximum), default==1
    middle => 5, # the expectation value, default==0
    random => 'rand', # use Perl's builtin (default)
  );
  
  foreach (1..100) {
      my $rnd = $cauchy->rand();
      # ...
  }
  
  # Use Math::Random::MT instead of bultin rand()
  use Math::Random::MT;
  my $mt = Math::Random::Mt->new($seed);
  $cauchy = Math::Random::Cauchy->new(
    random => sub { $mt->rand() };
  );

DESCRIPTION

This module transforms uniformly spaced random numbers into random numbers that follow the Cauchy Probability Density Function (PDF).

A more general transformation method is implemented in Math::Random::AcceptReject.

The algorithm is from Blobel et al as quoted in the SEE ALSO section below.

METHODS

new

Creates a new random number generator. Takes named arguments.

Optional parameters:

  random: The random number generator. Defaults to using Perl's
          rand() function. May be set to either 'rand' for the
          default or a subroutine reference for custom random
          number generators. Expected to return one or more(!)
          random numbers per call.
  fwhm:   Full width, half maximum. Defaults to 1.
  middle: Expectation value for x. Defaults to 0.

rand

Returns the next random number of Cauchy PDF.

SEE ALSO

Math::Random::MT, Math::Random, Math::Random::OO, Math::TrulyRandom, Math::Random::AcceptReject

The examples in the examples/ subdirectory of this distribution.

The algorithm was taken from: (German)

Blobel, V., and Lohrmann, E. Statistische und numerische Methoden der Datenanalyse. Stuttgart, Leipzig: Teubner, 1998

AUTHOR

Steffen Mueller, <smueller@cpan.org>

COPYRIGHT AND LICENSE

Copyright (C) 2007 by Steffen Mueller

This library is free software; you can redistribute it and/or modify it under the same terms as Perl itself, either Perl version 5.6 or, at your option, any later version of Perl 5 you may have available.