transportplotok.m

来自「Matlab codes for financial models」· M 代码 · 共 19 行

% 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])