  This article presents an original code for Big Bang Nucleosynthesis in a
baryon inhomogeneous model of the universe. In this code neutron diffusion
between high and low baryon density regions is calculated simultaneously with
the nuclear reactions and weak decays that compose the nucleosynthesis process.
The size of the model determines the time when neutron diffusion becomes
significant. This article describes in detail how the time of neutron diffusion
relative to the time of nucleosynthesis affects the final abundances of He4,
deuterium and Li7. These results will be compared with the most recent
observational constraints of He4, deuterium and Li7. This inhomogeneous model
has He4 and deuterium constraints in concordance for baryon to photon ratio eta
= (4.3 - 12.3) X 10^{-10} Li7 constraints are brought into concordance with the
other isotope constraints by including a depletion factor as high as 5.9. These
ranges for the baryon to photon ratio and for the depletion factor are larger
than the ranges from a Standard Big Bang Nucleosynthesis model.
