In [ ]:
library(igraph)
ba.net <- barabasi.game(100000)
ba.net.deg.dist <- degree.distribution(ba.net)
ba.net.deg <- degree(ba.net)
plot(0:max(ba.net.deg),
      ba.net.deg.dist,
    type="l",
    xlab=expression(k),
     ylab=expression(p[k]),
     yaxt="n",
    log="xy")

In [ ]:
ba.net.small <- barabasi.game(100)
plot.igraph(ba.net.small, vertex.size=2)