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="REZULTATI.DAT")
       TC=126
       PC=3400
       O=0.04
       R=8.314
       V=1500
       M=4000
       N=M/28.0134
       VM=V/N
       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)
        WRITE(6,200) "TEMPERATURA:","TLAK(1):","TLAK(2):",'RAZLIKA'
        !WRITE(20,200) "TEMPERATURA:","TLAK(1):","TLAK(2):",'RAZLIKA'
       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,*) TI(I), PI(I),PI2(I),POV(I)
        !WRITE(20,*) TI(I), PI(I),PI2(I),POV(I)
       ENDDO
!100    FORMAT(5X,F10.0,1X,"K",15X,F7.4,1X,"BAR",10X,F7.4,1X,"BAR")
200    FORMAT(A15,3X,A10,6X,A10,3X,A10)

       PAUSE
       END PROGRAM MAIN