Estoy trabajando con pyomo y ya tengo un modelo definido, con una función objetivo que lo acompaña. Una vez resuelto el modelo, la función objetivo tiene ciertos parámetros adjuntos. Entonces, si tuviera una variable de índice múltiple[x1, x2, x3], mi función objetivo cuadrática supondría un aspecto similar a este:(x1^2 + 13*x2^2 + 10*x3^2) + (2*x1 +......) .

Mi pregunta es: dado que realmente puedo acceder a esta expresión en formato de cadena desde el objetivo, ¿hay alguna forma de obtener la segunda derivada de esta función con respecto a todas las variables?