Functionality for efficient computation of properties of system lifetimes was developed, based on the Mathematica framework. The model of these systems consists of a system structure and the components independent lifetime distributions. The components are assumed to be non-repairable. In this work a very general implementation was created, allowing a large number of lifetime distributions from Mathematica for all the component distributions. All system structures with a monotone increasing structure function can be used. Special effort has been made to compute fast results when using the exponential distribution for component distributions. Standby systems have also been modeled in similar generality. Both warm and cold standby components are supported. During development, a large collection of examples were also used to test functionality and efficiency. A number of these examples are presented. The implementation was evaluated on large real world system examples, and was found to be efficient. New results are presented for standby systems, especially for the case of mixed warm and cold standby components.