def get_density(element):
return pt.elements.isotope(element).density<EOL>
return the density (g.cm-3) of the element

    Paramters:
    ==========
    element: string. Name of element

    Returns:
    ========
    the density of the element in g.cm-3 units