function dxdt=vteamRHSRegions(t,x,p)
Wx4=p.Gon*p.Goff/(p.Gon-(p.Gon-p.Goff)*x(4));
Wx5=p.Gon*p.Goff/(p.Gon-(p.Gon-p.Goff)*x(5));
dxdt1=-(p.alpha/p.G)*(Wx4+Wx5+p.GN2)*x(1)-p.alpha*x(3);
dxdt2=p.gamma*x(2)+x(3);
dxdt3=p.beta*(x(1)-x(2)-x(3));
if p.Region==1 || p.Region==4 || p.Region==5
  dxdt4=(p.eta1/p.i0)*Wx4*x(1);
else
  dxdt4=0;
end
if p.Region==1 || p.Region==2 || p.Region==3
  dxdt5=(p.eta2/p.i0)*Wx5*x(1);
else
  dxdt5=0;
end
dxdt=[dxdt1 dxdt2 dxdt3 dxdt4 dxdt5]';
