getRDD function

Get density dependent reproduction rate