This function calculates the degree distribution of the network. First it
fits exponential, power law and truncated power law distribution models,
and calculates the AIC values to select the best fit, and finally it plots
the degree distribution in a log log scale showing the three fitted models
mentioned above against the observed distribution.