PROGRAM MAIN
        INTEGER I
        REAL R,V,m,mmasa,mnozina,Vm,Tc,Pc,W,a1,a2,a3,a4,b,c1,c2
        Open (25,File='tlaki.TXT', STATUS='REPLACE')
        R=8.314
        V=1.5
        m=4000
        mmasa= 28
        mnozina=m/mmasa
        Vm=2.5*1000.0/mnozina
        Tc=126
        Pc=3400
        W=0.04

        CALL spremenljivke (R,Tc,Pc,W,a1,a2,a3,a4,b,c1,c2)
        Write (25,*) 'Tlak butana po Lee-Edmisterjevi in splosni plinsk'
        write (25,*) 't/oC', 'T/K', 'p/bar', 'pi/bar', 'p razlika'
!100     format (1X, A4, 6X, A3, 7X, A5, 6X, A6, 5X, A9)
         print*, "Gowno"
        do I=1,16,1
        Tk=563+I*10
        Tcelz=Tk-273
        a=a1-a2*Tk+a3*Tk**(-1.0)+a4*Tk**(-5.0)
        c=c1*Tk**(-0.5)+c2*Tk**(-2.0)
        p1=R*Tk/(Vm-b)-a/(Vm*(Vm-b))+c/(Vm*(Vm-b)*(Vm+b))
        p2=mnozina*R*Tk/(2.5*1000)
        praz=p2-p1
        write (25,*) Tcelz, Tk, p1, p2, praz

        end do
!200     format (1X, F5.1, 5X, F5.1, 5X, F6.3, 5X, F6.3, 5X, F5.3)
        close(25)
        PAUSE

        end  PROGRAM MAIN


        SUBROUTINE spremenljivke (R,Tc,Pc,W,a1,a2,a3,a4,b,c1,c2)
        a1=R**2.0*Tc**2.0/Pc*(0.25913-0.031314*W)
        a2=R**2.0*Tc/Pc*(0.0249+0.15369*W)
        a3=R**2.0*Tc**3.0/Pc*(0.2015+0.21642*W)
        a4=R**2.0*Tc**7.0/Pc*0.042*W
        b=R*Tc/Pc*0.0982
        c1=R**3.0*Tc**3.5/Pc**2.0*(0.059904)*(1.0-W)
        c2=R**3.0*Tc**5.0/Pc**2.0*(0.018126+0.091994*W)


        END SUBROUTINE spremenljivke