https://wiki.iac.isu.edu/index.php?action=history&feed=atom&title=Program1 2026-07-10T04:34:03Z Revision history for this page on the wiki MediaWiki 1.35.2 https://wiki.iac.isu.edu/index.php?title=Program1&diff=75780&oldid=prev 2012-07-24T23:31:04Z

<p>Created page with &#039;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load yield.m Y=yield(:,2); % Y(A) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% f…&#039;</p> <p><b>New page</b></p><div>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br /> <br /> load yield.m<br /> Y=yield(:,2); % Y(A)<br /> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br /> <br /> <br /> for j=1:61 % Possible masses of the heavy fragment<br /> A=j+119;%heavy fragment mass<br /> for i=1:61 %Possible numbers of protons and neutron for any given heavy fragment A<br /> Z(i)=i; %protons<br /> N(i)=A-Z(i); %neutrons<br /> load table.m % Table of parameters: Atomic mass, four parameters (Dx Fn Sz Fz) and the most probable charge<br /> Zp=table(:,6); % last column of the table<br /> Sz=table(:,4); % fourth column of the table<br /> fn=table(:,3);<br /> fz=table(:,5);<br /> V(i,j)=(Z(i)-Zp(j)+0.5)./(Sz(j)*sqrt(2));<br /> W(i,j)=(Z(i)-Zp(j)-0.5)./(Sz(j)*sqrt(2));<br /> E(i,j)=erf(V(i,j))-erf(W(i,j));<br /> <br /> if rem(Z(i),2)==0<br /> if rem(N(i),2)==0<br /> fa(i,j) = fz(i)*fn(j);<br /> end<br /> <br /> end<br /> if rem (Z(i),2)==0<br /> if rem(N(i),2)&gt;0;<br /> fa(i,j)= fz(i)/fn(j)<br /> end<br /> end<br /> if rem(Z(i),2)&gt;0<br /> if rem(N(i),2)==0<br /> fa(i,j)= fn(i)/fz(j);<br /> end<br /> end<br /> if rem(Z(i),2)&gt;0<br /> if rem(N(i),2)&gt;0<br /> fa(i,j)= 1/(fz(i)*fn(j));<br /> end<br /> <br /> end<br /> <br /> <br /> <br /> end<br /> % P(Z|A)=1/2*F(A)*[erf(V)-erf(W)] We still need to add the normalization factor.<br /> hold on<br /> end<br /> P=0.5*fa.*E;<br /> % Taking the sum.<br /> FI=(P)'; % Finishing it up - Y(Z,A)<br /> for h=1:61<br /> <br /> S(h)=sum(FI(h,:));<br /> FIf(h,:)=FI(h,:)./S(h)<br /> end<br /> <br /> <br /> <br /> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br /> load en.m<br /> en_av=en(:,2);<br /> sigma=en(:,3);<br /> for f=1:61<br /> for l=120:240<br /> TKE=l<br /> En(f,l) = 1/sqrt(2*pi*sigma(f).^2)*exp((-1./(2*sigma(f).^2))*(TKE-(en_av(f))).^2);<br /> end<br /> end<br /> En(:,1:119)=[];<br /> En2=flipud(En);<br /> Enf=[En2;En];</div>

Starvale