function []=ueb4(rho_l,m_l,rho_r,m_r,c)
// 
// Scilab-File zur Illustration des
// 4. Uebungsblattes der Vorlesung
// "Theorie und Numerik hyperbolischer Erhaltungsgleichungen"
// M.Ehrhardt, WS 2003/04
// 
// Phasenraum der isothermen Eulergleichungen
// Plot der 1-,2- Schocks/Verduennungswellen
// 
// linker Zustand (rho_l,m_l)
// rechter Zustand (rho_r,m_r)
// Schallgeschwindigkeit c
// 
// Beispielaufruf: ueb4(1,1,3,2,1)
// 
 
holdflag = mtlb_ishold();
// 1 falls "hold on", 0 falls "hold off""
// Plotte linken Zustand u_l (gruen)
mtlb_plot(rho_l,m_l,'go');
mtlb_hold('on');
 
// Plotte rechten Zustand u_r (magenta)
mtlb_plot(rho_r,m_r,'mo');
 
// Berechne 1-Welle
[rho_m,rho_p,m_m,m_p] = kurven(rho_l,m_l,c,1);
S_1 = mtlb_plot(rho_p,m_p,'b-');
V_1 = mtlb_plot(rho_m,m_m,'r--');
 
// Berechne 2-Welle
[rho_m,rho_p,m_m,m_p] = kurven(rho_r,m_r,c,2);
S_2 = mtlb_plot(rho_p,m_p,'b:');
V_2 = mtlb_plot(rho_m,m_m,'r-.');
 
cx = (max(min(rho_l,rho_r)-2,0)+max(rho_l,rho_r)+2)/2;
cy = (max(m_l,m_r)+min(m_l,m_r))/2;
w = max(max(rho_l,rho_r)+2-max(min(rho_l,rho_r)-2,0),max(m_l,m_r)-min(m_l,m_r)+4);
cx = max(cx,w/2);
 
square(cx-w/2,cy-w/2,cx+w/2,cy+w/2) 
 
xstring(rho_l+w/50,m_l-w/50,'u_l');
xstring(rho_r+w/50,m_r-w/50,'u_r');
xtitle(['4.Uebung Theorie und Numerik von hyperbol. Erhaltungsgleichungen';'isotherme Eulergleichungen: Schocks/Verduennungswellen im Phasenraum'],' ',' ');
xtitle(' ','Dichte rho',' ');
xtitle(' ',' ','Impuls(-dichte) m');
 
// "hold" wieder in alten Zusand setzen
if holdflag==%f then
  mtlb_hold('off');
end
 
 



