Computer Program
REAL ALPHA, Ti, K1, K2, K3, K4, W, H, N, A, B
WRITE(*,*)'INPUT INTEGRAION LIMITS'
READ(*,*) A, B
WRITE(*,*)'INPUT NUMBER OF ITERATIONS N'
READ(*,*) N
WRITE(*,*)'ENTER INITIAL CONDITION Y(0)'
READ(*,*) ALPHA
H = (B - A) / N
W = ALPHA
DO I=0, N
WRITE(*,*)' FOR I =',I
Ti = A + H * I
K1 = H * (-1 * W + Ti + 1)
K2 = H * (-1 * (W + (K1/2)) + (Ti + H/2) + 1)
K3 = H * (-1 * (W + (K2/2)) + (Ti + H/2) + 1)
K4 = H * (-1 * (W + K3) + (Ti + H) + 1)
W = W + (K1 + 2 * K2 + 2 * K3 + K4) / 6
WRITE(*,*)' K1 =',K1
WRITE(*,*)' K2 =',K2
WRITE(*,*)' K3 =',K3
WRITE(*,*)' K4 =',K4
WRITE(*,*)' W =',W
ENDDO
END
docsity.com
