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

def mask_average(dset,mask):
    '''Returns average of voxels in ``dset`` within non-zero voxels of ``mask``'''
    o = nl.run(['3dmaskave','-q','-mask',mask,dset])
    if o:
        return float(o.output.split()[-1])