function dxSplitdt = ODESplit(x,mbypass,m4,T4,m1)
% Currently not used
Ms = x(1); T0 = x(2); % Dynamic states

dMsdt = m4 - m1 - mbypass; % Dynamic mass balance over the split to compute Ms
dT0dt = (m4/Ms)*(T4-T0); % Dynamic energy balance to compute T0

dxSplitdt = [dMsdt; dT0dt];
end

