Benchmark solution ===================== The reference solution corresponds to model B.