# States are the number of zebrafish embryo in each states, they can be output States = {P_N_N, P_N_H, P_N_E, P_N_C, P_N_EH, P_E_C, P_E_H, P_E_EH, P_E_E, P_H_EH, P_H_H, P_EH_H, P_EH_EH}; Inputs = {q_E_H_on}; Outputs = {q_N_H, Mass_Balance}; # Parameters # Hill coefficient n = 12; # **** problem, change to nHill to fix... **** q_N_H_max; x; t_on = 2; q_N_N; q_N_E; q_N_C; q_E_C; q_E_E; q_E_EH; q_H_EH; q_H_H; q_EH_H; q_EH_EH; N_N_t0; N_H_t0; N_E_t0; N_C_t0; N_EH_t0; N_transitions; Dynamics { # ***** this is where n is used q_N_H = (q_E_H_on < 0.5 ? 0 : q_N_H_max * pow(t,n) / (pow(x,n) + pow(t,n))); dt(P_N_N) = - P_N_N * (q_N_H + q_N_E + q_N_C); dt(P_N_H) = P_N_N * q_N_H - P_N_H * q_H_EH + P_N_EH * q_EH_H; dt(P_N_E) = P_N_N * q_N_E - P_N_E * (q_E_C + q_E_EH); dt(P_N_C) = P_N_N * q_N_C + P_N_E * q_E_C; dt(P_N_EH) = P_N_H * q_H_EH + P_N_E * q_E_EH - P_N_EH * q_EH_H; dt(P_H_H) = - P_H_H * q_H_EH + P_H_EH * q_EH_H; dt(P_H_EH) = P_H_H * q_H_EH - P_H_EH * q_EH_H; dt(P_E_H) = - P_E_H * q_H_EH + P_E_EH * q_EH_H; dt(P_E_E) = - P_E_E * (q_E_C + q_E_EH); dt(P_E_C) = P_E_E * q_E_C; dt(P_E_EH) = P_E_H * q_H_EH + P_E_E * q_E_EH - P_E_EH * q_EH_H; dt(P_EH_H) = - P_EH_H * q_H_EH + P_EH_EH * q_EH_H; dt(P_EH_EH) = P_EH_H * q_H_EH - P_EH_EH * q_EH_H; } CalcOutputs { Mass_Balance = P_N_N + P_N_H + P_N_E + P_N_C + P_N_EH + P_E_C + P_E_EH + P_E_E + P_E_H + P_H_EH + P_H_H + P_EH_H + P_EH_EH; } End.