dirac.m

Image segmentation script (code)

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 ) );