I wrote a little demo to compute and render Diffusion Curves using matlab.
The demo lets the user draw a set of curves,
meshes the curves with a properly slitted and constrained Delaunay triangulation
and then applies random colors varying along each of the curves to compute diffusion curves using an FEM solver.
Try it for yourself, unzip the code and issue:
diffusion_curves_demo