Gibbs energy minimization   |POLVER05_3
f(v1) = -6.089+ln(v1/xs)+v11 #
f(v2) = -17.164+ln(v2/xs)+2*v11 #
f(v3) = -34.054+ln(v3/xs)+2*v11+v13 #
f(v4) = -5.914+ln(v4/xs)+v12 #
f(v5) = -24.721+ln(v5/xs)+2*v12 #
f(v6) = -14.986+ln(v6/xs)+v11+v12 #
f(v7)=-24.1+ln(v7/xs)+v12+v13 #
f(v8)=-10.708+ln(v8/xs)+v13 #
f(v9)=-26.662+ln(v9/xs)+2*v13 #
f(v10)=-22.197+ln(v10/xs)+v11+v13 #
f(v11)=v1+2*v2+2*v3+v6+v10-2 #
f(v12)=v4+2*v5+v6+v7-1 #
f(v13)=v3+v7+v8+2*v9+v10-1 #
xs=v1+v2+v3+v4+v5+v6+v7+v8+v9+v10 #
v1(0)=0.05
v2(0)=0.2
v3(0)=0.8
v4(0)=0.001
v5(0)=0.5
v6(0)=0.0007
v7(0)=0.03
v8(0)=0.02
v9(0)=0.1
v10(0)=0.1
v11(0)=10
v12(0)=10
v13(0)=10