-adjusted fit_linear_diffusion_data() to return data arrays for plotting the linear fit to the data (to match similar functionality in the non-linear function)
-updated the docs to reflect above behaviour