% 
% Written by:
% -- 
% John L. Weatherwax, Ph.D.         2003-02-11
% MIT Lincoln Laboratory
% 244 Wood Street
% Lexington, MA 02420-9176 USA
% 
% email: weatherwax@ll.mit.edu
% Voice: (781) 981-5370       Fax: (781) 981-0721
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

% (a):
q_l = [ 0 1 ];
q_r = [ 1 1 ];
A = [ 0 4; 1 0 ];

[ states, waves ] = riemann2x2( A, q_l, q_r )

disp( 'Press return for the next experiment' );
pause

% (b):
q_l = [ 1 1 ];
q_r = [ 0 1 ];
A = [ 0 4; 1 0 ];

[ states, waves ] = riemann2x2( A, q_l, q_r )

disp( 'Press return for the next experiment' );
pause

% (c):
q_l = [ 1 0 ];
q_r = [ 4 0 ];
A = [ 0 9; 1 0 ];

[ states, waves ] = riemann2x2( A, q_l, q_r )

disp( 'Press return for the next experiment' );
pause

% (d):
q_l = [ 1 0 ];
q_r = [ 2 0 ];
A = [ 1 1; 1 1 ];

[ states, waves ] = riemann2x2( A, q_l, q_r )

disp( 'Press return for the next experiment' );
pause

% (e):
q_l = [ 0 1 ];
q_r = [ 1 0 ];
A = [ 2 0; 0 2 ];

[ states, waves ] = riemann2x2( A, q_l, q_r )

disp( 'Press return for the next experiment' );
pause

% (f):
q_l = [ 0 1 ];
q_r = [ 1 0 ];
A = [ 2 1; 10^(-4) 2 ];

[ states, waves ] = riemann2x2( A, q_l, q_r )