


%
p=DDGhole_driftdiffusion(mesh,Dsides,nin,pin,V,un,tn,tp,n0,p0,weight)
IN:
V = electric potential
mesh = integration domain
nin = electron density in the past
pin = hole density in the past + initial guess
n0,p0 = equilibrium densities
tn,tp = carrier lifetimes
weight = BDF weights
up = mobility
OUT:
p = updated hole density
%