for j = 1:24
    WsAug28(j) = sum(Current.Energy(find(Current.Month == 8 & Current.Year == 2013 & Current.Day == 28 & Current.Hour == j)));
end
kWhAug28 = WsAug28/60/60/1000;
hours = 1:24;
figure(2);
plot(hours,kWhAug28)
xlabel('Hour on August 28th'); ylabel('Hourly Energy Usage (kWh)'); title('Hourly Energy Usage on August 28th, 2013'); axis tight