Lista 03 Questão 7

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

# -*- coding: utf-8 -*-
"""
Created on Sun Oct 08 21:31:06 2017
@author: Vinícius
"""
from math import pi
from raiz import bissecao
def f(x):
 e0 = 8.9*10**(-12)
 F = 1.25
 q = 2*10**(-5)
 Q = 2*10**(-5)
 a = 0.85
 return ((q*Q*x)/((4*pi*e0)*((x**2 + a**2)**1.5))-F)
 
x1 = bissecao(f, 0.0, 1.0, 1E-6, True)
x2 = bissecao(f, 1.0, 2.0, 1E-6, True)
print x1, x2