//=====================================================================
// File: MaxwellBloch.java
//
//=====================================================================

//=====================================================================
// MaxwellBloch Class
//
//=====================================================================
public class MaxwellBloch extends NL3System {


  double alpha1=0.1, alpha2=0.1, beta=0.3, lambda=22;
  //-----------------------------------------------
  // The functions f(x,y,z), g(x,y,z) and h(x,y,z)
  //-----------------------------------------------
  
  protected double f(double x, double y, double z, double t)
    { return( beta*(y-x)); }

  protected double g(double x, double y, double z, double t)
    { return( alpha1*(x*z-y)); }

  protected double h(double x, double y, double z, double t)
    { return( alpha2*(lambda+1-z-lambda*x*y)); }

  MaxwellBloch() {
    super();
  }

  public void setParams(double ualpha1, double ualpha2, double ubeta,
                        double ulambda) {
     alpha1=ualpha1; alpha2=ualpha2; beta=ubeta;
     lambda=ulambda;
  }

};