<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:

def dipole_moment(r_array, charge_array):
    '''Return the dipole moment of a neutral system.
    '''
    return np.sum(r_array * charge_array[:, np.newaxis], axis=0)