degree of p: 4
positive roots in (0, L_tot): 3 ; roots in (0,inf): 3
squarefree (gcd(p,p')==1): True
  root in [0.0096, 0.0096] mult 1
  root in [0.08672435511, 0.08672436015] mult 1
  root in [410, 410] mult 1
