function [] = prob_3_3_b() % % Written by: % -- % John L. Weatherwax 2006-10-16 % % email: wax@alum.mit.edu % % Please send comments and especially bug reports to the % above email address. % %----- close all; % Compute the shock/CD speeds: % lambda1 = 1; lambda2 = 2; lambda3 = 3; % Draw the x-t diagram solution to this Riemann problem (at time="time"): % time=1.0; figure; hold on; line( [ 0; lambda1*time ], [ 0 ; time ] ); line( [ 0; lambda2*time ], [ 0 ; time ] ); line( [ 0; lambda3*time ], [ 0 ; time ] ); text( .25, .8, 'q_l' ); text( 1.0, 0.8, 'q_{ml}' ); text( 1.5, 0.6, 'q_{mr}' ); text( 2.0, 0.3, 'q_r' ); xlabel( 'x' ); ylabel( 't' ); %grid on; saveas( gcf, '../LeVeque_WU_FVM/Graphics/prob_3_3_b_xt.eps', 'psc2' );