

%Calculate for 55Fe
Z=26;
A=55;
M=51.1618; %GeV/c^2
[Eb, Weizaecker]=CalcEb(Z,A,M);
%results:
% binding energy Eb is 480.8 MeV
% This is an even-odd or odd-even nucleus. Weizaecker binding energy is 478.6 MeV
% Wiezaecher is off by  0.46 percent

%Calculate for 57Co
Z=27;
A=57;
M=53.0225; %GeV/c^2
[Eb, Weizaecker]=CalcEb(Z,A,M);
%results:
% binding energy Eb is 497.9 MeV
% This is an even-odd or odd-even nucleus. Weizaecker binding energy is 495.7 MeV
% Wiezaecher is off by  0.44 percent

%Calculate for 58Ni
Z=28;
A=58;
M=53.9526; %GeV/c^2
[Eb, Weizaecker]=CalcEb(Z,A,M);
%results:
% binding energy Eb is 506.1 MeV
% This is an even-even nucleus. Weizaecker binding energy is 502.6 MeV
% Wiezaecher is off by  0.68 percent





