Let's consider a modified example from Sage reference manual
f(x,y)=x^2*y+y^2+y
solutions=solve(list(f.diff()),[x,y])
the_solution=solutions[2]
H=f.diff(2); # Hessian matrix
How can i properly substitute the_solution into H?
I have tried:
H(x,y).subs(the_solution) - does not work.
This will work for H(x,y).subs(x==0)
H(x,y).subs_expr(*the_solution) - does not work.
This will work for f(x,y).subs_expr(*the_solution)