With[{data = {4515, 5974, 7717, 9692, 11948, 14380, 17205, 20440,
24324, 28018, 31161, 34546, 37198, 40171, 42638, 44653, 46472,
48467}},
NonlinearModelFit[data,
Exp[a x^2 + b x + c], {{a, 0}, {b, 0.3}, {c, 8}}, x,
Weights -> 1/data]]["RSquared"]
(* 0.99982 *)