function [dxdt]=ODEs_CL_GenI(x,d,p,u)

h_2 = x(1);
W_a4 = x(2);
W_b4 = x(3);

q_1 = d(1);
q_2 = d(2);

q_1e_nom = p(21);
h_2_nom = p(30);

Cv_1 = q_1e_nom/h_2_nom^(0.5);
W_a1 = p(23);
W_b1 = p(24);
W_a2 = p(25);
W_b2 = p(26);
W_a3 = p(27);
W_b3 = p(28);
%q_3 =  p(20);
A_1 = p(5);
A_2 = p(6);
h_1 = p(29);

q_3 = u;

q_1e = Cv_1 * h_2^0.5;

dxdt(1) = (q_1 - q_1e)/A_2;
dxdt(2) = (q_1e*(W_a1-W_a4)+q_2*(W_a2-W_a4)+q_3*(W_a3-W_a4))/(A_1*h_1);
dxdt(3) = (q_1e*(W_b1-W_b4)+q_2*(W_b2-W_b4)+q_3*(W_b3-W_b4))/(A_1*h_1);

end