# ReduceLODE(F, L, UP, A, LO)¶

reduceLODE(op, g) returns [m, v] such that any solution in A of op z = g is of the form z = (z_1, ..., z_m) . (b_1, ..., b_m) where the b_i's are the basis of A over F returned by basis() from A, and the z_i's satisfy the differential system M.z = v.
reduceLODE(op, [g1, ..., gn]) returns [m, [v1, ..., vn]] such that any solution in A of op z = c1*g1 + ... + cn*gn where ci are constants satisfy the differential system M.z = c1*v1 + ... + cn*vn