1, J.Chem.Phys. paper, 2007: Electrodiffusion: A continuum modeling framework for biomolecular systems with realistic spatiotemporal resolution,
BZ Lu, et. al., J. Chem. Phys., Vol. 127:135102, (2007)

	1), in Eq. (5), it is better to add \lambda(r) for consistence.
	2), Eq. 6, \Omega_m -> \Omega.
	3), The intrinsic reaction rate k(r) in the boundary conditions
	    has not been noted, and this symbol is the same
	    as the reaction rate coefficient k used later.
	    It is better to use \alpha.

	4) P.7, 
	"The condition jn=0 is actually a Neumann condition for the
	pure diffusion equation or, more generally, a type of Robin
	condition", this sentence not proper. Robin condition
	is not a more general BD than Neumann BD...
	
2, CiCP paper: Recent Progress inNumericalMethods for the Poisson-
  	Boltzmann Equation in Biophysical Applications,
	BZ. Lu, et al.,, 3:973-1009, (2008).

	(1), in Eq. (1.1), it should be minus, instead of plus, in fornt of
		\lambda.