LISTING OF PROGRAM

1 PROGRAM

2 $DEFINE ORDINATE ARRAY AND COMPUTE ABSCISSAE

3 X=INTS(5500,17000)/10000

4 V=5*(1 - EXP(-(X-1)) )**2

5 $SET ENERGY AND COMPUTE LIMITS

6 E=1.0

7 P2=E-V

8 CONSTRAIN(X,V,P2:P2.GT.0)

9 XMIN=MIN(X)

10 XMAX=MAX(X)

11 P=SQRT(P2)

12 $EVALUATE NORMALIZATION INTEGRAL

13 DX=X(2)-X(1)

14 NORM=SUM(1/P)*DX

15 PRB=1/(P*NORM)

16 $COMPUTE AVERAGES

17 X1AV=SUM(PRB*X)*DX

18 X2AV=SUM(PRB*X**2)*DX

19 XM1AV=SUM(PRB/X)*DX

20 $PRINT RESULTS

21 PRINT XMIN XMAX

22 PRINT NORM

23 PRINT X1AV

24 PRINT X2AV

25 PRINT XM1AV

26 PTS=NOELS(X)

27 PRINT PTS

28 END

EXECUTION STARTED

XMIN = .6304 XMAX = 1.5927

NORM = 1.5616

X1AV = 1.1672

X2AV = 1.475

XM1AV = .94537

PTS = 9624