程序已经写好了。我不是学编程的，但设计要用到这个程序。求助各位编程高手们帮我解决它。感激不尽啊。有问题可以直接问我，544701808。文件怎么上传不上呢？
贴在这吧。     real x(3,1),sss(5,5),qk(5),rk(5),t0(1),y(3,1)
         real P1(3,1),Yi(3,1),yir(3),yrc(3),k(3,1),t
         real r(3),l(3),nxx,nix,pp,arph(3,1),xkk(3,1)
         real nx(3),xnix(3,5),xx(5),XW(3,1),add,newxt(3,1),newxj(3,1)
         real q(3),sj(5),ff(5,5),f(5),xxi(3,5),xsj(3,5)
         real fi(3,5),bbbb(3,5),dddd(3),yy(3,1)
         real x1,eeee(3),xsss(3),xd(3,1),xm0,xm1,xm2,ph(3,1)
         integer nn(3,5)
         WRITE（*，*）‘INPUT X’
         READ(*，*)X
         write(*,*)‘input t0(1)’
         read(*,*) t0(1)
         m=5
         n=3
         ji=1
         nnn=1
         er1=0.001
         p0=760
         mm=5
          data sss/0.0,-146.8,340.7,62.32,362.3,167.0,0.0,4102,268.2,377.6,
    -144.3,-269.7,0.0,349.2,692.7,537.4,603.8,730.4,0.0,-14.0,903.8,5693.0,634.2,-66.17,0.0/
         data xd/0.003,0.007,0.99/
        data xw/0.0068,0.93,0.002/
         data rk/0.5313,1.2663,1.1167,1.3013,0.92/
         data qk/0.4,0.986,0.867,1.224,1.40/
         data nn/4,0,0,2,0,0,0,1,0,0,2,0,0,0,1/
    35    t=t0(1)
    c    find yic activity coefficiention
         rx=0
         jjj=ji
         t00=t0(ji)+273
         do 20 kk=1,5
         r(j)=r(j)+nn(j,kk)*rk(kk)
    20   continue
         rx=rx+r(j)+nn(j,kk)*rk(kk)
    10   continue
         qx=0
         do 30  j=1,3
         q(j)=0
         do 40 kk=1,5
         q(j)=q(j)+nn(j,kk)*qk(kk)
    40   continue
         qx=qx+q(j)*l(j)
    30   continue
        do 50 j=1,3
        l(j)=5*(r(j)-q(j))-r(j)+1
    50  continue
        x1=0
        do 60 j=1,3
          x1=x1+x(j,1)*l(j)
    60  continue
        do 111  i=1,3
        zz=r(i)/rx
        si=qi/qx
        yic=log(zz)+5.0*q(i)*log(si/zz)+l(i)-zz*xl
    111  continue
    c    find yir activity coefficiention
    c    find yir/fk activity coefficent
        do 110 jj=1,5
        nxx=0
        do 70 j=1,3
       nx(j)=0
        do 80 kk=1,5
         nx(j)=nx(j)+nn(j,kk)*x(j,1)
    80  continue
    nxx=nxx+nx(j)
    70   continue
        qx=0
        do 90  kk=1,5
        nix=0
       do 100 j=1,3
       xnix(j,kk)=nn(j,kk)*x(j,1)
       nix=nix+xnix(j,kk)
    100  continue
       xx(kk)=nix/nxx
        qx=qx+qk(kk)*xx(kk)
    90  continue
        qjx=qk(jj)*xx(jj)
       sj(jj)=qjx/qx
    110  continue
        do 200  kk=1,5
        do 210 jj=1,5
        ff(jj,kk)=exp(-sss(jj,kk)/t00
    210  continue
    200  continue
        do 300 kk=1,5
        siff=0.0
        do 120 n2=1,5
        sjf=sj(n2)*ff(KK,n2)
        ssjf=0.0
        do 130 n1=1,5
        sssjf=ssjf+sjf/ssjf
    120  continue
        ssss=0.0
        do 140 n3=1,5
        ssss=ssss+sj(n3)*ff(n3,kk)
    140  continue
        f(kk)+qk(kk)*(1.0-log(ssss)-sjff)
    300  continue
    c   find yir/fki activity cofficent
        do 150
        dsss=0
        do 151 j=1,5
         xxi(n4,j)=xnix(n4,j)/nx(n4)
        dsss=dsss+qk(j)*xxi(n4,j)
    151  continue
        dddd(n4)=dsss
    150  continue
        Do 160  n4=1,3
        Do 160  jj=1,5
       aaaa=qk(jj)*xxi(n4,jj)
       xsj(n4,jj)=aaaa/dddd(n4)
    160  continue
       Do 190 n4=1,3
       Do 190 jj=1,5
       Bbbb(n4,jj)=0
       Do 191 n6=1,5
       Bbbb(n4,jj)=bbbb(n4,jj)+xsj(n4,n6)*ff(n6,jj)
    191  continue
    190  continue
       Do 400 n4=1,3
       Do 400 kk=1,5
       Xsss(n4)=0
       Do 170 jj=1,5
         Xsss(n4)=xsss(n4)+xsj(n4,jj)*ff(jj,kk)
    170  continue
       Eeee(n4)=0
       Do 180 jj=1,5
       Cccc=xsj(n4,jj)*ff(kk,jj)
       Eeee(n4)=eeee(n4)+cccc/bbbb(n4,jj)
    180  continue
       Fi(n4,kk)=qk(kk)*(1-log(xsss(n4))-eeee(n4))
    400  continue
    C    find yir activity cofficent result
         500 i=1,3
         Yir(i)=yir(i)+nn(i,kk)*(f(kk)-fi(i,kk))
    501  continue
    500  continue
    C   find yi activity cofficent result
        Do 600 i=1,3
    600  continue
        A=16.323-3366.99/(t00-59.46)
        P1(1,ji)=exp(a)
        B-18.306-3816.44/(t00-59.46)
        P1(3,ji)=exp(b)
        C=7.62029-2297.25/(t0(ji)=273.321)
        Cc=10*(c)
        P1(2,ji)=cc*(1000/133.33)
        Pp=p1(2,ji)*yi(2,ji)
        Do 999 i=1,3
        Arph(i,ji)=p1(i,1)*yi(i,ji)/pp
    999 continue
        Do 888 i=1,3
       Xkk(i,ji)=arph(i,ji)*x(i,ji)
    888 continue
        Add=arph(1,ji)*x(1,ji)+arph(2,ji)+arph(2,ji)*x(2,ji)+arph(3,ji)*x(3,ji)
        Do 777 j=1,3
        Yy(i,ji)=arph(i,ji)*x(i,ji)/add
    777 continue
        Data newxj/0,0,0,0,0,0/
        Data newxt/0,0,0,0,0,0/
        Do 900 j=1,3
        Newxj(j,1)=4.67*yy(j,1)-3.667*xd(j,1)
    900  continue
        Do 901 i=1,3
        Newxt(i,1)=0.22*xw(i,1)+0.78*yy(i,1)
    901  continue
        Do 444 i=1,3
        Ph(i,ji)=arph(i,ji)*x(i,ji)
    444  continue
        Xm0=0.0
        Xm1=0.0
        Xm2=0
        Do 333  i=1,3
       Xm0=xm0+newxj(i,ji)
       Xm1=xm1+newxt(i,ji)
        Xm2=xm2+yy(i,ji)
    333  continue
       Write(*,*)'x=',x
       Write(*,*)'t0=',t0
       Write(*,*)'yi=',yi
       Write(*,*)'p1=',p1
       Write(*,*)'pp=',pp
       Write(*,*)'arph=',arph
       Write(*,*)'ph=',ph
       Write(*,*)'yy=',yy
       Write(*,*)'xt=',newxt
       Write(*,*)'xj=',newxj
       Write(*,*)'xm0(newj)=',xm0
       Write(*,*)'xm1(newt)=',xm1
       Write(*,*)'zm2(yy)=',xm2
       end