// Created by libAntimony v2.8.0
model *Tham2008_PDmodel_TumourShrinkage()
// Compartments and Species:
compartment COMpartment_;
species $Ce in COMpartment_;
// Assignment Rules:
Exposure := piecewise(Dose, time < Cycle_Int*N_Cycle && Dose_Int1 < rem_time && rem_time < Dose_Length, Dose, time < Cycle_Int*N_Cycle && Dose_Int2 < rem_time && rem_time < Dose_Int2 + Dose_Length, 0);
Keq := ln(2)/Teq;
rem_time := (time*conversion_factor - floor((time*conversion_factor)/(Cycle_Int*conversion_factor))*Cycle_Int*conversion_factor)/conversion_factor;
RateIn := Size_0*Kover;
Effect := 1 - Ce/(AE50 + Ce);
Kover := 1/T_Turnover;
// Rate Rules:
Ce' = Exposure/1 - Ce*Keq;
Size' = (RateIn*Effect - Kover*Size)*Size;
// Species initializations:
Ce = 0;
// Compartment initializations:
COMpartment_ = 1;
Exposure has mg;
Keq has per_week;
rem_time has week;
// Variable initializations:
conversion_factor = 604800;
conversion_factor has sec_per_week;
Cycle_Int = 3;
Cycle_Int has week;
Dose = 5203.84;
Dose has mg;
N_Cycle = 6;
N_Cycle has dimensionless;
Dose_Int1 = 0;
Dose_Int1 has week;
Dose_Length = 0.44359;
Dose_Length has week;
Dose_Int2 = 1;
Dose_Int2 has week;
Size = 6.66;
Size has cm;
RateIn has per_week;
Effect has dimensionless;
Kover has per_cm_per_week;
AE50 = 10600;
AE50 has mg;
Teq = 7.67;
Teq has week;
Size_0 = 6.66;
Size_0 has cm;
T_Turnover = 21.8;
T_Turnover has cm_week;
// Other declarations:
var Exposure, Keq, rem_time, Size, RateIn, Effect, Kover;
const COMpartment_, conversion_factor, Cycle_Int, Dose, N_Cycle, Dose_Int1;
const Dose_Length, Dose_Int2, AE50, Teq, Size_0, T_Turnover;
// Unit definitions:
unit effect_unit = dimensionless;
unit week = 604800 second;
unit sec_per_week = second / 604800 second;
unit cm = 1e-2 metre;
unit mg = 1e-3 gram;
unit cm_week = 1e-2 metre * 604800 second;
unit per_week = 1 / 604800 second;
unit per_cm_per_week = 1 / (1e-2 metre * 604800 second);
unit m2 = metre^2;
// Display Names:
effect_unit is "normalised";
end