Listing 2

// circuit.c
#include "dvec.h"
#include "dgemat.h"
#include "dludecmp.h"
#include "rstream.h"

enum bool {FALSE, TRUE };

int     main(int argc, char **argv)
   {
   DoubleVec   Voltage;
   DGEMatrix   Impedance;
   cin >> Voltage >> Impedance;    // read voltage vector, impedance matrix
   DoubleVec   Current = solve(Impedance, Voltage); // solve equation
   cout << Current;                // print current vector
   }