Durand-Kerner method