function dy = odefun(t,y)
dy = zeros(2,1);
rho_salt = 2170; %kg/m^3
rho_water = 1000;
V_jt = 8.333*10^(-5); %m^3/s
V_it = V_jt;
Q = 2500;
Co = 122.652; %kg/m^3
Vol = 0.106; % m^3
m_it =…Statistics and Probability