delta.m

水平集分割算法的C-V模型

function Delta_h = Delta(phi, epsilon)
%   Delta(phi, epsilon) compute the smooth Dirac function
Delta_h=(epsilon/pi)./(epsilon^2+ phi.^2);