#include "mat.h"

int 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;
   eos::state st; cin >> st;
   cout << "state = " << st << endl;
   cout << "p = " << an_eos->p(st) << '\n';
   cout << "c^2 = " << an_eos->csqv(st) << '\n';
   cout << "T = " << an_eos->t(st) << '\n';
}