Matlab code for Chan and Vese paper