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

real (8) coef(0:n)

zpolval =coef(n)

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

end
