% matrix of [shunt volt mean, shunt volt std, supply volt, supply std]
shunt=1.020;  %ohms
data = [
    80.599e-3   49e-6     4.0162  2e-3      % Empty FPGA
    89.47e-3    36e-6     4.0243  2.214e-3  % Empty socket test
    89.849e-3   33.1e-6   4.026   6e-3      % OISC8 mult 16bit loop
    89.968e-3   35.6e-6   4.0222  2.47e-3   % RISC8 mult 16bit loop
];

I=data(:,1)*shunt;  % current vector
P=(data(:,3)-data(:,1)).*I;  % power in W
Pstd=data(:,2).*data(:,4);

bar(1:4,P)                
hold on
er = errorbar(1:4,P,-Pstd./2,+Pstd./2);