


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


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


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


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

  LL := []:

   for j from 1 to 5 do

       f := randpoly([x, t1, t2], degree = i, dense)*1/randpoly([x, t1, t2], degree = i,dense):

       fd := Derivative(K, f):

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

    end do:

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

save L, TestData2:
