Function poly_eval 

pub fn poly_eval(a: &[BigInt], x: &BigInt, p: &BigInt) -> BigInt

Evaluate a polynomial at x mod p (Horner).