def compute_dK_mean_dh(A, R, dA_dh, dR_dh, n_manning): return 0.5/n_manning*(dA_dh*R**(2/3) + 2/3*A/R**(1/3)*dR_dh) # friction term is high and constant. This makes ...

k.append(f(x[n - 1] + c[i] * h, u[n - 1] + h * numpy.sum([a[i][j] * k[j] for j in range(i)]))) u.append(u[n - 1] + h * numpy.sum([b[j] * k[j] for j in range(s)])) ...