Program1

From New IAC Wiki
Revision as of 23:31, 24 July 2012 by Starvale (talk | contribs) (Created page with '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load yield.m Y=yield(:,2); % Y(A) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% f…')
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

load yield.m Y=yield(:,2); % Y(A) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


for j=1:61 % Possible masses of the heavy fragment

   A=j+119;%heavy fragment mass

for i=1:61 %Possible numbers of protons and neutron for any given heavy fragment A

   Z(i)=i; %protons
   N(i)=A-Z(i); %neutrons

load table.m % Table of parameters: Atomic mass, four parameters (Dx Fn Sz Fz) and the most probable charge Zp=table(:,6); % last column of the table Sz=table(:,4); % fourth column of the table fn=table(:,3); fz=table(:,5); V(i,j)=(Z(i)-Zp(j)+0.5)./(Sz(j)*sqrt(2)); W(i,j)=(Z(i)-Zp(j)-0.5)./(Sz(j)*sqrt(2)); E(i,j)=erf(V(i,j))-erf(W(i,j));

   if rem(Z(i),2)==0
        if rem(N(i),2)==0
             fa(i,j) = fz(i)*fn(j);
           end
      
   end
   if rem (Z(i),2)==0
                  if rem(N(i),2)>0;
             fa(i,j)= fz(i)/fn(j)
           end
       end
     if rem(Z(i),2)>0
                   if rem(N(i),2)==0
               fa(i,j)= fn(i)/fz(j);
           end
        end
   if rem(Z(i),2)>0
                  if rem(N(i),2)>0
               fa(i,j)= 1/(fz(i)*fn(j));
           end
    
   end
   
       
               

end

% P(Z|A)=1/2*F(A)*[erf(V)-erf(W)] We still need to add the normalization 										factor.

hold on end P=0.5*fa.*E; % Taking the sum. FI=(P)'; % Finishing it up - Y(Z,A) for h=1:61

   S(h)=sum(FI(h,:));
   FIf(h,:)=FI(h,:)./S(h)

end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% load en.m en_av=en(:,2); sigma=en(:,3); for f=1:61 for l=120:240 TKE=l En(f,l) = 1/sqrt(2*pi*sigma(f).^2)*exp((-1./(2*sigma(f).^2))*(TKE-(en_av(f))).^2); end end En(:,1:119)=[]; En2=flipud(En); Enf=[En2;En];