function xk1 = Holos_microreactorStateFcn(xk, u)
uk = u(1);
Ts = u(2);
xk1 = Holos_microreactorDT0(xk, uk, Ts);
