Hi Elena,
You need to normalize the distribution so its integral is 1.
dist = ProbabilityDistribution[
0.001257*x^4 Exp[-0.285714*x],
{x, 0, 60},
Method -> "Normalize"]
RandomVariate[dist, 15]
(*
{27.4943,16.8757,5.30406,6.00207,11.8668,28.2957,21.9066,19.6846,12.266,18.8735,19.3128,
15.404,21.1927,13.8153,19.9035}
*)