function res = CGC_2TD_Eqns_V(t,SOL,DSOL)

global Vp Vi tpti td E Gin di Vg

I = SOL(1); % Insulin variable from SOL vector

G = SOL(2); % Glucose variable from SOL vector

T1 = DSOL(2,1); %tau1

T0 = DSOL(1,2); %tau0

dIdt = f1(T1)-di*I;

dGdt = Gin-f2(G)-f3(G)*f4(I)+f5(T0);

res = [dIdt; dGdt]; % return solution vector

end