That is difficult to read. It appears to be this:
p = @(delta,sigma) (exp(-sigma./delta)./delta) .* (sigma >= 0);
It will default to being 0 for sigma < 0, so it is not necessary to specifically include that condition.
Call the ‘p’ anonymous function with the two arguments as you would any other function.
