

K := x, [t1, t2, t3], [2*x/(x^2+1), x*t2, x*t2*t3];



read "Liouville.txt";
with(Liouville);
with(DField);




L := [];
                            
for i from 10 to 20 do

  LL := []:

   for j from 1 to 5 do

       f := randpoly([x, t1, t2, t3], degree = i, dense):

       fd := Derivative(K, f):

       LL := [op(LL), fd]:

    end do:

   print(i);
    
   L := [op(L), LL]:
end do:

save L, TestData4:
