2       ic
 5.e0    beta for initial pulse

 1.e0      rhol
 1.e0      cl
 4.e0      rhor
 0.5d0     cr