

function Thetas = initializeThetas( sizes, EPSILON = 0.03 )

   L = length(sizes);  % nr de camadas

   for i = 1:L-1,
      Thetas{i} = rand(sizes(i+1), sizes(i)+1) * 2 * EPSILON - EPSILON;
   end

end
