% TransportPlot.m 
function TransportPlot(xmin, dx, xmax, dt, tmax, c)
sol = transport(xmin, dx, xmax, dt, tmax, c, 'f0transp');
subplot(2,2,1)
plot(xmin:dx:xmax, sol(:,1))
axis([xmin xmax -0.1 1.1])
subplot(2,2,2)
plot(xmin:dx:xmax, sol(:,50))
axis([xmin xmax -0.1 1.1])
subplot(2,2,3)
plot(xmin:dx:xmax, sol(:,100))
axis([xmin xmax -0.1 1.1])
subplot(2,2,4)
plot(xmin:dx:xmax, sol(:,200))
axis([xmin xmax -0.1 1.1])