clear
clc
disp("example 8.4")
flow=[1500 1000 500 500 500 1200 2900 2900 1000 400 600 1600]
cod=1000//constant demand
[m n]=size(flow)
mf(1)=1500
for i=2:n
    mf(i)=mf(i-1)+flow(i)
end
plot(mf)
dd=1:cod:mf(n)
avg=sum(flow)/12
if cod<avg then
    for x=1:6
       t=flow(1,x)
        if t>cod|t==avg then
        
        t=0
    else
        t=cod-t
       end
     flow1(1,x)=t;
 end
 
 else
     for x=1:12
       t=flow(1,x)
    a=cod
        if t>a|t==avg then
        t=0
    else
        t=t-cod
       end
     flow1(1,x)=t;
    end
end

sto=sum(flow1)
printf("storage capacity of plant is %dsec-m-month",sto)