Computing Galois Groups of Completely Reducible Differential Equations

We give an algorithm to calculate a presentation of the Picard–Vessiot extension associated to a completely reducible linear differential equation (i.e. an equation whose Galois group is reductive). Using this, we show how to compute the Galois group of such an equation as well as properties of the Galois groups of general equations.