! Rad tripla:
! f =(x-1)**3*(x-2) =x**4 -5*x**3 +9*x**2 -7*x +2
! n =4; coef(n:0:-1) =(/1, -5, 9, -7, 2/);

! f =(x-1.1)**2*(x-2.1)
! n=3; coef(n:0:-1) =(/1,	-4.3,	5.83,	-2.541/)

! Zeng 2003: n = 10
!Zeng 2003: (x-3)^2*(x-2)^3*(x-1)^5
!n = 10; 1 -17 127 -549 1521 -2823 3557 -3007 1634 -516 72