Only order 2 is implemented so far.  However, the most interesting
case is order 3, because for order 2 one can do something similar
to the differential case, namely: for L = a2 tau^2 + a1 tau + a0
compute the rational solutions of
 symmetric_product(  symmetric_power(L,2),  tau - a2/a0  )