  In the present work we are going to add a correction to the BHF potential
calculation by introducing a two- body density dependent potential that acts as
a three body interaction. Adding the result of this potential to the BHF
potential would give us new saturation curves which coincide with the results
obtained from the experiment. We compute the energy per particle versus density
curves for symmetric nuclear matter, pure neutron matter and the symmetry
energy using modern nucleon nucleon potentials like CD Bonn.
