# Problem 14:
#
print( polyroot( c(-16, -8, 1, 24) ) )

x = seq( -1, +2, length.out=100 )
y = 24 * x^3 + x^2 - 8 * x - 16
plot( x, y, type='l' )
grid()