<!-- hamlet
<!-- copytar
examples/                                                                                              755    1107     770            0  5635343070   5712                                                                                                                                                                                                                                                                                                                                                                      examples/roundoff1.f                                                                                   644    1107     770         1116  5506142063  10044                                                                                                                                                                                                                                                                                                                                                                      C Illustrates roundoff error for REAL numbers

C Created:  Joe Zachary, November 15, 1992
C Modified:


      PROGRAM ROUND1

      IMPLICIT NONE
      INTEGER I, N
      REAL PRODUCT, RATIO

      PRINT *, 'Please enter a positive integer value for N'
      READ *, N

      PRODUCT = 1

      DO I = N, 2, -2
         RATIO = REAL(I) / REAL(I-1)
         PRODUCT = PRODUCT * SQRT(RATIO)
      ENDDO

      DO I = N, 2, -2
         RATIO = REAL(I-1) / REAL(I)
         PRODUCT = PRODUCT * SQRT(RATIO)
      ENDDO

      PRINT 30, PRODUCT

 30   FORMAT (F9.7)

      STOP
      END

      
 7  cs365  	   !  cs366 Az   N>  cs367      7  cs376 pO     cs431 n1    >  cs451 rt  0 6G  cs677  	  @ d>  cs506 L  P o+  cs507 d  ` 	d3  cs508 o-  p z  cs509 n1   
X  cs511 rt   n  cs531  	   D  cs537    C  cs539     D  cs541 am   8  cs542      cs546     ɥ  cs547 n2    d%  cs100 rt   !e  cs561      *  cs56examples/roundoff2.f                                                                                   644    1107     770         1134  5506142063  10045                                                                                                                                                                                                                                                                                                                                                                      C Illustrates roundoff error for REAL numbers

C Created:  Joe Zachary, November 15, 1992
C Modified:


      PROGRAM ROUND1

      IMPLICIT NONE
      INTEGER I, N
      DOUBLE PRECISION PRODUCT, RATIO

      PRINT *, 'Please enter a positive integer value for N'
      READ *, N

      PRODUCT = 1

      DO I = N, 2, -2
         RATIO = DBLE(I) / DBLE(I-1)
         PRODUCT = PRODUCT * SQRT(RATIO)
      ENDDO

      DO I = N, 2, -2
         RATIO = DBLE(I-1) / DBLE(I)
         PRODUCT = PRODUCT * SQRT(RATIO)
      ENDDO

      PRINT 30, PRODUCT

 30   FORMAT (F16.14)

      STOP
      END

      
5  	   !  cs366 Az   N>  cs367      7  cs376 pO     cs431 n1    >  cs451 rt  0 6G  cs677  	  @ d>  cs506 L  P o+  cs507 d  ` 	d3  cs508 o-  p z  cs509 n1   
X  cs511 rt   n  cs531  	   D  cs537    C  cs539     D  cs541 am   8  cs542      cs546     ɥ  cs547 n2    d%  cs100 rt   !e  cs561      *  cs56examples/csqrt.f                                                                                       644    1107     770          455  5552343675   7277                                                                                                                                                                                                                                                                                                                                                                      C Illustrates complex numbers

C Created:  Joe Zachary, November 15, 1992
C Modified:


      PROGRAM COMPLX

      IMPLICIT NONE
      REAL X
      COMPLEX  C

      PRINT *, 'Please enter a REAL number'
      READ *, X

      C = SQRT(CMPLX(X))
      
      PRINT *, C

      STOP
      END

      
      DO I = N, 2, -2
         RATIO = DBLE(I) / DBLE(I-1)
         PRODUCT = PRODUCT * SQRT(RATIO)
      ENDDO

      DO I = N, 2, -2
         RATIO = DBLE(I-1) / DBLE(I)
         PRODUCT = PRODUCT * SQRT(RATIOexamples/rsqrt.f                                                                                       644    1107     770          471  5552343540   7303                                                                                                                                                                                                                                                                                                                                                                      C Illustrates complex numbers

C Created:  Joe Zachary, November 15, 1992
C Modified:


      PROGRAM COMPLX

      IMPLICIT NONE
      REAL X
      REAL C

      PRINT *, 'Please enter a REAL number'
      READ *, X

      C = SQRT(X)
      
      PRINT 100, C

 100  FORMAT (F7.3)

      STOP
      END

      
 N, 2, -2
         RATIO = DBLE(I) / DBLE(I-1)
         PRODUCT = PRODUCT * SQRT(RATIO)
      ENDDO

      DO I = N, 2, -2
         RATIO = DBLE(I-1) / DBLE(I)
         PRODUCT = PRODUCT * SQRT(RATIOexamples/soccer.txt                                                                                    644    1107     770          515  5635343071  10000                                                                                                                                                                                                                                                                                                                                                                      1 2 4 1
2 1 1 1
2 4 3 1
3 1 1 0
3 3 4 1
4 2 2 0
5 1 7 0
5 1 8 0
6 1 5 0
6 2 7 1
8 2 6 1
8 2 7 1
9 2 10 2
10 3 12 1
11 3 9 0
11 1 10 1
11 2 12 0
12 6 9 1
13 1 14 0
13 1 15 1
13 3 16 2
15 3 14 1
15 2 16 2
16 0 14 0
17 4 18 0
17 2 19 1
19 2 18 0
19 3 20 0
20 2 17 0
20 4 18 0
21 1 22 0
21 0 23 0
22 1 23 0
22 1 24 1
23 1 24 0
24 2 21 1
ATIO = DBLE(I) / DBLE(I-1)
         PRODUCT = PRODUCT * SQRT(RATIO)
      ENDDO

      DO I = N, 2, -2
         RATIO = DBLE(I-1) / DBLE(I)
         PRODUCT = PRODUCT * SQRT(RATIO                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                