function Q = water_sediment_flux(Q0, S0, t0, T1, dt)
% 计算水沙通量
% 输入参数：
% Q0 - 初始流量 (m^3/s)
% S0 - 初始沙子浓度 (kg/m^3)
% t0 - 初始时间 (s)
% T1 - 最终时间 (s)
% dt - 时间步长 (s)
%
% 输出参数：
% Q - 最终流量 (m^3/s)

% 计算时间步数
N = int((T1 - t0) / dt) + 1;

% 初始化流量和沙子浓度数组
Q = zeros(1, N);
S = zeros(1, N);
S(1) = S0;
Q(1) = Q0;

% 迭代计算水沙通量
for i = 2:N
    dt = (i-1)*dt; % 更新时间步长
    dS = S(i-1) * dt; % 沙子浓度变化率
    dQ = Q(i-1) * dt; % 流量变化率
    Q(i) = Q(i-1) + dQ; % 更新流量
    S(i) = S(i-1) + dS; % 更新沙子浓度
end
end