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)