function x = i4_normal_ab ( a, b )

%*****************************************************************************80
%
%% i4_normal_ab() returns a scaled pseudonormal I4.
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    04 November 2021
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real A, the mean of the normal PDF.
%
%    real B, the standard deviation of the normal PDF.
%
%  Output:
%
%    integer X, a sample of the standard normal PDF.
%
  r1 = rand ( );
  r2 = rand ( );
  x = sqrt ( - 2.0 * log ( r1 ) ) * cos ( 2.0 * pi * r2 );

  x = round ( a + b * x );

  return
end
