Estou trabalhando com pyomo e já tenho um modelo definido, com uma função objetiva para acompanhá-lo. Depois que o modelo é resolvido, a função objetivo possui determinados parâmetros anexados. Então, se eu tivesse uma variável com vários índices[x1, x2, x3], minha função objetiva quadrática suporia algo como isto:(x1^2 + 13*x2^2 + 10*x3^2) + (2*x1 +......) .

Minha pergunta é: dado que eu posso realmente acessar essa expressão no formato de string a partir do objetivo, existe alguma maneira de obter a segunda derivada dessa função em relação a todas as variáveis?