dirac.m

这是一个利用水平集方法进行分割的方法

function y = dirac( x, epsilon )
% DIRAC Dirac function of x
%    DIRAC( x, epsilon ) Computes the derivative of the heaviside
%    function of x with respect to x. Regularized based on epsilon.

y = ( 1 ./ pi ) .* ( epsilon ./ ( epsilon.^2 + x.^2 ) );