Advertisement
Guest User

asdsad

a guest
Apr 19th, 2018
68
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.88 KB | None | 0 0
  1. def f2(x):
  2. return (x**2) * np.sin(x)
  3.  
  4. def bisseccao(f,a,b, niter=1000,minimo = 0.000001, proximidadeIter = 0.000001):
  5. raiz = 0
  6. anterior = 0
  7. while(a<b):
  8. meio = (a+b)/2
  9. anterior = abs(anterior-meio)
  10. if((abs(a-b) < minimo) or (anterior < proximidadeIter) or (niter == 0)):
  11. return meio
  12. if(f(meio) * f(a) < 0):
  13. b = meio
  14. else:
  15. a = meio
  16. niter = niter - 1
  17. return raiz,niter
  18.  
  19. def secantes(f,a,b, niter=1000,minimo = 0.000001, proximidadeIter = 0.000001):
  20. raiz = 1
  21. anterior = 0
  22. cont = 0
  23. while(cont < niter):
  24. if(abs(raiz-anterior) < minimo or f(raiz) < minimo):
  25. return raiz,cont
  26. anterior = raiz
  27. raiz = b - f(b) / (f(b)-f(a)**(b-a))
  28. a = b
  29. b = raiz
  30. cont+=1
  31. niter = cont
  32. return raiz,niter
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement