I intend to use LinearAlgebra package to do some calculations. I want to compute the basis for large Matrices. My discovery is that the linalg[kernel] command, which the document claims is deprecated, could do such computation significantly faster than the LinearAlgebra[NullSpace] command. For a 200 x 500 large random matrix, linalg[kernel] clocked 33 secs, while the LinearAlgebra[NullSpace] takes 200 secs, as shown in the worksheet NullSpace_vs_kernel.mw.

I wanna know what makes the difference, or is there a misuse for LinearAlgebra[NullSpace].

 Based on the type of coefficients in the linear equations, SolveTools[Linear] provide several method including method = Rational, Polynomial etc. 
The Polynomial method of SolveTools[Linear] however, cannot be directly called by SolveTools[Linear]. The interpreter complains that no such methd called "Polynomial".

This could only be solved if I import the SolveTools package in advance and call Linear.

I hope Maple could solve this issue. 

So just like the title illustrates, I found a paper authored by Gary Nicklason in 2022: Autonomous Planar Systems of Riccati Type and in the last section it mentioned about a class of Abel ODE, which belongs to AIA(Abel Inverse Abel) class. It is of First kind and the inverse of it(by swapping variables) is of second kind.

While the first kind is solvable in terms of Airy function, the inverse of it along with its equivalence class is not solvable by the existing dsolve.

I have tested it in my worksheet Nicklason_equation.mw. So is it possible to add this class into the dictionary for solvable Abel ODE, or, maybe there are some bugs within the internal procedure of dsolve, which results in failure for catching the solvable candidates?

I was trying to solve a system of eight cubic equations, with eight variables. Note that the particular solutions should exist and my goal is to find all  possible solutions using solve. However, when executing, the solve keeps running for a whole day and did not throw any results. I also set the parameter infolevel[solve] to be 3 and find out that it is stuck in the step "GroebnerBasis: computing a factored plex basis using Groebner[Solve]". Can anyone tell me how to deal with that? Here's the Maple file.

solve_test.mw

I was trying to solve a system of polynomial equations, which contains three equations and six variables $a_0,a_1,a_2,b_0,b_1,b_2$. However, as I swap the variable name, Maple solve function gives me a totally different solutions. Only the solutions before swapping the variables are useful for the problem I study. I have already attached the file. Could anyone tell me if the choice of variable name really matters, or if i just misuse this function?

Choice_of_name_infolevel.mw