11 December  2023   8:27:17.607 AM      
 
truncated_normal_rule():
  Fortran77 version
 
  For the (truncated) Gaussian probability density function
    pdf(x) = exp(-0.5*((x-MU)/SIGMA)^2) / SIGMA / sqrt ( 2 * pi )
  compute an N-point quadrature rule for approximating
    Integral ( A <= x <= B ) f(x) pdf(x) dx
 
  OPTION determines the truncation interval [A,B]:
  0: (-oo,+oo)
  1: [A,+oo)
  2: (-oo,B]
  3: [A,B]
 
  The user specifies OPTION, N, MU, SIGMA, A, B and FILENAME.
 
  FILENAME is used to generate 3 files:
 
    filename_w.txt - the weight file
    filename_x.txt - the abscissa file.
    filename_r.txt - the region file, listing A and B.
 
  Enter OPTION, 0/1/2/3:
 
  Enter the rule order N:
 
  Enter MU, the mean value of the normal distribution:
 
  Enter SIGMA, the standard deviation:
 
  Enter FILENAME, the "root name" of the quadrature files).
 
  OPTION =        0
  N =        5
  MU =    1.00000    
  SIGMA =    2.00000    
  A = -oo
  B = +oo
  FILENAME = "option0".
 
  Creating quadrature files.
 
  "Root" file name is   "option0".
 
  Weight file will be   "option0_w.txt".
  Abscissa file will be "option0_x.txt".
  Region file will be   "option0_r.txt".
 
truncated_normal_rule():
  Normal end of execution.
 
11 December  2023   8:27:17.616 AM      
11 December  2023   8:27:17.618 AM      
 
truncated_normal_rule():
  Fortran77 version
 
  For the (truncated) Gaussian probability density function
    pdf(x) = exp(-0.5*((x-MU)/SIGMA)^2) / SIGMA / sqrt ( 2 * pi )
  compute an N-point quadrature rule for approximating
    Integral ( A <= x <= B ) f(x) pdf(x) dx
 
  OPTION determines the truncation interval [A,B]:
  0: (-oo,+oo)
  1: [A,+oo)
  2: (-oo,B]
  3: [A,B]
 
  The user specifies OPTION, N, MU, SIGMA, A, B and FILENAME.
 
  FILENAME is used to generate 3 files:
 
    filename_w.txt - the weight file
    filename_x.txt - the abscissa file.
    filename_r.txt - the region file, listing A and B.
 
  Enter OPTION, 0/1/2/3:
 
  Enter the rule order N:
 
  Enter MU, the mean value of the normal distribution:
 
  Enter SIGMA, the standard deviation:
 
  Enter the left endpoint A:
 
  Enter FILENAME, the "root name" of the quadrature files).
 
  OPTION =        1
  N =        9
  MU =    2.00000    
  SIGMA =   0.500000    
  A =    0.00000    
  B = +oo
  FILENAME = "option1".
 
  Creating quadrature files.
 
  "Root" file name is   "option1".
 
  Weight file will be   "option1_w.txt".
  Abscissa file will be "option1_x.txt".
  Region file will be   "option1_r.txt".
 
truncated_normal_rule():
  Normal end of execution.
 
11 December  2023   8:27:17.619 AM      
11 December  2023   8:27:17.621 AM      
 
truncated_normal_rule():
  Fortran77 version
 
  For the (truncated) Gaussian probability density function
    pdf(x) = exp(-0.5*((x-MU)/SIGMA)^2) / SIGMA / sqrt ( 2 * pi )
  compute an N-point quadrature rule for approximating
    Integral ( A <= x <= B ) f(x) pdf(x) dx
 
  OPTION determines the truncation interval [A,B]:
  0: (-oo,+oo)
  1: [A,+oo)
  2: (-oo,B]
  3: [A,B]
 
  The user specifies OPTION, N, MU, SIGMA, A, B and FILENAME.
 
  FILENAME is used to generate 3 files:
 
    filename_w.txt - the weight file
    filename_x.txt - the abscissa file.
    filename_r.txt - the region file, listing A and B.
 
  Enter OPTION, 0/1/2/3:
 
  Enter the rule order N:
 
  Enter MU, the mean value of the normal distribution:
 
  Enter SIGMA, the standard deviation:
 
  Enter the right endpoint B:
 
  Enter FILENAME, the "root name" of the quadrature files).
 
  OPTION =        2
  N =        9
  MU =    2.00000    
  SIGMA =   0.500000    
  A = -oo
  B =    3.00000    
  FILENAME = "option2".
 
  Creating quadrature files.
 
  "Root" file name is   "option2".
 
  Weight file will be   "option2_w.txt".
  Abscissa file will be "option2_x.txt".
  Region file will be   "option2_r.txt".
 
truncated_normal_rule():
  Normal end of execution.
 
11 December  2023   8:27:17.621 AM      
11 December  2023   8:27:17.623 AM      
 
truncated_normal_rule():
  Fortran77 version
 
  For the (truncated) Gaussian probability density function
    pdf(x) = exp(-0.5*((x-MU)/SIGMA)^2) / SIGMA / sqrt ( 2 * pi )
  compute an N-point quadrature rule for approximating
    Integral ( A <= x <= B ) f(x) pdf(x) dx
 
  OPTION determines the truncation interval [A,B]:
  0: (-oo,+oo)
  1: [A,+oo)
  2: (-oo,B]
  3: [A,B]
 
  The user specifies OPTION, N, MU, SIGMA, A, B and FILENAME.
 
  FILENAME is used to generate 3 files:
 
    filename_w.txt - the weight file
    filename_x.txt - the abscissa file.
    filename_r.txt - the region file, listing A and B.
 
  Enter OPTION, 0/1/2/3:
 
  Enter the rule order N:
 
  Enter MU, the mean value of the normal distribution:
 
  Enter SIGMA, the standard deviation:
 
  Enter the left endpoint A:
 
  Enter the right endpoint B:
 
  Enter FILENAME, the "root name" of the quadrature files).
 
  OPTION =        3
  N =        5
  MU =    100.000    
  SIGMA =    25.0000    
  A =    50.0000    
  B =    150.000    
  FILENAME = "option3".
 
  Creating quadrature files.
 
  "Root" file name is   "option3".
 
  Weight file will be   "option3_w.txt".
  Abscissa file will be "option3_x.txt".
  Region file will be   "option3_r.txt".
 
truncated_normal_rule():
  Normal end of execution.
 
11 December  2023   8:27:17.624 AM