I want to change the F value at those timepoints how do i integrate it here is my code thank u:)
When i simulate the F values are 0
`def bior(t, y):
V,F,u,ud,Qs,Qmab,Qlac, xv, xd, s, mab, lac = y
T_values = [2.03, 5.02,6.03,7.03,7.88,9.03, 10.02,12.02]
F_values = [0.2,0.194,0.148,0.156,0.128,0.119,0.232,0.235]
# Interpolate the flow value at the current time
Fin = Fout = np.interp(t, T_values, F_values)
dV = Fin-Fout
du = umaxap*(s/(kglc+s))*(Kilac/(Kilac+lac))
dud = udmax*(lac/(Kilac+lac)) #udmax / (1 +((Kilac/lac)**1))
dQs= (u-ud)/Yxglc
dQmab=(u-ud)/Yxmab
dQlac= ylp*Qs
dxv =(umaxap*(s/(kglc+s))*(Kilac/(Kilac+lac)))*xv-(udmax / (1 +((Kilac/lac)**1)))*xv
dxd = (udmax / (1 +((Kilac/lac)**1))) * xv - Klys * xd
ds = -(Qs)*xv+((Fin*Cs)/V)
dmab = (umaxap*(s/(kglc+s))*(Kilac/(Kilac+lac)))*xv-(udmax / (1 +((Kilac/lac)**1)))/Yxmab* xv
dlac = ylp*Qs* xv
is
return [dV,F,du,dud,dQs,dQmab,dQlac,dxv,dxd,ds,dmab,dlac]
y0 = [15, 0,0,0,0,0,0, 1502000000.00, 19000000.00, 5, 0, 0.000001]
t= (0, 325)`
New contributor
user23827055 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.