'for' e 'if'

[Antonio José]

Linhas 117-133. Ele conta como 16 equações e não como 4, parece que não esta fazendo o somatório de forma correta.

[Rafael de Pelegrini Soares]

Nas linhas que comentaste tens um for duplo:

for i in [1:NComp]
    
    for j in [1:NComp]

Como NComp = 4 para o teu caso:

PP as Plugin (Type="PP", Components = ["hydrogen", "nitrogen", "carbon dioxide", "carbon monoxide"]

Então são esperadas mesmo 16 equações e não 8. Loop duplo em 4 elementos, 4*4 = 16.

Acredito que a equação que desejavas escrever é:

for i in [1:NComp]
	Dm(i)=(1-z(i))/(sum(z(:)/Dj(i,:)));
end

Daí terás 4 equações no lugar de 16, acho que era o que querias.