function [] = prob_13_11_b()
% 
% Written by:
% -- 
% John L. Weatherwax                2006-08-28
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

% Holds the states (v_star, u_star): 
q_star = [ .1 .25 .5 .75 1 ];
u_star = q_star; 

q = linspace( 0.01, 1, 100 ); 

figure; hold on; 

% plot one-integral/Hugoniot curves:
%
for ii=1:length(q_star),
  qu_prod = q_star(ii)*u_star(ii); 
  u1 = qu_prod./q; 
  plot( q, u1, '-b' ); 
end

% plot two-integal/Hugoniot curves:
%
for ii=1:length(q_star),
  u2 = repmat(u_star(ii),1,length(q)); 
  plot( q, u2, '-r' ); 
end

axis( [ 0 1 0 1.25 ] ); 
xlabel( 'q' ); ylabel( 'u' ); grid on; 
saveas( gcf, 'prob_13_11_pt_b.eps', 'psc2' );