A (u 1 , u 2 ,...)-parking function of length n is a sequence (x 1 , x 2 ,..., x n ) whose order statistics (the sequence (x (1) , x (2) ,..., x (n) ) obtained by rearranging the original sequence in non-decreasing order) satisfy x (i) ≤ u i . The Goncarov polynomials g n (x; a 0 , a 1 ,..., a n−1 ) are polynomials biorthogonal to the linear functionals ε(a i )D i , where ε(a) is evaluation at a and D is differentiation. In this paper, we give explicit formulas for the first and second moments of sums of u-parking functions using Goncarov polynomials by solving a linear recursion based on a decomposition of the set of sequences of positive integers. We also give a combinatorial proof of one of the formulas for the expected sum. We specialize these formulas to the classical case when u i = a + (i − 1)b and obtain, by transformations with Abel identities, different but equivalent formulas for expected sums. These formulas are used to verify the classical case of the conjecture that the expected sums are increasing functions of the gaps u i+1 − u i . Finally, we give analogues of our results for real-valued parking functions.