magma_factorization_finiteField
R<x> := PolynomialRing(Integers()); b := x^5 - x^3 + 2*x^2 - 2; F<f> := PolynomialRing(GF(5)); s := [ w[1] : w in Factorization( F ! b ) ];