## Abstract

A distribution G on (0, ∞) is called matrix-exponential if the density has the form αe^{Tz}s where a is a row vector, T a square matrix and s a column vector. Equivalently, the Laplace transform is rational. For such distributions, we develop an operator calculus, where the key step is manipulation of analytic functions f(z) extended to matrix arguments. The technique is illustrated via an inventory model moving according to a reflected Brownian motion with negative drift, such that an order of size Q is placed when the stock process down-crosses some level q. Explicit formulas for the stationary density are found under the assumption that the leadtime Z has a matrix-exponential distribution, and involve expressions of the form f(T) where f(z) = √1-2z.

Original language | English |
---|---|

Pages (from-to) | 166-176 |

Number of pages | 11 |

Journal | Mathematics of Operations Research |

Volume | 23 |

Issue number | 1 |

DOIs | |

State | Published - Feb 1998 |

## Keywords

- (s, S) model
- Brownian motion
- Computer algebra
- EOQ model
- Inventory system
- Matrix-exponential distribution
- Operator calculus
- Phase-type distribution
- Stochastic decomposition
- Storage model

## ASJC Scopus subject areas

- Mathematics (all)
- Computer Science Applications
- Management Science and Operations Research