// Scilab Code Ex5.1 :: Page-5.2 (2009)
clc;clear;
m_p = 1.007826;     // Mass of a proton, amu
m_n = 1.008665;     // Mass of a neutron, amu
M_He = 4.002604;    // Measured mass of He nucleuc, amu
delta_m = 2*m_p+2*m_n - M_He;   // Mass defect of He, amu
printf("\nThe mass defect of He = %f amu", delta_m);

// Result
//  The mass defect of He = 0.030378 amu