

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



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, TestData6:
