complex function K1(z)
use work

complex z, q, p1, polval


q =polval(n_q, coef_q(0:n_q), z)

if(q /=0) then
p1 = polval(n-1, coef_1, z)
    K1 =-q/p1

else
K1 =0
endif


end