Dear Gap forum, Does anyone have an implementation in GAP of Lagrange’s method <https://encyclopediaofmath.org/wiki/Lagrange_method> for reducing a quadratic form with coefficients in a field (of characteristic other than 2) to a sum of squares? Thanks, Daniel Ruberman