
function k=texgausspdf(t,p)
% EXGAUSSPDF The ex-Gaussian pdf
% Syntax k=exgauss(t,p)
mu=p(1);
 sigma=p(2);
% sigma=36;
tau=p(3);
k=exp(-t./tau + mu./tau + sigma.^2./2./tau.^2).*normcdf((t-mu-sigma.^2./tau)./sigma)./tau;
k(k==Inf)=zeros(length(k(k==Inf)),1);
if tau<0 | sigma<0,
   k=zeros(length(k),1);
end   
end
