cylindrical-coordinates code snippets in matlab

solve liquid phase mass balance and energy on a spatially 2d domain in cylindrical coordinates in matlab