subroutine compute_polval(n, coef, z, polval)
real (4) coef(0:n)
complex polval, z

polval =coef(n)

do  k =n-1,0,-1
polval =coef(k) +z*polval
enddo

end
