Yes, the code works with the correct a_p. This will "re-find" the 1476 in a few minutes:
sage: from psage.modform.hilbert.sqrt5.ellcurve import *
sage: E = EllipticCurve([1,-a-1,a+1, 3904*a-6577, 53806*a-88477])
sage: E2 = find_curve_from_Lfunction(E.conductor(), compute_aplist(E, 30000), 30000, 1, limit1=3, limit2=(0,20), limit3=(0,5), ncpus=12, verbose=1)
The other curve we don't know looks pretty good, though. The periods look good, but it just doesn't find the curve. I'm going to look at it for a few minutes right now