Fernando Duarte

I am trying to run the example from the [Readme.md](https://github.com/SciML/Quadrature.jl#readme) file but computing hessians instead of gradients, and getting an error for both Zygote and ForwardDiff. ```Julia using Quadrature, ForwardDiff,...

```julia using ModelingToolkit A = [[1,0] [0,1]] # identity matrix @variables z[1:4] Symbolics.scalarize(A*z[3:4]) # gives [z_1,z_2] but I was expecting [z_3,z_4]? ``` with ```julia Pkg.status() Status `~/SageMaker/networks_hub/SAN/code/network_simulations/mwe/Project.toml` [961ee093] ModelingToolkit v5.21.0...