! Polinomul dat (p0), si imbricat (p1) - definit ca functie
subroutine sub_pol(x, p0, p1)

! cp0: (x-1)^3
!p0 =x**3 -3*x**2 +3*x -1
!p1 =1 +x*(-3 +x*(3 -x))

! cp1: (x-1) ... (x-7)
p0 =x**7 -28.*x**6 +322.*x**5 -1960.*x**4 +6769.*x**3 -13132.*x**2 &
+13068.*x -5040.
p1 =-5040+x*(13068+x*(-13132+x*(6769+x*(-1960+x*(322+x*(-28+x*1))))))

!cp2: (x-1.1)^3*(x-2.1)
!p0 =x**4 -5.4*x**3 +10.56*x**2 -8.954*x +2.7951
!p1 =2.7951 +x*( -8.954 +x*(10.56 +x*(-5.4 +x)))


END

