//-------------------------------------------------------------------------
//
// Family of curves of Te = epsilon (Kpm - ln(epsilon))
// for 0<epsilon<=1 and some fixed values of Kpm.
//
// Note that Te = - epsilon ln(epsilon/exp(Kpm)).
//
// Author: Jose Paulo V. S. da Cunha
//
// Rio de Janeiro, January 05, 2004.
//
//-------------------------------------------------------------------------

// epsilon range:

epsilon=0.002:0.002:1;

// Computes the family of curves:

Te = [epsilon.*(-log(epsilon));
      epsilon.*(-log(epsilon/exp(0.2)));
      epsilon.*(-log(epsilon/exp(0.4)));
      epsilon.*(-log(epsilon/exp(0.6)));
      epsilon.*(-log(epsilon/exp(0.8)));
      epsilon.*(-log(epsilon/exp(1.0)));
      epsilon.*(-log(epsilon/exp(1.2))) ]';

// Plots the family of curves:

xbasc();
xset("default");
xset("font size",6);
xset("line style",1);
xset("thickness",2);

plot2d(epsilon,Te,rect=[0,0,1,1.2],nax=[2,5,2,6]);

// Plot grids and comments:

xset("line style",2);
xset("thickness",1);
xset("mark size",0);
xgrid(1);

xtitle('','epsilon','Te');