Focus on the wrong equations. They may explain the BK violation.
0 = A- rho_a*A(-1) - shock_A;
implies a steady state of 0 for A. That does not work with
0 = Y - A*H;
Focus on the wrong equations. They may explain the BK violation.
0 = A- rho_a*A(-1) - shock_A;
implies a steady state of 0 for A. That does not work with
0 = Y - A*H;