I am doing a non-linear thermal model. The heat of fusion is of course a large non-linearity. When doing a dynamic model I can do a crude approximation with a jump in the sensible heat over a small temperature range. That's kind of ugly though. Is there a cleaner way to easily do this, or should I do the crude thing? Thanks, David