better fix

Author: ChrisRackauckas

test/reduction.jl

@@ -175,8 +175,8 @@ nlprob.f(residual, reducedsol.u, pp)
 N = 5
 @variables xs[1:N]
 A = reshape(1:(N^2), N, N)
-eqs = xs .~ A * xs
-@named sys′ = NonlinearSystem(collect(eqs), [xs], [])
+eqs = xs ~ A * xs
+@named sys′ = NonlinearSystem(eqs, [xs], [])
 sys = structural_simplify(sys′)
 
 # issue 958