Plot[frac(n! HarmonicNumber(n!, 3/2)),{n,10,20}]