Project 1 Text
Project 1 Text
PROGRAM MAIN
INTEGER N, X
EXTERNAL SUB1
COMMON /GLOBALS/ N
X = 0
PRINT *, 'Enter number of repeats'
READ (*,*) N
CALL SUB1(X,SUB1)
END
SUBROUTINE SUB1(X,DUMSUB)
INTEGER N, X, Y
EXTERNAL DUMSUB
COMMON /GLOBALS/ N
DATA Y /0/
IF(X .LT. N)THEN
X = X + 1
Y = Y + 1
PRINT *, 'x = ', X, ', y = ', Y
CALL DUMSUB(X,DUMSUB)
END IF
END
2PROGRAM MAIN
INTEGER I, I_START, I_END, I_INC
REAL A(100)
I_START = 1
I_END = 100
I_INC = 1
DO I = I_START, I_END, I_INC
A(I) = 0.0E0
END DO
END
3C
C NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, (c) 1985
C
C FILE: first.f
C
C FIRST PROGRAMMING EXPERIMENT
C
DATA N/25/, H/1.0/, X/0.5/
F = SIN(X)
G = COS(X)
DO 2 I = 1,N
H = 0.25*H
D = SIN(X + H) - F
Q = D/H
E = ABS(G - Q)
PRINT *,H,D,Q,E
2 CONTINUE
STOP
END
4C
C NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, (c) 1985
C
C FILE: xsinx.f
C
C EXAMPLE OF PROGRAMMING F(X) = X - SIN(X) CAREFULLY (F)
C
DOUBLE PRECISION DX,DY,DZ
DX = 1.0D0/15.0D0
DY = DSIN(DX)
DZ = DX - DY
PRINT *,DX,DY,DZ
X = 16.0
DO 2 N = 1,52
X = 0.25*X
Y = F(X)
PRINT *,N,X,Y
2 CONTINUE
X = SIN(1.0)
Y = F(1.0)
PRINT *,X,Y
STOP
END
FUNCTION F(X)
IF(ABS(X) .GE. 1.9) THEN
F = X - SIN(X)
ELSE
T = X**3/6.0
F = T
DO 2 N = 1,9
T = -T*X*X/REAL((2*N+2)*(2*N+3))
F = F + T
CONTINUE
END IF
RETURN
END