Well, there's a closed formula right? I wonder if it's faster to carry out the computation in Z[sqrt(5]).