Problem 2.21

Estimate the multiplicity function to the two-Einstein-solid function

Clear[f, n, z]

f[z_, n_] := (((4 z (1 - z)))^n)

Plot[Evaluate[Table[f[z, n], {n, {10, 100, 1000}}]], {z, 0, 1}, PlotRange -> {0, 1}, PlotLegends -> {"10", "100", "1000", "10000"}]