# Test file for HypergeomDeg2:
# Solution of L[12..14] involves gauge transformation.

L:= {Dx^2+(1/2)*(8*x^7-39*x^6-20*x^5+843*x^4-2168*x^3-2225*x^2+11284*x-6755)*Dx/((x^2+2*x-11)*(2*x^2-x-12)*(x^2-7)*(x^2-4*x+7))+
     (225/392)*(x^2-4*x+7)^2/((2*x^2-x-12)*(x^2+2*x-11)^2*(x^2-7)), Dx^2+(3*x^6-23*x^4-9*x^2-63)*Dx/(x*(x^2+1)*(x-3)*(x+3)*(x^2-7))+
     (900/49)*x^2/((x-3)*(x+3)*(x^2+1)^2*(x^2-7)), Dx^2+(1/2)*(15*x^6+58*x^5+275*x^4-68*x^3-2951*x^2-2926*x+805)*Dx/((x+1)*(5*x^2-x-36)*(x^2-7)*(x^2+2*x+7))+
     (225/784)*(x^2+2*x+7)^2/((x^2-7)*(5*x^2-x-36)*(x+1)^2), Dx^2+(2/3)*(6*x^5-4*x^4-7*x^3+x^2-13*x-1)*Dx/((x+1)*(2*x^2-x+2)*(x^2+1)*(x-1))+
     (935/882)*(x-1)^2/((2*x^2-x+2)*(x+1)^2*(x^2+1)), Dx^2+(1/6)*(9*x^4+139*x^3+664*x^2+1164*x+576)*Dx/((x+2)*(x+3)*x*(x+4)*(x+6))-
     (1/144)*(x+6)^2/((x+3)*x*(x+2)^2*(x+4)), Dx^2+(1/6)*(2*x^2+8*x+71)*Dx/((x+2)*(x+5)*(x-1))-(169/16)/((x+2)^2*(x+5)*(x-1)),
     Dx^2+(1/2)*(8*x^3-x^2-118*x+21)*Dx/((x+2)*(x-1)*(2*x^2-10*x-1))+(507/16)/((2*x^2-10*x-1)*(x+2)^2), Dx^2+(1/12)*(15*x^2+68*x-36)*Dx/((x+2)*(x-2)*x)-
     (169/72)/(x*(x+2)^2), Dx^2+(1/3)*(30*x^2+104*x-83)*Dx/((2*x+3)*(x-2)*(4*x-1))-(1183/36)/((4*x-1)*(2*x+3)^2), Dx^2+(1/3)*(108*x^5+729*x^4+2016*x^3+2853*x^2+
     2064*x+610)*Dx/((2*x+3)*(3*x^2+12*x+11)*(x^2-2)*(3*x+4))-(55/36)*(3*x+4)^2/((3*x^2+12*x+11)*(2*x+3)^2*(x^2-2)), Dx^2+(1/3)*(32*x^4+1156*x^3+6720*x^2+
     15939*x+29403)*Dx/((2*x+3)*(4*x+11)*x*(x+9)*(2*x+33))-(55/144)*(2*x+33)^2/((4*x+11)*x*(2*x+3)^2*(x+9)), Dx^2+(1/3)*(32*x^4+1156*x^3+6720*x^2+15939*x+
     29403)*Dx/((2*x+3)*(4*x+11)*x*(x+9)*(2*x+33))+(17/144)*(2*x+33)^2/((4*x+11)*x*(2*x+3)^2*(x+9)), Dx^2+(1/3)*(22*x^4+132*x^3+285*x^2+261*x+
     108)*Dx/((x^2+3*x+4)*x*(x+3)*(2*x+3))-(17/144)*(2*x+3)^2/((x^2+3*x+4)*x*(x+3)), Dx^2-(2/7)*(6*x^4+94*x^3+582*x^2-2395*x-
     1912)*Dx/((x-1)*(2*x^2-7*x-10)*(x+4)*(x+9))+(8525/2352)*(x+9)^2/((2*x^2-7*x-10)*(x-1)^2*(x+4)), Dx^2+(1/7)*(14*x^3+423*x^2-1288*x-
     2774)*Dx/((x-1)*(x+4)*(x^2+28*x+46))+(51475/588)/((x^2+28*x+46)*(x-1)^2), Dx^2-(1/7)*(38*x^2-76*x+9)*Dx/((x-1)*(2*x^2-4*x-1))-
     (2059/588)/((x-1)^2*(2*x^2-4*x-1)), Dx^2+(1/21)*(71*x^2-142*x-28)*Dx/((x-2)*x*(x-1))+(2059/1764)/(x*(x-2)), Dx^2+3*(16*x^2-16*x-19)*Dx/((3*x-1)*(x+2)*(6*x+5))+
     (2059/28)/((6*x+5)*(3*x-1)^2)}:

read"C:/Users/jvijay/Desktop/HypergeomDeg2/HypergeomDeg2.txt":

n:= nops(L);
Ans:= [];
for i from 1 to nops(L) do
  lprint(i);
  L1:= L[i];
  Fld:= indets(L1,{RootOf, nonreal, radical});
  SingsL1:= printsing(L1,Fld);
  Soln:= HypergeomDeg2(L[i],x);  
  if Soln = NULL then return "something went wrong!"; fi;
  Check:= simplify(eval(diffop2de(L1, y(x)), y(x) = Soln[1]));
  if Check <> 0 then return "something went wrong"; fi;
  Ans:= [op(Ans), Check]; nops(Ans);
od;