for i=1:101
    Ei(i)=heat_added(i+1)-heat_added(i);
end

