


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


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




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