Removed printing to screen
Working model, implements SAHA DIST. FUNCTION: Gives population density at a given r