I am working on a problem that combines Negative Exponential distributions and their impact on the throughput of a system. The output I want to create is a range of Jobs Per Hour given as a probability distribution (and/or 95% confidence interval). The system to be modeled is simple, with two machines connected by a single buffer. Inputs (variables) include cycle time, the Mean Time Before Failure (MTBF) and Mean Time to Repair (MTTR) of each machine, as well as Buffer Size (can be zero) and Time Delay (per part, modeling an index time per position) in the Buffer between the machines. The Negative Exponential's Probability Distribution Function for the MTBF and MTTR values are F(x) = (1/MTBF) x e^((-1/MTBF) x Time). Time is given in minutes. Example, (40/60) minute cycle time (40 seconds), Machine 1 MTBF = 200 minutes / MTTR = 5 minutes, Buffer Size between 0 and 10 units (with a 10 sec delay per part), Machine 2 MTBF = 150 minutes / MTTR = 7 minutes. In this system with no downtime, the throughput is 3600 seconds / 40 second cycle time = 90 Jobs Per Hour Gross. Using the failure inputs for the two machines, the throughput is reduced. By adding a buffer, the throughput recovers. Any insight on a robust approach to model these variables and creating a probability distribution for a range of throughput values would be greatly appreciated.