Factor de efectividad fenomenos

Encontrar el factor de efectividad no isotérmico de una reacción de tipo
2A B
BALANCE DE MATERIA
dNardr + 2*Narr + K01*e^-E/RT *( Ca2- 1Kc*(CT-Ca)) =0
Ni-Ni-1r + 2*Nir + K01*e^-E/RT *( Ci2- 1Kc*(CT-Ci)) =0

f(N1) = ((N1 - N0) / dr) + 2 * N1 / (1 * dr) + K01 * exp(-9700 / ((T1 - 32) * 5 / 9 + 460)) * (C1 ^ 2 - 1 * (CT - C1) / (KC))
N1(0) = 1
f(N2) = ((N2 - N1) / dr) + 2* N2 / (2 * dr) + K01 * exp(-9700 / ((T2 - 32) * 5 / 9 + 460)) * (C2 ^ 2 - 1 * (CT - C2) / (KC))
N2(0) = 1
f(N3) = ((N3 - N2) / dr) + 2 * N3 / (3 * dr) + K01 * exp(-9700 / ((T3 - 32) * 5 / 9 + 460)) * (C3 ^ 2 - 1 * (CT - C3) / (KC))
N3(0) = 1
f(N4) = ((N4 - N3) / dr) + 2 * N4 / (4 * dr) + K01 * exp(-9700 / ((T4 - 32) * 5 / 9 + 460)) * (C4 ^ 2 - 1 * (CT - C4) / (KC))
N4(0) = 1
f(N5) = ((N5 -N4) / dr) + 2 * N5 / (5 * dr) + K01 * exp(-9700 / ((T5 - 32) * 5 / 9 + 460)) * (C5 ^ 2 - 1 * (CT - C5) / (KC))
N5(0) = 1
f(N6) = ((N6 - N5) / dr) + 2 * N6 / (6 * dr) + K01 * exp(-9700 / ((T6 - 32) * 5 / 9 + 460)) * (C6 ^ 2 - 1 * (CT - C6) / (KC))
N6(0) = 1
f(N7) = ((N7 - N6) / dr) + 2 * N7 / (7 * dr) + K01 * exp(-9700 / ((T7 - 32) * 5 / 9 + 460)) * (C7 ^ 2 - 1 * (CT - C7) / (KC))
N7(0) = 1f(N8) = ((N8 - N7) / dr) + 2 * N8 / (8 * dr) + K01 * exp(-9700 / ((T8 - 32) * 5 / 9 + 460)) * (C8 ^ 2 - 1 * (CT - C8) / (KC))
N8(0) = 1
f(N9) = ((N9 - N8) / dr) + 2 * N9 / (9 * dr) + K01 * exp(-9700 / ((T9 - 32) * 5 / 9 + 460)) * (C9 ^ 2 - 1 * (CT - C9) / (KC))
N9(0) = 1
f(N10) = ((N10 - N9) / dr) + 2 * N10 / (10 * dr) + K01 * exp(-9700 / ((T10 - 32) * 5 / 9 + 460)) * (C1 ^ 2 - 1 * (CT - C10) /(KC))
N10(0) = 1

N0 = 0

LEY DE FICK
dCadr+ (1- Ca2*CT )* NarDab = 0
Ci-Ci-1▲▲ ▲r +(1- Ca2*CT )* NiDab=0
f(C0) = (C1 - C0) / dr + (1 - C1 / (2 * CT)) * N1 / DAB
C0(0) = 10
f(C1) = (C2 - C1) / dr + (1 - C2 / (2 * CT)) * N2 / DAB
C1(0) = 10
f(C2) = (C3 - C0) / dr + (1 - C3 / (2 * CT)) * N3 / DAB
C2(0) = 10
f(C3) = (C4 - C0) / dr + (1 - C4 / (2 * CT)) * N4 / DAB
C3(0) = 10
f(C4) =(C5 - C0) / dr + (1 - C5 / (2 * CT)) * N5 / DAB
C4(0) = 10
f(C5) = (C6 - C0) / dr + (1 - C6 / (2 * CT)) * N6 / DAB
C5(0) = 10
f(C6) = (C7 - C0) / dr + (1 - C7 / (2 * CT)) * N7 / DAB
C6(0) = 10
f(C7) = (C8 - C0) / dr + (1 - C8 / (2 * CT)) * N8 / DAB
C7(0) = 10
f(C8) = (C9 - C8) / dr + (1 - C9 / (2 * CT)) * N9 / DAB
C8(0) = 10
f(C9) = (C10 - C9) / dr + (1 - C10 / (2 * CT)) * N10 / DABC9(0) = 10
C10 = CT - CB0

BALANCE DE ENERGIA
dQrdr + 2*Qrr – ▲Hr*e^-E/RT *( Ca2- 1K*(CT-Ca)) =0
Qi-Qi-1dr - 2*K*(Ti-Ti-1)idr^2 – ▲Hr*e^-E/RTi *( Ci2- 1KC*(CT-Ci)) =0
f(Q1) = (Q1 - Q0) / dr - (2 * K / (1 * dr)) * (T1 - T0) / dr + Hr * K01 * exp(-9700 / ((T1 - 32) * 5 / 9 + 460)) * (C1 ^ 2 - 1 * (CT - C1) / (KC))
Q1(0) = 10
f(Q2) = (Q2 - Q1) / dr - (2 * K / (2 * dr)) * (T1 - T0) / dr + Hr * K01* exp(-9700 / ((T2 - 32) * 5 / 9 + 460)) * (C2 ^ 2 - 1 * (CT - C2) / (KC))
Q2(0) = 10
f(Q3) = (Q3 - Q2) / dr - (2 * K / (3 * dr)) * (T1 - T0) / dr + Hr * K01 * exp(-9700 / ((T3 - 32) * 5 / 9 + 460)) * (C3 ^ 2 - 1 * (CT - C3) / (KC))
Q3(0) = 10
f(Q4) = (Q4 - Q3) / dr - (2 * K / (4 * dr)) * (T1 - T0) / dr + Hr * K01 * exp(-9700 / ((T4 - 32) * 5 / 9 + 460)) * (C4 ^ 2 - 1 * (CT - C4) / (KC))Q4(0) = 10
f(Q5) = (Q5 - Q4) / dr - (2 * K / (5 * dr)) * (T1 - T0) / dr + Hr * K01 * exp(-9700 / ((T5 - 32) * 5 / 9 + 460)) * (C5 ^ 2 - 1 * (CT - C5) / (KC))
Q5(0) = 10
f(Q6) = (Q6 - Q5) / dr - (2 * K / (6 * dr)) * (T1 - T0) / dr + Hr * K01 * exp(-9700 / ((T6 - 32) * 5 / 9 + 460)) * (C6 ^ 2 - 1 * (CT - C6) / (KC))
Q6(0) = 10
f(Q7) = (Q7 - Q6) / dr - (2 * K / (7 * dr)) * (T1 - T0) / dr + Hr * K01* exp(-9700 / ((T7 - 32) * 5 / 9 + 460)) * (C7 ^ 2 - 1 * (CT - C7) / (KC))
Q7(0) = 10
f(Q8) = (Q8 - Q7) / dr - (2 * K / (8 * dr)) * (T1 - T0) / dr + Hr * K01 * exp(-9700 / ((T8 - 32) * 5 / 9 + 460)) * (C8 ^ 2 - 1 * (CT - C8) / (KC))
Q8(0) = 10
f(Q9) = (Q9 - Q8) / dr - (2 * K / (9 * dr)) * (T1 - T0) / dr + Hr * K01 * exp(-9700 / ((T9 - 32) * 5 / 9 + 460)) * (C9 ^ 2 - 1 * (CT - C9) / (KC))...