def square_root (n,d)
  g1 = (n* 1.0) / 2
  g2 = (g1 + (n /g1)) /2
  while (g1 - g2).abs >= d
    g1 = g2
    g2 = (g1 + (n / g1)) / 2
  end
  g2
end