pell.eq_pell source