Hi, I also have to do calculation with the Laplace transform involving the heaviside function. Is it possible to have your code for the lap and ilap functions please? That would be very useful... Thanks a lot. Ben