function [ e, mu, azm, x, y, z ] = neutron_source ( )

%*****************************************************************************80
%
%% neutron_source() generates a new particle from the neutron source.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    02 September 2012
%
%  Author:
%
%    Original FORTRAN77 version by Kahaner, Moler, Nash.
%    This version by John Burkardt.
%
%  Reference:
%
%    David Kahaner, Cleve Moler, Steven Nash,
%    Numerical Methods and Software,
%    Prentice Hall, 1989,
%    ISBN: 0-13-627258-4,
%    LC: TA345.K34.
%
%  Output:
%
%    real E, the initial energy of the particle.
%
%    real MU, the cosine of the angle between the
%    particle's direction and the X axis.
%
%    real AZM, the azimuthal angle of the particle's
%    direction.
%
%    real X, Y, Z, the initial coordinates of the particle.
%
  mu = rand ( 1, 1 );
  azm = 2.0 * pi * rand ( 1, 1 );

  x = 0.0;
  y = 0.0;
  z = 0.0;

  e = energy ( );

  return
end

