View source for Program2

load Y.mat;
load FI_norm.mat;
load En.mat;

for b=1:61
    Pr2(b,:)=Y(b)*(En(b,:))
end
Pr4=flipud(Pr2);
Pr4(61,:)=[];
Pr=[Pr4; Pr2];
R1=rot90(Pr);
R2=flipud(R1);
[C,h] = contourf(R2,10);
caxis([0.0001 0.2])
%shading interp
axis([1 121 1 121])
set(gca,'XTick',[1 21 41 61 81 101 121])
set(gca,'XTickLabel',{'60','80','100','120','140','160', '180'})
set(gca,'YTick',[1 21 41 61 81 101 121])
set(gca,'YTickLabel',{'120','140','160','180','200','220', '240'})
xlabel('Atomic mass, amu','FontSize',16)
ylabel('Total kinetic energy, MeV','FontSize',16)
title('Fragment distribution for photofission of U-238','FontSize',20)