close all
clear all
clc

PI= 150;
PII= 230;
PIII= 190;

% pomerný prírastok palivových nákladov TE:
% bt= 4 + 0.15*Pt + 0.003*Pt^2

% Spotrebná charakteristika vody pre VE:
% V= 10 + 0.2*Pv + 0.03*Pv^2

% celková spotreba vody na VE za reg. obdobie je: [m3] 
V1= 900;


% I interval

PtI_0= 90;
PVI_0= 60;

btI_0= 4 + 0.15*PtI_0 + 0.003*PtI_0^2 %menim podla zadania
vI_0= 0.2 + 2*0.03*PVI_0 %menim po derivacii

alfa= btI_0/vI_0
WI_0= PtI_0 + PVI_0 - PI % < epsilon ---- platí


% II interval

PtII_0= 90;
PVII_0= 60;

WII_0= PtII_0 + PVII_0 - PII % abs(WII_0) > epsilon --- neplatí

% btII_0 = btI_0 .... deltabtII_0 volíme napr 20 
deltabtII_0= 20; %tipovane čislo (nemoze byt vacsie ako -80 vysledok preto tipujeme 20 lebo 40+20 je menej ako 80)
btII_1= btI_0 + deltabtII_0

vII_1= btII_1/alfa

PtIIa= [0.003 0.15 4-btII_1]; %menime podla btI(0)
PtIIa_xy= roots(PtIIa);
PtII_1= PtIIa_xy(2,1)

PVII_1= (vII_1-0.2)/0.06 %menim podla V rovnice a derivujem

WII_1= PtII_1 + PVII_1 - PII % abs(WII_1) > epsilon --- neplatí

deltabtII_1= (WII_1/(WII_0-WII_1))*deltabtII_0

btII_2= btII_1 + deltabtII_1
PtIIb= [0.003 0.15 4-btII_2] %tu sa mení kod Bta Btb Btc (bt)
PtIIb_xy= roots(PtIIb)
PtII_2= PtIIb_xy(2,1)

vII_2= btII_2/alfa
PVII_2= (vII_2-0.2)/0.06 %menim podla V a derivujem

WII_2= PtII_2 + PVII_2 - PII % abs(WII_1) > epsilon --- neplatí

deltabtII_3= (WII_2/(WII_1-WII_2))*deltabtII_1

btII_3= btII_2 + deltabtII_3
PtIIc= [0.003 0.15 4-btII_3] %tu sa mení kod btI(0)
PtIIc_xy= roots(PtIIc)
PtII_3= PtIIc_xy(2,1)

vII_3= btII_3/alfa
PVII_3= (vII_3-0.2)/0.06 %podla V derivujem

WII_3= PtII_3 + PVII_3 - PII % abs(WII_1) > epsilon --- neplatí


% III interval