function [rho_m,rho_p,m_m,m_p]=kurven(rho_l,m_l,c,p)
rho_m=[];rho_p=[];m_m=[];m_p=[];
// berechnet 1-,2- Schocks/Verduennungswellen
// der isothermen Eulergleichungen
// 
// Numerierung bezieht sich auf Buch von Leveque
// 
if p==1 then
  cc = c;
elseif p==2 then
  cc = -c;
else
  error('p: invalid value!');
end
 
xi = 0:0.1:10;
%v=size(xi)
rho_p = rho_l*(ones(%v(1),%v(2))+xi);
// (7.11)
// (7.12) bzw. (7.13)
%v=size(xi)
%v=size(xi)
m_p = m_l*(ones(%v(1),%v(2))+xi)-cc*rho_l*xi .* sqrt(ones(%v(1),%v(2))+xi);
 
if p==1 then
  rho_m = rho_l*exp(-xi/cc);
  // (8.13)
  %v1=size(xi)
  m_m = rho_l*(xi+m_l/rho_l*ones(%v1(1),%v1(2))) .* exp(-xi/cc);
  // (8.15)
elseif p==2 then
  // Analoge Formeln (1),(2) aus der 3.Aufgabe der 2.Uebung:
  rho_m = rho_l*exp(xi/cc);
  %v1=size(xi)
  m_m = rho_l*(-xi+m_l/rho_l*ones(%v1(1),%v1(2))) .* exp(xi/cc);
else
  error('p: invalid value!');
end






