Ecuacion de continuidad

import math
import numpy as np
import matplotlib.pyplot as plt

r2=0.01
A2=math.pi*r2*r2
v2=10

N=10
A1=np.zeros(N)
v1=np.zeros(N)
v2=10

aux=0.2
for i in range(0,N):
  A1[i]=aux*A2
  v1[i]=(v2*A2)/A1[i]
  aux=aux+0.2

plt.plot(A1/A2,v1/v2,'-o',label='Teoria')
plt.title("Ecuacion de continuidad")
plt.xlabel("A1/A2")
plt.ylabel("v1/v2")
plt.grid()
plt.show()