function g = deltah2(q)
global m n rho mu grav
global Rey ff

Di=[0,0.4,0.3,0.3,0.25,0.3,0.3,0.25,0.25,0.3,0,...
    0.35,0.3,0.25,0.2,0.25,0.2,0.25,0.2,0.25,0.2];
L=1000*[0,3.5,1.5,2,1.5,2,3,1.5,1.5,1.5,0,0.5,1.8,...
  1.5,2.2,1.5,2.2,1,1.5,1,1.5];
 
eps=2e-4*ones(1,m);
rho=997;mu=0.890e-3;grav=9.80665;
qb1 = [0.15 0.25 0.35]; % [m^3/s] vazoes da bomba 1
hb1 = [50 47 42]; % [m] alturas da bomba 1
b1 = polyfit(qb1,hb1,2); % curva da bomba 1
qb2 = [0.20 0.25 0.30]; % [m^3/s] vazoes da bomba 2
hb2 = [30 28 25]; % [m] alturas da bomba 2
b2 = polyfit(qb2,hb2,2); % curva da bomba 2
Rey=[];ff=[];
for i=1:m
 if (i!=1 && i!=11) 
  u=q(i)/(pi*Di(i)*Di(i)/4);
  Rey(i,1)=abs(u)*Di(i)*rho/mu;
  if (Rey(i,1) < 100)
   Re=100
  endif
  ff(i,1)=friction(Rey(i,1),eps(i)/Di(i));
  g(i)=-ff(i,1)*L(i)/Di(i)*abs(u)*u/(2*grav);
 elseif (i==1)
  g(i)=polyval(b1,q(i)); ff(i)=0; Rey(i,1)=0;
 elseif (i==11)
  g(i)=polyval(b2,q(i)); ff(i)=0; Rey(i,1)=0;
 else
  i
  stop
 endif
endfor
endfunction