#include <mat.h>

void main()
{
   /*
   int neos;
   cin >> neos;
   eos eoslist[neos];
   int ieos; for (ieos=0;ieos<neos;ieos++) eoslist[ieos] = readeos(cin);

   int ncalcs;
   cin >> ncalcs;
   int icalc; for (icalc=0;icalc<ncalcs;icalc++) {
      int mat; cin >> mat;
      double d,e; cin >> d >> e;
      cout << "p = " << eoslist[mat].p(d,e) << '\n';
   }
   */
   eos *an_eos;
   cin >> an_eos;
   cout << an_eos << endl;
   Array<eos::state> states; cin >> states;
   for (int i=0;i<states.n;i++) {
      cout << "state = " << states[i] << '\n';
      cout << "p = " << an_eos->p(states[i]) << '\n';
      cout << "c^2 = " << an_eos->csqv(states[i]) << '\n';
      cout << "T = " << an_eos->t(states[i]) << '\n';
   }
}