www.pudn.com > levelsetImage.rar > 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 ) );