import numpy as np
import matplotlib.pyplot as plt

#Nuclei names
nuclei= ['H-2','He-3', 'He-4', 'Li-6', 'C-12', 'Kr-84', 'Sn-119', 'Ti-205', 'U-235', 'U-238']

#Atomic mass values for each nuclei in atomic mass unit
#Source: Atomic Mass Data Center (AMDC)

atomic_masses=[2,3,4,6,12,84,119,205,235,238]

#Binding energies in MeV from AMDC
binding_energies= [2.224,7.718,28.296,31.994,92.162,732.265,1011.434,1615.071,1778.567,1801.690]

#Calculate E/A values 
E_A = np.array(binding_energies)/np.array(atomic_masses)

#Plot the graph
plt.figure(figsize=(10,6))
plt.scatter(range(1,len(nuclei)+1), E_A, color ='skyblue', label='Data')
plt.plot(range(1,len(nuclei)+1), E_A, color='pink', label='Line')

plt.xticks(range(1,len(nuclei)+1), nuclei, rotation='vertical')
plt.xlabel('Mass number, A')
plt.ylabel('Binding Energy per Nucleon, E/A (MeV)')
plt.title('Binding Energy per Nucleon for Different Atomic Nuclei')
plt.legend()
plt.tight_layout()
plt.show()