clc T=300 //K k=8.617*10^-5 //eV/K q=1.6*10**-19 //C epsilonx=8.854*10^-14 //F/cm L=10*10^-4//cm^-2 v=10^7//sec n0=10^12/L disp(n0,"n0 in cm^-3 is=") t=L/v disp(t,"t in sec is=")