Is it possible to modify lagrange_polynomial so that it accepts points with variable entries? I was able to make it work on a small example, by disabling the coercion functions in the source file "polynomial_ring.py", but then the program failed on a finite field.