“…error of the calculated vapor pressure c c**********************************************************************c parameter(nn=3,mm=100,ml=100) implicit double precision(a-h,o-z) integer maxfev,mode,nprint,nfev,ldfjac, 1 ipvt(nn) double precision double precision 1 ftol,xtol,gtol,epsfcn,factor x(nn),f(mm),diag(nn),fjac(mm,nn),qtf(nn), wal(nn),wa2(nn),wa3(nn),wa4(mm),diag(nn), pp(ml),x2(ml),y2(ml),y3(ml),pcal(ml),y2cal(ml), gamaA (ml), gamaB (ml), x3 (ml), PA0,PBO,k1,K2,beta,vA,vB,T 2 3 4 double precision dpmpar character*100 infile,infile2,outfile,infofile character*75 name external func common/plock1/pp,x2,y2 common/block2/pcal,y2cal,gamaA,gamaB common/block3/T,pa0,pb0 common/block4/kl,K2,vA,vB,beta common/block7/infile2 open(unit=9,file='run.ctl' ,status='unknown') read(9,*)infile read(9,*)infile2 read(9,*)outfile read(9,*)infofile close (9) **************************************************************** number of parameters to be optimized n=nn total nuber of functions to be minimized m=l convergenz and toleranz parameters of LMDIF ftol=dsqrt(dpmpar(1)) xtol=dsqrt(dpmpar (1) func(l,n,x,f,iflag) c***~***************************************************************c c This subroutine calculates the residuals of the function c c to be minimized c c*******************************************************************c …”