docs: def initial_finall_mass_relation(self,marker='o',linestyle='--'):
  		'''
  			INtiial to final mass relation
  		'''
  
  		final_m=[]
                  ini_m=[]
                  for i in range(len(self.runs_H5_surf)):
  			sefiles=se(self.runs_H5_out[i])
                          ini_m.append(sefiles.get("mini")) 
  			h1=sefiles.get(int(sefiles.se.cycles[-2]),'H-1')
  			mass=sefiles.get(int(sefiles.se.cycles[-2]),'mass')
  			idx=-1
  			for k in range(len(h1)):
  				if h1[k]>0.1:
  					idx=k
  					break
  			final_m.append(mass[idx])
  		label='Z='+str(sefiles.get('zini'))
  		plt.plot(ini_m,final_m,label=label,marker=marker,linestyle=linestyle)
  		plt.xlabel('$M_{Initial} [M_{\odot}]$',size=23)
  		plt.ylabel('$M_{Final} [M_{\odot}]$',size=23)