# Reading the file hypergeometricsols_source in Maple produces the file "hypergeometricsols.m"

# Reading that latter file loads the implementation:
read "hypergeometricsols.m";
read Examples:
Tests := [ Linp_walks, Linp_alg, Linp_nonlog, Linp_log ]:
_Envdiffopdomain := [Dx,x]:

Test3 := false; # Set to false if you want to run all tests (about 10 minutes)

if Test3 = true then
	# test only 3 examples of each type
	Tests := map(i -> i[1..3], Tests)
fi:
Tests := map(op, Tests):

for L in Tests do
	H := `DEtools/hypergeometricsols`(L);
	if H = [] then
		error "not solved"
	fi
od;

# We should also check for some examples if the output is really a solution of the input:

# Modify one of examples from Linp_walks:
L := DEtools[symmetric_product]( Linp_walks[-1], Dx+x+3-3/x^2);
L := DEtools[diffop2de](L, y(x));
S := `DEtools/hypergeometricsols`(L);
simplify(eval(L, y(x) = S[1]));
simplify(eval(L, y(x) = S[2]));