We can do this in Julia quite well. ModelingToolkit.jl for example takes in numerical Julia functions and symbolically calculates and builds the function for co

We can do this in Julia quite well. ModelingToolkit.jl for example takes in numerical Julia functions and symbolically calculates and builds the function for computing the Jacobian. Example: https://github.com/JuliaDiffEq/ModelingToolkit.jl/pull/155 https://github.com/JuliaDiffEq/ModelingToolkit.jl/pull/155 The JuliaDiffEq team will be pushing symbolics for automatic transformations and automated model order reductions over the next few years. It's already working really well.