(PARI) {a(n) = if(n==0, 1, 5^(n*(n-1)/2) + 5^(n*(n+1)/2))}
for(n=0, 15, print1(a(n), ", "))
(PARI) /* As Coefficients in a Jacobi Theta Function: */
{a(n) = polcoeff( sum(m=-n, n, 5^(m*(m+1)/2)*x^(m^2) +x*O(x^(n^2))), n^2)}
for(n=0, 15, print1(a(n), ", "))
(PARI) /* By the Jacobi Triple Product identity: */
{a(n) = polcoeff( prod(m=1, n^2\2+1, (1 - 5^m*x^(2*m)) * (1 + 5^m*x^(2*m-1)) * (1 + 5^(m-1)*x^(2*m-1)) +x*O(x^(n^2))), n^2)}
for(n=0, 15, print1(a(n), ", "))