Mp=1.00728;//mass of proton in amu//
Mn=1.00866;//mass of neutronin amu//
MH=2.01355;//isotopic mass of H atom in amu//
dM=((1*Mp)+(1*Mn)-MH);//dM value of H atom in amu//
printf('dM value of H atom=dM=%famu',dM);
BE=dM*931;//binding energy of H atom in MeV//
printf('\nBinding energy of H atom=BE=%fMeV',BE);