I am trying to solve using the Fenics solver the heat transfer equation considering temperature-dependent parameters. However, I am getting the following error:

It seems that the user-defined function
density needs to be included somehow in the solver. Any thoughts on how this can be fixed? I am using FEATool Multiphysics as a Matlab toolbox on Linux.