code: def plot_images(self, outfile):
          
          sum_struct = self.__images[0].sites
          for image in self.__images:
              for site_i in self.__relax_sites:
                  sum_struct.append(PeriodicSite(image.sites[site_i].specie,
                                                 image.sites[site_i].frac_coords,
                                                 self.__images[0].lattice,
                                                 to_unit_cell=True,
                                                 coords_are_cartesian=False))
          sum_struct = Structure.from_sites(sum_struct, validate_proximity=False)
          p = Poscar(sum_struct)
          p.write_file(outfile)
docs: Generates a POSCAR with the calculated diffusion path with respect to the first endpoint.
          :param outfile: Output file for the POSCAR