Why is it not working?

조회 수: 3(최근 30일)
Rebecca D'Onofrio
Rebecca D'Onofrio 2021년 9월 7일
편집: Jan 2021년 9월 7일
I wrote this code, and it does not give any errors, but still I think it is not working since I cannot see any results! Is there any problem?
B2=6E-5;
n=0;
while abs(K-1) > 10^-4
A=zeros(26,26);
for ii=1:26
if ii==jj
A(ii,jj)=D(ii).*B2+(CAPTURE(ii)+FISSION(ii)+ELREM(ii)+INELASTIC(ii));
end
end
A(ii,jj)=A(ii,jj)-SCATTERING(ii,jj);
FLUSSI=A^-1.*MATFIS;
SOURCE=0;
SOURCE=SOURCE+(F39.*FLUSSI+F40.*FLUSSI);
K=SOURCE;
deltaB=10^-6;
B2=B2+deltaB;
n=n+1;
end
R=pi/(sqrt(B2))
  댓글 수: 2
Mathieu NOE
Mathieu NOE 2021년 9월 7일
well , most of the variables used in the main loop are not defined / initialized...

댓글을 달려면 로그인하십시오.

채택된 답변

Simon Chan
Simon Chan 2021년 9월 7일
Index jj is not defined inside the for-loop, i guess the code should be modified as the following
for ii=1:26
for jj = 1:26
if ii==jj
A(ii,jj)=D(ii).*B2+(CAPTURE(ii)+FISSION(ii)+ELREM(ii)+INELASTIC(ii));
end
end
end

추가 답변(0개)

태그

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by