def write_axsf(ofile, c):
    """ Write the Crys in .axsf file format"""
    print('ANIMSTEPS 1', file=ofile)
    print('CRYSTAL', file=ofile)
    print('PRIMVEC 1', file=ofile)
    for k in range(3):
        print('     %8.5f %8.5f %8.5f' % (c.matrix[k][0], c.matrix[k][1], c
            .matrix[k][2]), file=ofile)
    print('PRIMCOORD 1', file=ofile)
    print('%d 1' % c.natom, file=ofile)
    for i in range(c.natom):
        print('%3s %8.3f %8.3f %8.3f ' % (c.atom_element[i], c.atom_xyz[i][
            0], c.atom_xyz[i][1], c.atom_xyz[i][2]), file=ofile)
    return
