// erro_quad.sci

function E = erro_quad (f, T)
  [nl, nc] = size(T);
  
  E = 0;
  for i = 1:nc do
    E = E + (f(T(1,i)) - T(2,i))^2;
  end
endfunction

