Задача с условиями и интегралами на Fortran

program F
    real :: A(1:5), A1(1:5), B(1:5), B1(1:5), T(1:5), T1(1:5),H(1:5),D(1:5),V(1:5),X(1:5),CC(1:5),CC1(1:5),T,G,M,L,Y,O1,K1,K2,F1,R,C,a,b,S,S1,e
    integer :: i,n
    open(1, file = 'A.txt')
    open(2, file = 'B.txt')
    open(3, file = 'T.txt')
    open(4, file = 'predeli.txt')
    T = 2
    M = 3
    Y = 5
    L = (Y**2)/9,81
    e = 2,718
    
    !Изначальные условия
do i = 1,5
     
        read (1,*) A
        read (2,*) B
        read (3,*) T
        read (4,*) X
    
        if (A(i).gt.0 .and. B(i).gt.0) go to 111
        else
            A1(i) = 0
            B1(i) = 0
            T1(i) = 0
        endif
        111 continue    
        if ((A(i)/B(i)).le.T(i)) then
            A1(i) = A(i)
            B1(i) = B(i)
            T1(i) = T(i)
        else
            T1(i) = T(i)
            B1(i) = B(i)
            A1(i) = B1(i) * T1(i)
        endif
        
        !V
        n=size(X)
        S=0.0
        do i=1, n-1
        S1=S
        S=S+0.5*(X(i+1)-X(i))*(A1(i+1)+A1(i))
        V(i)=S1
        end do
        V(n)=S
        
        !H
        n=size(X)
        CC(i) = (1/12.) * (B1(i))**3
        S=0.0
        do i=1, n-1
        S1=S
        S=S+0.5*(X(i+1)-X(i))*(CC(i+1)+CC(i))
        H(i)=S1
        end do
        H(n)=S * (1/V(n))
        
        !D
        n=size(X)
        CC1(i) = ((-T1(i))/2.) * A1(i)
        S=0.0
        do i=1, n-1
        S1=S
        S=S+0.5*(X(i+1)-X(i))*(CC1(i+1)+CC1(i))
        D(i)=S1
        end do
        D(n)=S * (1/V(n)) + T
        
        !G и O1
        G = D(n) + H(n) + M
        O1 = G - T
        
        !K1
        do i = X
        K1(i) = ((sin(L *(B1(i)/2.)))/(L*B1(i))/2.) * (((1+L*T1(i)) * e ** (-L*T1(i) - 1)) / ((L**2) * T1(i)))
        end do
        
        !K2
        do i = X
        K2(i) = (-((e ** (-L*T1(i) - 1)) / ((L**2) * T1(i))))* (cos(L*(B1(i)/2.)) - ((sin(L*(B1(i)/2.)))/(L*(B1(i)/2.))))
        end do
        
        !F1
        do i = X
        F1(i) = -((1 - (e**(-L*T1(i))))/(L*T1(i))) * ((sin(L*(B1(i)/2.)))/((L*B1(i))/2.))
        end do
        
        !C
        do i = X
        if (A1(i) .eq. 0 .and. B1(i) .eq. 0) then
        C(i) = 0
        else
            C(i) = A1(i) * (K1(i) + K2(i) + F1(i) * O1)
        endif
        end do
        
        !R
        do i = X
        n=size(X)
        S=0.0
        do i=1, n-1
        S1=S
        S=S+0.5*(X(i+1)-X(i))*(C(i+1)+C(i))
        R(i)=S1
        end do
        R(n)=abs(S / (V(n) * M))
        end do
        print *, R
end do
end program F

program F
    Implicit none
 
    integer ::  k
    real,parameter :: TT = 1.00, M = 0.23
 
    real :: O1, V
    real,dimension(21) :: A, B, T, X, K1, K2, F1, C, KKF
    real,dimension(25) :: Y, L
 
    open(1, file = 'A.txt')
    read(1,*) A
    close(1)
 
    open(2, file = 'B.txt')
    read(2, *) B
    close(2)
 
    open(3, file = 'T.txt')
    read(3, *) T
    close(3)
 
    open(4, file = 'predeli.txt')
    read(4, *) X
    close(4)
 
    where(A>0 .and. B>0)
        where(A/B <= T)
            T=A/B
        elsewhere
            A=B*T
        end where
    end where
    !Для нахождения V
    V = integrate(X, A)
 
    !Найденные H, D, G
    !Решение для 3-й картинки
    O1 = V/2.0 * (integrate(X,-T*A)  + integrate(X, B**3)* 1.0/6.0 ) - M 
 
    open(5, file = 'Y.txt')
    read (5,*) Y
    close(5)
 
    L = Y**2.0/9.81
 
  do k = 1,25
 
        K1 = (sin((L(k)*B)/2.0))/((L(k)*B)/2.0)
        K2 = (exp(-L(k)*T - 1.0) / (L(k)**2 * T)) * ((1 + L(k) * T) - (1.0/(exp(-L(k)*T))))
        F1 = ((exp(-L(k)*T)) * cos(L(k)*B/2.0)) / ((L(k)) **2) * T 
        KKF = K1 * K2 - F1
 
    !C
    C=0
    where(A/=0 .and. B/=0) C = A * (KKF * O1)
 
    !Ответ
 
   print *, Y(k), abs(integrate(X, C)/V*M)
 
  end do
 
contains
 
  function integrate(X,Y) result(S)
      real, dimension(1:), intent(in):: X,Y
    real :: S
    integer:: n
    n=size(X)
        S= dot_product(X(2:n)-X(1:n-1),(Y(1:n-1))/2)
  end function integrate
 
end program F

program F
 
    real,dimension(21) :: A, B, T, X, K1, K2, C, F1, p, pl, pll
    real,dimension(25) :: Y, L
 
    character*60 c1,c2,c3,c4,c5,c6
      c1='d:\Code\A.txt'
      c2='d:\Code\B.txt'
      c3='d:\Code\T.txt'
      c4='d:\Code\predeli.txt'
      c5='d:\Code\Y.txt'
      c6='d:\Code\redcoeff.txt'
      open(unit=1,file=c1,status='old')
      open(unit=2,file=c2,status='old')
      open(unit=3,file=c3,status='old')
      open(unit=4,file=c4,status='old')
      open(unit=5,file=c5,status='old')
      open(unit=6,file=c6,status='unknown')
      read(1,*) A
      close(1)
      read(2, *) B
      close(2)
      read(3, *) T
      close(3)
      read(4, *) X
      close(4)
      read (5,*) Y        
      close(5)
     
    TT = 1.00  
    GM=0.23     
    T=A/B          
    V = integrate(X,A)    
 
    
    AKB=TT+integrate(X,-T/2.*A)/V 
    BM=integrate(X,(B**3)/12.)/V
    AKG =AKB+BM-GM     
 
    OG=AKG-TT
    L=Y**2./9.81     
    write(6,*) '      i          ans'
  do k = 1,25
     P=L(k)*B/2
     PL=L(k)*T
     PLL=T*(l(k))**2
     K1 = sin(P)/P*((1.+PL)*exp(-PL)-1.)/PLL
     K2 = -exp(-PL)/PLL*(cos(P)-sin(P)/P)
     F1 = -(1.-exp(-PL))/PL*sin(P)/P
     C=0.
     where(A/=0.and.B/=0) C=A*(K1+K2+F1*OG)
     r=abs(integrate(X,C)/(V*GM))     
     print *, Y(k),r
     write(6,*) Y(k),r
     end do
     close(6)
    contains
 
  function integrate(X,Y) result(S)
    real, dimension(1:), intent(in):: X,Y
    n=size(X)
    S= dot_product(X(2:n)-X(1:n-1),(Y(2:n)-Y(1:n-1)))/2)
    end do
    end function integrate
 
    end program F