Anthônio Nunes Moreira Netto - anthonionetto@fisica.ufc.br
Departamento de Física, Centro de Ciências, Universidade Federal do Ceará
In [1]:
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
In [2]:
g = 9.8
dt = 0.5
n = 20
v0 = 0.
tmax = n*dt
v = [];
v.append(v0)
In [3]:
t = np.linspace(0, tmax, n) # Eixo-x (Tempo).
theory = v0-g*t # Solução analítica.
for i in range(1,n):
v.append(v[i-1]-g*dt)
In [4]:
plt.plot(t,theory,'b-',t,v,'ro')
plt.axis([0.,tmax,theory[n-1],0.])
plt.legend(['Teórico','Numérico'],loc=0)
plt.title('Exercício 1.1')
plt.xlabel('Tempo (s)')
plt.ylabel('Velocidade (m/s)')
plt.show()
In [ ]: