LargeComponents <- function(n,lambda){

G<-sample_gnp(n,lambda/n)

C<-components(G)

Cluster_Sizes<-sort(C$csize,decreasing=TRUE)[1:5]

if (lambda<1)
{
rescaling<-log(n)
print(1/(lambda-1-log(lambda)))
}

if (lambda>1)
{
rescaling<-n
zeta_lambda<-1-uniroot(function(x) exp(lambda*(x-1)) - x, lower = 0, upper = 0.99, tol = 1e-9)$root
print(zeta_lambda)
}

return(Cluster_Sizes/rescaling)
}