Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def f2(x):
- return (x**2) * np.sin(x)
- def bisseccao(f,a,b, niter=1000,minimo = 0.000001, proximidadeIter = 0.000001):
- raiz = 0
- anterior = 0
- while(a<b):
- meio = (a+b)/2
- anterior = abs(anterior-meio)
- if((abs(a-b) < minimo) or (anterior < proximidadeIter) or (niter == 0)):
- return meio
- if(f(meio) * f(a) < 0):
- b = meio
- else:
- a = meio
- niter = niter - 1
- return raiz,niter
- def secantes(f,a,b, niter=1000,minimo = 0.000001, proximidadeIter = 0.000001):
- raiz = 1
- anterior = 0
- cont = 0
- while(cont < niter):
- if(abs(raiz-anterior) < minimo or f(raiz) < minimo):
- return raiz,cont
- anterior = raiz
- raiz = b - f(b) / (f(b)-f(a)**(b-a))
- a = b
- b = raiz
- cont+=1
- niter = cont
- return raiz,niter
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement