function zpolval(n, coef, z)
complex(16) zpolval, z

real(16) coef(0:n)

zpolval =coef(n)

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

end
