function [ prob ] = exponential( u, t )
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here

if (t>=0)
    prob = u * exp(-u * t); 
else
    disp('t<0')
end
end