I am implementing the kinetic model proposed by this paper (modelling diffusion of molecules into different compartments of a biological cell) in R. However; I am not sure if the following code is correct (the main struggle is the exp() function)
A2_t_1_calc <- function(t, C0plus, C0minus, Eplus, Eminus, A1_0, D21, D32, G12, G23, V3, P42) {
y = C0plus * (V3 * Eplus + D32 * G23)*exp(Eplus * t) + C0minus * (V3 * Eminus + D32 * G23)*exp(Eminus * t) + ((D21 * A1_0)/(D21 * G12 + P42))
}
For this formula:
As written, I am unsure if I "translated" the exp function correctly. Therefore, I appreciate the correction and help!