Buenas noches.
No logramos obtener los resultados esperados de ningún modo, por lo que vamos a poner en duda cosas que antes obviábamos...
El c_e (ya sea OD_e o DBO_e) ¿es el mismo para cada paso de la iteración?, es decir, ¿es siempre el c_e que nos dan cómo dato? Porque pensaba que quizás el dato de c_e sea solo para el paso inicial, y luego sería el valor del contaminante en el paso n-1.
Ejemplifico la idea... Si estoy calculando OD por Euler explícito, sería:
OD[n+1] = OD[n] + h*f(OD_e=OD_e(dato)) ---> en el paso 0
OD[n+1] = OD[n] + h*f(OD_e=OD[n-1]) ---> en los pasos siguientes
Obviamente f dependería de todas las demás variables del paso n (OD[n], DBO[n], V[n], etc)