The development of a program based on 2 × 2 matrices for computing the linear steady-state performance of networks against chosen frequencies of a sinusoidal input signal is described. Limitations and projected extensions of the program are briefly discussed.