clear all
close all
clc

set(gca,'FontSize',18)

linewidth = 2;

rho = 0:0.01:1;
plot(rho,1-rho,'Linewidth',linewidth)
xlabel('\rho / \rho_{max}')
ylabel('v/v_0')

print -djpeg velocity

figure
set(gca,'FontSize',18)
plot(rho,rho.*(1-rho),'Linewidth',linewidth)
xlabel('\rho / \rho_{max}')
ylabel('q /(v_0/\rho_{max})')
print -djpeg q_flux

figure
set(gca,'FontSize',18)
plot(rho,1-2*rho,'Linewidth',linewidth)
xlabel('\rho / \rho_{max}')
ylabel('c/(v_0/\rho_{max})')
print -djpeg c_speed