Exercício 1.2

Anthônio Nunes Moreira Netto - anthonionetto@fisica.ufc.br
Departamento de Física, Centro de Ciências, Universidade Federal do Ceará

Módulos



In [19]:

    
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np

Variáveis



In [20]:

    
v = 40.
dt = 0.5
n = 40
x0 = 0.
tmax = n*dt
x = [];
x.append(x0)

Cálculo

Analiticamente:
$x(t)=x_0 + v\,t$

Numericamente:
$x_i=x_{i-1}+v\,\Delta{t}$
$t_i=t_{i-1}+\Delta{t}$



In [21]:

    
t = np.linspace(0, tmax, n)  # Eixo-x (Tempo).
theory = x0+v*t   # Solução analítica.

for i in range(1,n):
    x.append(x[i-1]+v*dt)

Resultados



In [22]:

    
plt.plot(t,theory,'b-',t,x,'ro')
plt.axis([0.,tmax,0.,theory[n-1]])
plt.legend(['Teórico','Numérico'],loc=0)
plt.title('Exercício 1.2')
plt.xlabel('Tempo (s)')
plt.ylabel('Posição (m)')
plt.show()



In [ ]: