! Cu "polval" pt. polinoame
! met. Halley
subroutine Deriv_Pol()
use work


do i = n-1, 0, -1
coef_1(i) =coef(i+1)*(i+1)
enddo

!f  =polval(n, coef, x)
!f1 =polval(n-1, coef_1, x)

END
