intrf.spad line 184 [edit on github]
undocumented
split(p, D) returns [n,s] such that p = n s, all the squarefree factors of n are normal w.r.t. D, and s is special w.r.t. D. D is the derivation to use.
split(p, D)
[n,s]
p = n s
n
w
r
t
D
s