BermanLambdaPhase.h

#ifndef BERMANLAMBDAPHASE_H
#define BERMANLAMBDAPHASE_H

#include "BermanPhase.h"

class BermanLambdaPhase: public Phase {
 public:
  BermanLambdaPhase();
  virtual ~BermanLambdaPhase();

  // setting methods ////////////////////////////////////
  virtual void setTk(double ltk);
  virtual void setPa(double lpa);
  virtual void resetTPbounds(double lowTbound,double hightTBound,
                             double lowPbound,double highPbond);

  // getting methods //////////////////////////////////
  virtual double getGibbs();  
  virtual double getEnthalpy();
  virtual double getEntropy();
  virtual double getCp();
  virtual double getVolume();
  virtual double getdVdT();
  virtual double getdVdP();

 protected:
  BermanPhase *alpha, *beta;
  double tt0;
  double l1,l2;
  double ttpk;
  double ttr0;
  virtual void init(BermanData* alph,BermanData* bet);
 private:
  const double pbr;
  double pb;
  double tt;
  double ttr;
  double ttd;
  double x1, x2, x3, x4,
    Dx1Dp, D2x1Dp2, D3x1Dp3, Dx2Dp, D2x2Dp2, D3x2Dp3,
    Dx3Dp, D2x3Dp2, D3x3Dp3, Dx4Dp, D2x4Dp2, D3x4Dp3;

  double G(double tk,double pb);
  double H(double tk,double pb);
  double S(double tk,double pb);
  double Cp(double tk,double pb);
  double V(double tk,double pb);
  double dVdT(double tk,double pb);
  double dVdP(double tk,double pb);
  void filx(void);
};

#endif