global P_pr

close all;

E=0;
E_brz=0;
for i=1:length(P_pr)
    if P_pr(i)>0
        E=E+1000*P_pr(i)*t_inc;
    end
    if P_pr(i)<0
        E_brz=E_brz-1000*P_pr(i)*t_inc;
    end
end
E_dis=E-E_brz;

P_E_brz=E_brz/E*100;
P_E_dis=E_dis/E*100;

disp('E=celkova energie dodana stroji za 1 cyklus')
disp(E)
disp('E_brz=celkova energie zmarena skrz brzdeni')
disp(E_brz)
disp('E_dis=celkova energie disipovana skrz odpor soustavy')
disp(E_dis)
disp('Podil brzdne energie: P_E_brz=')
disp(P_E_brz)
disp('%')
disp('Podil energie disipovane ztratami soustavy')
disp(P_E_dis)
disp('%')