#!/bin/csh
#
#  S M A L L   S I M U L A T I O N   P R O G R A M  
#
echo Copy source code to demo_triangle.f
#
cat >demo_triangle.f <<EOF
      PROGRAM DEMO_HBOOK
      IMPLICIT NONE
      INTEGER PAW_MEMORY_SIZE, NBIN, I
      REAL PAW_SPACE, LOW, UPP, SUM, RNDM, WEIGHT
      PARAMETER (PAW_MEMORY_SIZE=10000)
      COMMON /PAWC/ PAW_SPACE(PAW_MEMORY_SIZE)
*
      CALL HLIMIT(PAW_MEMORY_SIZE)      ! Tell HBOOK the size of work space 
*
      NBIN = 60 
      LOW = 0.
      UPP = 2.
*
* Initialize 1D histogram with ID 10 and name 'One dim...'
*
      CALL HBOOK1(10,'Random + Random (0,1)',NBIN,LOW,UPP,0.)
*
* Loop to fill histograms entries
*
      DO I=1,100000
*
* Compute the sum of two random variates unif(0,1)            
*
        SUM = RNDM(I) + RNDM(I+1)
*
* Put one entry to 1D histogram
*
        WEIGHT = 1.
CCC        CALL HFILL(10,SUM,0,WEIGHT)
        CALL HF1(10,SUM,WEIGHT)
      ENDDO
*
* Now the histogram outputs
*
      CALL HISTDO                         ! On program listing
CCC      CALL HRPUT(0,'triangle_hist','N')   ! On a file for PAW
*
      STOP
      END
EOF
#
echo Compile and link
#
f77 demo_triangle.f -o demo_triangle.exe \
`cernlib`
#
echo Execute
#
demo_triangle.exe > demo_triangle.listing
#
echo The histogram listing is in file demo_triangle.listing