function [] = caoEnergyBudgets(whichAuthor)

if strcomp(whichAuthor,'Lennie')
    %% First he re-calculated some numbers based on the lower 
    %% neuron density in human cortex, and the larger number of
    %% synapses per cell
    EPSP = 1.2e9; % molecules of ATP per postsynaptic potential
    propagation = 9.2e8; % molecules of ATP to propagate a spike
                         % in unmyelinated axon (100mm per cell)
    glutamateCycling = 2.3e4; % mol. ATP per vesicle, to release,
                              % uptake, and recycling
    synapsesPerCell = 17500;
    failureRate = .5;
    gluCycling  = glutamateCycling*synapsesPerCell* ...
        failureRate;  %% total cost of glutamate cycling
    presynapticCa = gluCycling; %% ?? didn't see this specified
    %% Lennie's ATP per spike number
    aggregateSpikeCost = EPSP + propagation + gluCycling;
    %% Then he worked backward from total glucose consumption in
    %% the brain, starting with "vegetative" metabolism - things
    %% that have nothing to do with the Na/K ATPase
    vegetativeMetabolism = 1.5e9;  %% 50% of energy consumption
                                   %% is cut by ouabain
    %% neuron and glial resting potentials are re-scaled from
    %% Atwell's numbers, depending on membrane surface area 
    %% and numbers of cells
    neuronRP = 8.6e8; % mol. ATP per sec.
    gliaRP = 3.1e8; %mol. ATP per sec.
    %% Finally, having calculated vegetative and resting metabolism
    %% and knowing how much energy is consumed, he found that 13% of
    %% the energy can be spent on spiking.  Given the number of 
    %% neurons and the cost per spike, he calculated a total average
    %% spiking rate
    spikeRate = 0.16;  % spikes/s on average
    spikeConsumption = spikeRate*aggregateSpikeCost;
    %% make a pair of pie charts
    subplot(1,2,1)
    pie([propagation EPSP neuronRP gluCycling presynapticCa])
    legend('propagation','EPSPs','restingPotentials','glutamate cycling','presynaptic Ca^+^+')
    title('Lennie')
    subplot(1,2,2)
    pie([spikeConsumption vegetativeMetabolism gliaRP neuronRP])
    legend('spikes','vegetative metabolism','glia resting potential','neuron resting potential')
elseif strcomp(whichAuthor,'Attwell')
    %% Pages and pages of calculations of energy consumption ...
    actionPotentials = .47;
    postsynapticReceptors = .34;
    restingPotentials = .13;
    gluCycling = .03;
    presynapticCa = 0.03;
    pie([actionPotentials postsynapticReceptors restingPotentials ...
            gluCycling presynapticCa])
    legend('propagation','EPSPs','restingPotentials','glutamate cycling','presynaptic Ca^+^+')
    title('Attwell and Laughlin')
else 
    error('not familiar with that energy budget')
end