diff=[c for c in roots if mask(c,a,b)!=3]
dc=base[:]+diff
db=base[:]
u,dc,db,dx,dpi=runnode(dc,db,a,b,2,u)
print('DIFF',obj(db,dx),integ(db,dx),[(c['cid'],c['cov'],dx[k]) for k,c in enumerate(db)])
print('DD',DD[0],'HS',HS)