Solving quasi-linear heat diffusion using DiffEqTools.jl and DifferentialEquations.jl

∂ut isn’t a function. The f should be a function.