program main
        implicit none
        real::R,V,VM,TC,PI(16),PC,A1,A2,A3,A4,B,C1,C2,CI(16),PI2(16)
        REAL::O,N,M,AI(16),TI(16),POV(16)
        integer::I
        OPEN(20,FILE="tlak.DAT")
        TC=126
        PC=3400
        O=0.04
        R=8.314
        V=1500
        M=4000
        N=M/28.0134
        VM=V/N
        WRITE(20,200) "TEMPERATURA:","TLAK(1):","TLAK(2):",'RAZLIKA'
        CALL ENACBE(R,TC,PC,A1,A2,A3,A4,B,C1,C2,O)
 200    FORMAT(A15,3X,A10,6X,A10,3X,A10)
        DO I=1,16,1
        TI(I)=563+I*10
        AI(I)=A1-A2*TI(I)+A3*TI(I)**(-1)+A4*TI(I)**(-5)
        CI(I)=C1*TI(I)**(-0.5)+C2*TI(I)**(-2)
        PI(I)=R*TI(I)/(VM-B)-AI(I)/(VM*(VM-B))+CI(I)/(VM*(VM-B)*(VM+B))
        PI(I)=PI(I)/100
        PI2(I)=N*R*TI(I)/150000
        POV(I)=PI(I)-PI2(I)
           !WRITE(6,100) TI(I), PI(I),PI2(I)
           WRITE(20,*) TI(I), PI(I),PI2(I),POV(I)
        ENDDO
        CLOSE (20)
        PAUSE
        END PROGRAM MAIN





        SUBROUTINE ENACBE(R,TC,PC,A1,A2,A3,A4,B,C1,C2,O)
                   A1=R**2*TC**2/PC*(0.25913-0.031314*O)
                   A2=R**2*TC/PC*(0.0249+0.15369*O)
                   A3=R**2*TC**3/PC*(0.2015+0.21642*O)
                   A4=R**2*TC**7/PC*(0.042*O)

                   B=R*TC/PC*0.0982

                   C1=R**3*TC**3.5/PC**2*(0.059904)*(1-O)
                   C2=R**3*TC**5/PC**2*(0.018126+0.091944*O)
        END SUBROUTINE ENACBE