Why is it not working?

2 vues (au cours des 30 derniers jours)
Rebecca D'Onofrio
Rebecca D'Onofrio le 7 Sep 2021
Modifié(e) : Jan le 7 Sep 2021
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 commentaires
Mathieu NOE
Mathieu NOE le 7 Sep 2021
Could not run the code
at least K and jj not defined
Mathieu NOE
Mathieu NOE le 7 Sep 2021
well , most of the variables used in the main loop are not defined / initialized...

Connectez-vous pour commenter.

Réponse acceptée

Simon Chan
Simon Chan le 7 Sep 2021
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

Plus de réponses (0)

Catégories

En savoir plus sur Matrix Indexing dans Help Center et File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by