Code: Select all
DERIVATIVE blrrate {
calcParameters(blr, aG, cAMP, aCaMK, Vcilia, CaCaM, Ca, IX)
blr' = k1 * od * (Rtot - blr) - r1 * blr
aG' = kG * (Gtot - aG) - rG
cAMP' = syncAMP - pd * cAMP
Ca' = inf * Icng - JNCX - (cc1 - cc2 * CaCaM)
CaCaM' = cc1lin * Ca - cc2 * CaCaM
aCaMK' = ck1lin * CaCaM - ck2 * aCaMK
IX' = cx1lin * Ca - cx2 * IX
Vcilia' = (Icng+IClCa+IL)/ CapCilia
}
When I debugged with printf statement, only the first 2 equations are getting solved. The rest of the variables do not have values.
t: 4.8, blr: 7.25551e-011, aG: 7.26423e-007, cAMP: -1.#IND, Ca: -1.#IND, CaCaM: -1.#IND, aCaMK: -1.#IND, IX: -1.#IND, Vcilia: -1.#IND
t: 4.825, blr: 6.10429e-011, aG: 6.52792e-007, cAMP: -1.#IND, Ca: -1.#IND, CaCaM: -1.#IND, aCaMK: -1.#IND, IX: -1.#IND, Vcilia: -1.#IND
t: 4.85, blr: 5.13573e-011, aG: 5.86624e-007, cAMP: -1.#IND, Ca: -1.#IND, CaCaM: -1.#IND, aCaMK: -1.#IND, IX: -1.#IND, Vcilia: -1.#IND
t: 4.875, blr: 4.32085e-011, aG: 5.27162e-007, cAMP: -1.#IND, Ca: -1.#IND, CaCaM: -1.#IND, aCaMK: -1.#IND, IX: -1.#IND, Vcilia: -1.#IND
t: 4.9, blr: 3.63526e-011, aG: 4.73726e-007, cAMP: -1.#IND, Ca: -1.#IND, CaCaM: -1.#IND, aCaMK: -1.#IND, IX: -1.#IND, Vcilia: -1.#IND
Could you please provide some pointers as to why this is happening? The same equations work fine in Matlab.