import numpy as np

r = np.array([1, 2, 3])
V = (4*np.pi/3)*(r**3)

print('r [m] =', r)
print('V [m^3] =', V)