Transform-based algorithms have wide applications in applied probability, but rarely provide computable error bounds to guarantee the accuracy. We propose an inversion algorithm for two-sided Laplace transforms with computable error bounds. The algorithm involves a discretization parameter C and a truncation parameter N. By choosing C and N using the error bounds, the algorithm can achieve any desired accuracy. In many cases, the bounds decay exponentially, leading to fast computation. Therefore, the algorithm is especially suitable to provide benchmarks. Examples from financial engineering, including valuation of cumulative distribution functions of asset returns and pricing of European and exotic options, show that our algorithm is fast and easy to implement.

In this paper we consider Lévy processes without negative jumps, reflected at the origin. Feedback information about the level of the Lévy process (‘workload level’) may lead to adaptation of the Lévy exponent. Examples of such models are queueing models in which the service speed or customer arrival rate changes depending on the workload level, and dam models in which the release rate depends on the buffer content. We first consider a class of models where information about the workload level is continuously available. In particular, we consider dam processes with a two-step release rule and M/G/1 queues in which the arrival rate, service speed, and/or jump size distribution may be adapted depending on whether the workload is above or below some level K. Secondly, we consider a class of models in which the workload can only be observed at Poisson instants. At these Poisson instants, the Lévy exponent may be adapted based on the amount of work present. For both classes of models, we determine the steady-state workload distribution.