text: Find the intervals with only one root using this function: def numberOfRoots(P, range_a, range_b):
    seq = SturmSequence(P)
  
    pos = 0
    neg = 0
    print(seq)
    for i in range(range_a, range_b):
        eval_poly = [np.polyval(poly, i) for poly in seq]  
        n = np.size(eval_poly) - 1 
        for j in range(0, n):
          if eval_poly[j] > 0 and eval_poly[j+1] < 0:
              pos = pos + 1
          if eval_poly[j] < 0 and eval_poly[j+1] > 0:
              neg = neg + 1       
                     
    return pos, neg