src/algebra/rdeefx.spad

561         ansp := eval(res1.primpart, nxk, kf)
562         [a1, r1, ansp, res1.sol?]
1399         ansp : F := 0
1405             ansp := ansp + csol(i)*first(cand_prims)
1407         [a1p, rr1f, ansp, true]
1473                        r1 : F, g : F, ansp : F, scoeffs : List(RF),
1493         res2.sol? => [a1 + a1p, g, ansp, true]
1504             neg_deg > 0 => [a1, r1, ansp, false]
1505             [a1 + a1p, g - rr1f, ansp, false]
1506         [a1 + res3.ans, g, ansp + res3.primpart, true]
1524           a1 : F, g : F, ansp : F, derivation : UP -> UP,
1531         (spec1p := res2.remainder) = 0 => return [a1, g, ansp, true]
1547         uvu0 case "failed" => [a1, r1, ansp, false]
1583         sol_rec.particular case "failed" => [a1, r1, ansp, false]
1587                 if sol0(1)^2 ~= sol0(2) then return [a1, r1, ansp, false]
1589             rcu case "failed" => [a1, r1, ansp, false]
1601                        ansp, scoeffs, cands, 0, derivation)
1610             cis case "failed" => [a1, r1, ansp, false]
1624             [a1, r1, ansp, false]
1627         [a1, r1, ansp, false]
1701         ansp : F := 0
1703             ansp := ansp1
1706                 ansp := ansp + eval(ansp1, rf1k, zz)
1715         [ansp, lp_cor, poly_cor, sp_cor]
1801         ansp : F := 0
1803             ansp := ansp1
1807                 ansp := ansp + eval(ansp1, rf1k, zz)
1809         [ansp, dvr2.remainder::RF/denom(lp_cor)::RF, dvr2.quotient, 0]
1817         ansp : F := 0
1828             ansp := ansp + ar1.primpart
1832         [ansp, lp_cor, poly_cor, sp_cor]
1853         ansp : F := 0
1862                 ansp := ansp + ar1.primpart
1866             [ansp, lp_cor, poly_cor, sp_cor]
1900         ansp : F := 0
1913                 ansp := ansp + ar1.primpart
1948             [ansp, lp_cor, poly_cor, sp_cor]
1953             ansp := ansp + ar1.primpart
1957         [ansp, lp_cor, poly_cor, sp_cor]
2010         ansp : F := 0
2039         ansp : F := 0
2050             ansp := ansp + res1.primpart
2053         [r1, rh, ansp, is_sol]
2258         ansp := ar1.primpart
2265         her.polypart = 0 and her.specpart = 0 => [a1, g, ansp, true]
2273             [a1, r1, ansp + res5.primpart, res5.sol?]
2284         rr1 = 0 => [a1, g, ansp, true]
2286             quadratic_log_erf_integrate(nf, nfp, rr1, a1, g, ansp,
2303             [a1 + res3.ans, r1 + res3.right, ansp + res3.primpart, res3.sol?]
2326                     return [a1, r1, ansp, res7.sol?]
2327             ff = 0 => return [a1, r1, ansp, true]
2340             ansp := ansp + res4.primpart
2342             return [a1, r1, ansp, res4.sol?]
2343         r1 = g => [a1, g, ansp, true]
2345         [a1 + res6.ans, r1 + res6.right, ansp, res6.sol?]