[GreyMat]

-transient state condition
- density is function of liquid fraction
- liquid fraction is function of temperature, hence density is the function of temperature
- coded a user define function for density
- it worked fine but soon temperature limited warning in console which results in divergence
- what could be the possible problem, any solution pls reply.

[AlexanderZ]

compile your code and read log -> this is what you better do each time after any modification

code seems ok, however you may do following:
1. change double to real type
2. put #define {a_0,a_1,...} part into the header it's not supposed to be inside the loop
3. of course you shold check you equation and parametes it should be smooth and defined on the whole requested range.
4. As you have rather tight range of 2.1 degree make sure you have small enougth timestep to simulate it