#INPUT
pc_filename = "carrierDensity.txt";
#END INPUT

closeall;

pc_mxCarrDens = readdata(pc_filename);
pc_sz = size(pc_mxCarrDens);

pc_position=linspace(0,1,pc_sz(2));
pc_time=linspace(0,pc_sz(2),pc_sz(1));

image(pc_time,pc_position,pc_mxCarrDens,"time","position (left port =0)","carrier density [1/m^3]");

plot(pc_position,pc_mxCarrDens(1,1:pc_sz(2)),"position (left port =0)","start carrier density [1/m^3]");

plot(pc_position,pc_mxCarrDens(pc_sz(1),1:pc_sz(2)),"position (left port =0)","end carrier density [1/m^3]"); 

clear(pc_filename,pc_mxCarrDens,pc_sz,pc_position,pc_time);