syz(Matrix) returns nonminimal syzygies

Nobody has claimed this yet.

This is related to the example that produced the bug in #3017. In the [documentation](https://macaulay2.com/doc/Macaulay2/share/doc/Macaulay2/Macaulay2Doc/html/_syz_lp__Matrix_rp.html) for `syz(Matrix)`, the output is "the matrix of minimal or trimmed generators for the syzygies among the columns of h". However, this is not true:
```
i1 : R = QQ[x]

o1 = R

o1 : PolynomialRing

i2 : h = transpose matrix{{x^3+1},{x^2+1}}

o2 = {-3} | x3+1 x2+1 |

1 2
o2 : Matrix R <-- R

i3 : syz h

o3 = | x2+1 x3+x2+x+1 |
| -x3-1 -x4-x3-x-1 |

2 2
o3 : Matrix R <-- R

i4 : mingens image oo

o4 = | -x2-1 |
| x3+1 |

2 1
o4 : Matrix R <-- R
```

I am not sure if the solution here is to fix the code, or to fix the documentation -- the [documentation](https://macaulay2.com/doc/Macaulay2/share/doc/Macaulay2/Macaulay2Doc/html/_syz_lp__Groebner__Basis_rp.html) for `syz(GroebnerBasis)` reiterates that the output might include nonminimal syzygies, but the documentation for `syz(Matrix)` is misleading.

Start by reproducing the linked `syz(Matrix)` example, then compare its documentation with the `syz(GroebnerBasis)` documentation. Determine whether the discrepancy requires a code change or corrected wording; done means the behavior and documentation agree about whether nonminimal syzygies can be returned.

Written by the indexing model from the issue text.