questão 16

Prévia do material em texto

23/10/2020 questão 16
localhost:8888/nbconvert/html/Desktop/questão 16.ipynb?download=false 1/2
In [ ]:
######## questão 16 ############ 
### Método de Newton ---- Implementação python 
In [ ]:
from __future__ import division 
import numpy as np 
import matplotlib.pyplot as plt 
import scipy as sci 
import sympy 
from sympy import* 
import math 
def f(x): 
 
 return 3**(3*x+1)-7*(5**(2*x)) 
def df(x): 
 h= np.exp(-19)# o -19 é para diminuir o erro da raiz ao calcular 
 df= (f(x+h)- f(x))/h 
 return df 
def newton (x0,Erro, itmax): 
 it=1 
 Er=1 
 x=x0 
 while(Er>=Erro and it< itmax): 
 xold=x 
 x= x-f(x)/df(x) 
 Er= np.abs(float((x-xold)/x)) 
 it= it+1 
 return (x, Er, it) 
23/10/2020 questão 16
localhost:8888/nbconvert/html/Desktop/questão 16.ipynb?download=false 2/2
In [ ]:
### Fazer a inspeção gráfica da função 
""" 
xx= np.linspace (45,54) 
plt.plot (xx, f(xx)) 
plt.grid();plt.show () 
""" 
#####Fazer a inspeção gráfica da derivada da função 
"""yy= np.linspace (-25,10) 
plt.plot (yy, df(yy)) 
plt.grid();plt.show () 
""" 
x0=1 
Erro= 10 **-6 
itmax=800 
res = newton(x0,Erro,itmax) 
print(" Pelo Método da Falsa Posição obtemos:") 
print("a)O valor da raíz é: ", res[0]) 
print("b)O erro relativo foi:", res[1]) 
print("c) O número de iterações:", res[2])