"There is already a "delta" function defined in the laplace package" I will have to take a look at the other dirac function but I believe that it is only recognised by laplace. I didn't know there were two unit_step's, I will look into that as well. Also, thanks for your support of my work. Rich