The arithmetic standard operators +, -, *, /, **, % can be applied to arrays.
The arithmetic operations between equal-sized arrays are applied elementwise.
This is called vectorization.
>>> x = np.array([1,5,2])
>>> y = np.array([1,2,3])
>>> x + y
array([2, 7, 5])
>>> x - y
array([ 0, 3, -1])
>>> x * y
array([ 1, 10, 6])
>>> x / y
array([ 1. , 2.5 , 0.66666667])
>>> x // y
array([1, 2, 0])
>>> x % y
array([0, 1, 2])
>>> x ** 2
array([ 1, 25, 4])
>>> x ** y
array([ 1, 25, 8])