function fx = fx_reactor(t,x)
% fx = fx_reactor(t,x)

k1=50; k2=100; k3=10;


%u=25;
if t < 0.15
    u=25;
elseif t > 0.15
    u=30;
end
v=10;

x1=x(1);
x2=x(2);

f1=-k1*x1-k3*x1^2+(v-x1)*u;
f2=k1*x1-k2*x(2)-x2*u;

fx=[f1;f2];