【正文】 十三：附錄 include include double a[107]={,,,,20,25,,40,54,,}。double f(double x){int i,j,k。for(k=0。k++){ if(x=a[k+1]amp。x=a[k]) break。i2。 for(j=1。j++) { if(j!=i) m=m*(xa[j+k])/(a[i+k]a[j+k])。 } return y。 for(x=30。x+=) { a=(f()xd)/()。 c=(f(x+)xd)/(x+)。amp。 printf(%f %f回流比=%f\n,x,b,r)。 double xW=,xD=,xF=。 y=R/(R+1)*x+(xDxW)*x/(R+1)/(xFxW)(xDxF)*xW/(R+1)/(xFxW)。}double jld(double x,double R) /*精餾段操作線方程*/{ double y。 y=R/(R+1)*x+(R+1)。}double lglr(double y,int N){ double xy[107][2]={{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{11,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{20,},{,},{,},{,},{,},{,},{,},{25,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{40,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{54,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,},{,78},{,},{,},{,},{,},{,},{,},{,},{,}}。 double x。i107。amp。 x= (yxy[i1][1])*(yxy[ i ][1])*(yxy[i+1][1]) / ( (xy[i2][1]xy[i1][1])*(xy[i2][1]xy[ i ][1])*(xy[i2][1]xy[i+1][1]) ) *xy[i2][0] + (yxy[i2][1])*(yxy[ i ][1])*(yxy[i+1][1]) / ( (xy[i1][1]xy[i2][1])*(xy[i1][1]xy[ i ][1])*(xy[i1][1]xy[i+1][1]) ) *xy[i1][0] + (yxy[i2][1])*(yxy[i1][1])*(yxy[i+1][1]) / ( (xy[ i ][1]xy[i2][1])*(xy[ i ][1]xy[i1][1])*(xy[ i ][1]xy[i+1][1]) ) *xy[ i ][0] + (yxy[i2][1])*(yxy[i1][1])*(yxy[ i ][1]) / ( (xy[i+1][1]xy[i2][1])*(xy[i+1][1]xy[i1][1])*(xy[i+1][1]xy[ i ][1]) ) *xy[i+1][0]。 return (x)