MHD Shock tub problem ( proposed by Brio and Wu in Journal of Computational Physics, vol. 75, 400-422, 1988 )

Computational domain: x [0, 1],
Initial condition:


with Bx=0.75 and r =2.

Mesh: 801 grid points
Computational time: t = 0.1