Due to a do... for cycle I obtain thes data:

....

> for i from 0 to 10 do x1:=0.1*i: T:=Tx1(x1);print('x1'=x1);y1:=eval((p1s*x1*gamma1)/p); print('y1'=y1);print('T'=T);end do:
                            x1 = 0.
                            y1 = 0.
                        T = 100.0000416
                            x1 = 0.1
                       y1 = 0.4454088320
                        T = 85.93510754
                            x1 = 0.2
                       y1 = 0.4597014971
                        T = 85.41761562
                            x1 = 0.3
                       y1 = 0.4343600675
                        T = 85.92005198
                            x1 = 0.4
                       y1 = 0.4122667776
                        T = 86.09782822
                            x1 = 0.5
                       y1 = 0.4080528845
                        T = 86.09991663
                            x1 = 0.6
                       y1 = 0.4306891136
                        T = 86.42991221
                            x1 = 0.7
                       y1 = 0.4904070733
                        T = 87.60592371
                            x1 = 0.8
                       y1 = 0.6013800104
                        T = 89.96757730
                            x1 = 0.9
                       y1 = 0.7762669037
                        T = 93.42524099
                            x1 = 1.0
                        y1 = 1.000000002
                        T = 97.13963466
                       y1 = 0.7762669038
                        T = 93.42524099
                            x1 = 1.0
                        y1 = 1.000000002
                        T = 97.13963466

Do you know how to export these data to Excel?Cause my problem is to plot T(x1),T(y1),y1(x1) and it is impossible for Maple to define T:=->y1 because the equation is very complex. [I can plot T(x1) cause it is easier]. This is my exercise if you need:

the sistem I have to solve is the sequent:

system of                        p*y1=p1s(T)*x1*gamma1(x1)
                                        {
                                          p*(1-y1)=p2s(T)*(1-y1)gamma2(x1)

So I wrote:

> restart;
> print("1 - ?");
                            "1 - ?"                #unknow fluid
> print("2 - H2O");
                           "2 - H2O"           #water
> A1:=7.74439:                           #experimental data
> B1:=1437.686:
> C1:=198.463:
> A2:=8.07126:
> B2:=1730.63:
> C2:=233.426:
> p:=760:
> A12:=2.7549:
> A21:=0.925:
> p1s:=10^(A1+(-B1/(T+C1))):
> p2s:=10^(A2+(-B2/(T+C2))):
> gamma1:=exp((A12+2*(A21-A12)*(x1))*(1-x1)^2):
> gamma2:=exp((A21+2*(A12-A21)*(1-x1))*(x1)^2):
> Teb1:=evalf(-(B1/(log[10](p)-A1))-C1):print(Tbp1=Teb1);                                     #Boiling Point Temperature of 1, useless
                       Tbp1 = 97.1396347
> Teb2:=evalf(-(B2/(log[10](p)-A2))-C2):print(Tbp2=Teb2);                                     #Boiling Point Temperature of 2, useless
                       Tbp2 = 100.0000416
> gx:= solve(((p1s*x1*gamma1)/p)+((p2s*(1-x1)*gamma2)/p)-1,T):                     #solution of the system
> Tx1:=x1->evalf(RootOf(10^(1/100000*(774439000*_Z+9928887257)/(1000*_Z+198463))*x1*exp(-1/10000*(-27549+36598*x1)*(-1+x1)^2)+10^(1/50000*(201781500*_Z+3835298419)/(500*_Z+116713))*exp(-1/5000*(-22924+18299*x1)*x1^2)-10^(1/50000*(201781500*_Z+3835298419)/(500*_Z+116713))*exp(-1/5000*(-22924+18299*x1)*x1^2)*x1-760)):
> hy:=solve(p1s*x1*gamma1 - p*y1,x1):
> #gy:=solve(p*(1-y1) + p2s*(h-1)*exp((A21+2*(A12-A21)*(1-h))*(h)^2),T):         # this is the complex equation ( I put # )

> Ty1:=y1->evalf(what I'd find from the gy):
> plot([Tx1,Ty1] ,0..1);

Sorry for the very complex reading but I'm an Italian Chemical Engineer student