7 | | Given the above properties of the Riemann problem, there is an algebraic expression of which the root is the p in the central region or "star-region" (denoted by '*'), between the L- and R- waves. The overall structure of the Riemann problem then is to solve this algebraic equation for p*. Once p* is known, u* follows immediately. The remaining rho*_L and rho*_R follow from expressions valid for the specific L- or R- wave present. |

| 7 | Given the above properties of the Riemann problem, an algebraic expression can be derived which gives p in the central "star-region" (denoted by '*'). The overall structure of the Riemann problem then is to solve this algebraic equation for p*. Once p* is known, u* follows immediately. The remaining rho*_L and rho*_R follow from expressions valid for the specific L- or R- wave present. |

| 8 | |

| 9 | The program thus will consist of 4 parts: 1. Identification of types of nonlinear waves present in the solution states, 2. Code to find p* using an iterative process, 3. U* given by Toro 4.9, 4. Use appropriate expressions to find rho*_L and rho*_R. |