


read "TestData4":

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


read "Liouville.txt";
with(Liouville);
with(DataList);
with(Risch);
NoNewConst(K);
SetTower(K);



T := [];
                            
for i from 1 to nops(L) do

  

   T1:= 0:

   for j from 1 to 5 do
       print([i,j]);

       fd := L[i,j]:

       TT1 := time():  
      (g, r) := ReductionPair(0, fd):  g:=normal(g):
      TT1 := time() - TT1: T1 := T1 + TT1: print(T1);
       if r <> 0 then
          print(incorect); return fd:
       else
          G := subs([t1 = log(x^2+1), t2 = exp(x^2/2), t3 = exp(exp(x^2/2))], g):
          fD := subs([t1 = log(x^2+1), t2 = exp(x^2/2), t3 = exp(exp(x^2/2))], fd):
          s := normal(diff(G,x) - fD):
          if s <> 0 then return fd: end if: 
       end if:


       
    end do:
    T := [op(T), evalf(T1/5)]; print(T);
end do:

save T, Timing4CR:
