function u = SRC_tank1(v)
alpha_1 = -0.017990;
alpha_2 = 0.022516;
v_split = 44.41; % Split value of v
u=[1 + alpha_2 * (v - v_split); 1 + alpha_1 * (v - v_split)];
% Nominal point: v = v_1 = (0.5-1)/alpha_2 + v_split = 22.2036
% u0 = [;]
% v
end