//Created by libAntimony v2.4
model shrestha_2010__environment(time_)

  // Variable initializations:
  time_ = ;
end

model shrestha_2010__x1(x1, k, lambda_Ca, lambda1, time_)

  // Rate Rules:
  x1' = k - (lambda_Ca * x1 + lambda1 * x1);

  // Variable initializations:
  x1 = 0.00;
  k = ;
  lambda_Ca = ;
  lambda1 = ;
  time_ = ;
end

model shrestha_2010__lambda_Ca(lambda_Ca, A, B, Ca, S, m_Ca)

  // Assignment Rules:
  lambda_Ca := (A - B) / (1 + power(Ca / S, m_Ca)) + B;

  // Variable initializations:
  A = ;
  B = ;
  Ca = ;
  S = ;
  m_Ca = ;
end

model shrestha_2010__m_Ca(m_Ca, m1, m2, beta, Ca, R)

  // Assignment Rules:
  m_Ca := m1 / (1 + exp(- beta * (R - Ca))) + m2;

  // Variable initializations:
  m1 = ;
  m2 = ;
  beta = ;
  Ca = ;
  R = ;
end

model shrestha_2010__x2(x2, x1, lambda_Ca, lambda2, time_)

  // Assignment Rules:
  PTH := x2 / 2.75;

  // Rate Rules:
  x2' = lambda_Ca * x1 - lambda2 * x2;

  // Variable initializations:
  x2 = 0.00;
  x1 = ;
  lambda_Ca = ;
  lambda2 = ;
  time_ = ;
end

model shrestha_2010__Ca(Ca, Ca_0, Ca_1, alpha, t0, time_)

  // Assignment Rules:
  Ca := piecewise( Ca_0 , time_ < t0 , Ca_0 - Ca_1 * (1 - exp(- alpha * (time_ - t0))) );

  // Variable initializations:
  Ca_0 = ;
  Ca_1 = ;
  alpha = ;
  t0 = ;
  time_ = ;
end

model shrestha_2010__k(k, lambda1, lambda2, x1_n, x2_n)

  // Assignment Rules:
  k := lambda2 * x2_n + lambda1 * x1_n;

  // Variable initializations:
  lambda1 = ;
  lambda2 = ;
  x1_n = ;
  x2_n = ;
end

model shrestha_2010__S(S, Ca_0, A, B, x1_n, x2_n, lambda2, m_Ca)

  // Assignment Rules:
  S := Ca_0 * power(-((x1_n * B - lambda2 * x2_n) / (x1_n * A - lambda2 * x2_n)), 1 / m_Ca);

  // Variable initializations:
  Ca_0 = ;
  A = ;
  B = ;
  x1_n = ;
  x2_n = ;
  lambda2 = ;
  m_Ca = ;
end

model shrestha_2010__A(A, lambda1, lambda2, x2_max, k)

  // Assignment Rules:
  A := (lambda1 * lambda2 * x2_max) / (k - lambda2 * x2_max);

  // Variable initializations:
  lambda1 = ;
  lambda2 = ;
  x2_max = ;
  k = ;
end

model shrestha_2010__B(B, lambda1, lambda2, x2_min, k)

  // Assignment Rules:
  B := (lambda1 * lambda2 * x2_min) / (k - lambda2 * x2_min);

  // Variable initializations:
  lambda1 = ;
  lambda2 = ;
  x2_min = ;
  k = ;
end

model shrestha_2010__model_parameters(lambda1, lambda2, m1, m2, R, beta, x1_n, x2_n, x2_min, x2_max, t0, alpha, Ca_0, Ca_1)

  // Variable initializations:
  lambda1 = 0.0125;
  lambda2 = 0.5595;
  m1 = 112.5200;
  m2 = 15.00;
  R = 1.2162;
  beta = 1e6;
  x1_n = 490.7800;
  x2_n = 6.6290;
  x2_min = 0.6697;
  x2_max = 14.0430;
  t0 = 575.0;
  alpha = 0.0442;
  Ca_0 = 1.255;
  Ca_1 = 0.1817;
end

model *shrestha_2010____main()

  // Sub-modules, and any changes to those submodules:
  environment: shrestha_2010__environment(time_);
  x1: shrestha_2010__x1(x10, k0, lambda_Ca0, lambda1, time_);
  lambda_Ca: shrestha_2010__lambda_Ca(lambda_Ca0, A0, B0, Ca0, S0, m_Ca0);
  m_Ca: shrestha_2010__m_Ca(m_Ca0, m1, m2, beta, Ca0, R);
  x2: shrestha_2010__x2(x20, x10, lambda_Ca0, lambda2, time_);
  Ca: shrestha_2010__Ca(Ca0, Ca_0, Ca_1, alpha, t0, time_);
  k: shrestha_2010__k(k0, lambda1, lambda2, x1_n, x2_n);
  S: shrestha_2010__S(S0, Ca_0, A0, B0, x1_n, x2_n, lambda2, m_Ca0);
  A: shrestha_2010__A(A0, lambda1, lambda2, x2_max, k0);
  B: shrestha_2010__B(B0, lambda1, lambda2, x2_min, k0);
  model_parameters: shrestha_2010__model_parameters(lambda1, lambda2, m1, m2, R, beta, x1_n, x2_n, x2_min, x2_max, t0, alpha, Ca_0, Ca_1);
end