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

  // Variable initializations:
  time_ = ;
end

model friel_1995__rate_constants(kappa_L1, kappa_P1, kappa_L2, kappa_P2, time_, Ca_i)

  // Assignment Rules:
  kappa_L1 := piecewise( 0.000005 , ( geq(time_, 0)) && (time_ < 40 ), 0.00002 );
  kappa_L2 := kappa_L2_0 + kappa_L2_1 / (1 + power(Kd_Ca / Ca_i, n));
  Kd_Ca := piecewise( 1000 , ( geq(time_, 0)) && (time_ < 80 ), 500 );

  // Variable initializations:
  kappa_P1 = 0.132;
  kappa_P2 = 3.78;
  kappa_L2_0 = 0.054;
  kappa_L2_1 = 2.4;
  n = 3;
  time_ = ;
  Ca_i = ;
end

model friel_1995__cytosolic_calcium(Ca_i, Ca_o, Ca_i_ss, gamma, kappa_L1, kappa_P1, kappa_L2, kappa_P2, Ca_s, time_)

  // Assignment Rules:
  Ca_i_ss := Ca_o / (1.0 + kappa_P1 / kappa_L1);

  // Rate Rules:
  Ca_i' = -(kappa_L1 + kappa_P1 + gamma * (kappa_L2 + kappa_P2)) * Ca_i + gamma * kappa_L2 * Ca_s + kappa_L1 * Ca_o;

  // Variable initializations:
  Ca_i = 75;
  Ca_o = 2000000.0;
  gamma = 0.24;
  kappa_L1 = ;
  kappa_P1 = ;
  kappa_L2 = ;
  kappa_P2 = ;
  Ca_s = ;
  time_ = ;
end

model friel_1995__subspace_calcium(Ca_s, Ca_s_ss, Ca_i, Ca_i_ss, kappa_L2, kappa_P2, time_)

  // Assignment Rules:
  Ca_s_ss := Ca_i_ss * (1.0 + kappa_P2 / kappa_L2);

  // Rate Rules:
  Ca_s' = (kappa_L2 + kappa_P2) * Ca_i - kappa_L2 * Ca_s;

  // Variable initializations:
  Ca_s = 5300.0;
  Ca_i = ;
  Ca_i_ss = ;
  kappa_L2 = ;
  kappa_P2 = ;
  time_ = ;
end

model *friel_1995____main()

  // Sub-modules, and any changes to those submodules:
  environment: friel_1995__environment(time_);
  rate_constants: friel_1995__rate_constants(kappa_L1, kappa_P1, kappa_L2, kappa_P2, time_, Ca_i);
  cytosolic_calcium: friel_1995__cytosolic_calcium(Ca_i, Ca_o, Ca_i_ss, gamma, kappa_L1, kappa_P1, kappa_L2, kappa_P2, Ca_s, time_);
  subspace_calcium: friel_1995__subspace_calcium(Ca_s, Ca_s_ss, Ca_i, Ca_i_ss, kappa_L2, kappa_P2, time_);
end