library(PAFit)


net <- generate_BA(N = 100, multiple_node = 20, m = 1)

system.time(b <- to_igraph(net))

#plot(network.extract(b,at = 6))