clc
//to calculate proton gain in mass
c=3*10^8 //light speed(m/s)
K=500*10^6*1.6*10^-19 //kinetic energy (J)
deltam=K/c^2
disp("proton gain in mass is delm="+string(deltam)+"kg")
//answer is given wrong in the book=8.89*10^28 kg