<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def fermi_dist(energy, beta):
 """ Fermi Dirac distribution"""

    exponent = np.asarray(beta*energy).clip(-600, 600)
    return 1./(np.exp(exponent) + 1)