Integral Metodo do Trapézio

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

import matplotlib.pyplot as plt
def trap(g, inicio, fim, n):
 h = (fim - inicio) / float(n)
 soma = (f(inicio) + f(fim)) / 2
 for i in range(n):
 x = inicio + i * h
 soma += g(x)
 return soma * h
def f(x):
 return (10*x) / ((x ** 2) + (0.1)**2) ** (0.5)
print("Tx =", trap(f, 0.1, 0.4, 1000))
X=[]
Y=[]
i=0.1
while (i<0.4):
 X.append(i)
 a= ((10*i)/(2.2*((i**2)+(0.1**2))**0.5))-(9.8*0.4)+((10*0.1*0.4)/(2.2*((i**2)+(0.1**2))**0.5))
 Y.append(a)
 i=i+(0.01)
print (X)
print (Y)
plt.plot(X,Y)
plt.xlabel('x')
plt.ylabel('aceleração')
plt.show()