biseccao.rb

def sign(x)
  return '+' if x >= 0
  return '-'
end

# Calcula f(x)
def f(x)
  x ** 3 + x - 1
end

# Precisao
e = 0.00000000000005

# Valores dos extremos do intervalo
a = 0.0
b = 1.0

# Numero da iteracao
k = 0

# Numero de iteracoes
i = (Math.log(b-a)/Math.log(2) - Math.log(e)/Math.log(2)).ceil

puts "Iterações: #{i}"

sleep(1)

puts "k\ta_i\tf(a_i)\tc_i\tf(c_i)\tb_i\tf(b_i)"

while (b - a).abs > e
  c = (a + b)/2.0

  puts "#{k}\t#{a.round(4)}\t#{sign(f(a))}\t#{c.round(4)}\t#{sign(f(c))}\t#{b.round(4)}\t#{sign(f(b))}"
  # sleep(0.25)

  if f(c) == 0
    break
  end

  if f(a) * f(c) < 0
    b = c
  else
    a = c
  end

  k += 1
end

raiz = (a + b)/2

puts "Raiz aproximada: #{raiz}"