                         thickness=thickness,
                                            density=density)
    if energies is None:
        energies = np.linspace(2, 150, 1001)
    energies_keV = energies * u.keV
    absorption = material.absorption(energies_keV)
    return energies, absorption


#def get_transmission(energies):
#    return trans.get_transmission(energies)
if __name__ == &#39;__main__&#39;:
    energies = np.linspace(2, 150, 1001)
    #tr=Transmission()
    #t=tr.get_transmission(energies)
