#include <stdio.h>
#include <math.h>
#include <stdlib.h>
main()
{ 
  double T,dT,xg,xl, T1, T2, L1, L2;
  double A,B;
   FILE *archivo;
   FILE *archivo2;
  printf("\n Este programa calcula y dibuja las líneas de separacion");
  printf("\n  de fases para una disolución ideal");
  printf("\n dando ls temperaturas de ebullicion y calores latentes de los");
  printf("\n dos componentes puros");
  printf("\nGuarda los resultados en el archivo de texto Disol_ideal.dat");
  printf("\n Para que dibuje es necesario tener el programa libre wgnuplot.exe");
  printf("\n en la misma carpeta que este ejecutable");
  printf("\n o bien añadir a la variable de entorno PATH"); 
  printf("\nla carpeta donde se encuentre wgnuplot.exe");
  printf("\n el archivo Disol_ideal.gnp es para uso del programa wgnuplot.exe");
 printf("\n\n\nTemperaturas de ebullición T1, T2 (dar T1<T2 en K): ");
 scanf("%lf %lf",&T1,&T2);
 printf("\nCalores latentes L1/R, L2/R (K): ");
 scanf("%lf %lf",&L1,&L2);
 archivo=fopen("Disol_ideal.dat","w"); 
 fprintf(archivo,"#Diagrama T/x de una disolucion ideal");
 fprintf(archivo,"\n#T1 = %.2lf K, L1 = %lf J/mol",T1,L1);
 fprintf(archivo,"\n#T2 = %.2lf K, L2 = %lf J/mol",T2,L2);
 fprintf(archivo,"\n#    T     xl      xg   ");
 fprintf(archivo,"\n#=======================");
 dT = (T2-T1)/1000.; // Toma 1001 puntos entre T1 y T2
 for(T=T1;T<T2+dT;T+=dT)
 {A=exp(L1*(T-T1)/T/T);
  B=exp(L2*(T-T2)/T/T);
  xl=(1-B)/(A-B);
  xg=A*xl;
 fprintf(archivo,"\n%.2lf %lf %lf",T,xl,xg);
 }
 fclose(archivo);
 archivo2=fopen("Disol_ideal.gnp","w");
 fprintf(archivo2,"plot 'Disol_ideal.dat' using 2:1 with lines"); 
 fprintf(archivo2,"\nreplot 'Disol_ideal.dat' using 3:1 with lines");
 fclose(archivo2);
 system("wgnuplot Disol_ideal.gnp -");
      
}