clear ;

in = 1;
for n=20:20:80
    
  % Resolver um bairro quadrado de dimensoes cada vez maiores
  m = n;
  CH = 1.0;
  CV = 2.0;
  [nv,nc,conec,conducs,coords] = RedeHidraBairro(n,m,CH,CV);
  
  tic ;
  A = Assembly(nv, nc, conec, conducs);
  QB = 10;
  nB = 1;
  natm = nv;
  Pref = 0.0;
  [Atilde b] = BuildSystem(nv,A,QB,nB,Pref,natm);
  tass(in) = toc;

  tic ;
  p = Atilde \ b;
  dims(in) = nv;
  tsol(in) = toc;
  in = in + 1;

end

loglog(dims,tsol, 'b-o');