answered 13 years ago
The nice solution:
sage: model.substitute_expression(*sol) x |--> -2.0000000000000031*x^3 + 14.500000000000021*x^2 - 31.500000000000046*x + 23.000000000000028