modified to first plot exact solution alone and then approximate solutions.  Fiddled with q_l and q_r