def get_coordinates(filename):
    """ Extract information about atoms and their coordinates from .xyz file."""
    with open(filename, "r") as f:
        data = []
        lines = f.readlines()
        for line in lines:
            data.append(line.split())
    N = int(data[0][0])
    data = data[2:]
    atom_types, coords = [], []
    for x in range(N):
        atom_types.append(data[x][0])
        coords.append([float(data[x][1]), float(data[x][2]), float(data[x][3])])
    return atom_types, coords