ODE45 in Maxwell-Stefan equation

I am trying to solve Maxwell Stefan's equation over a membrane to get the transient mole fraction distribution over the membrane thickness 'z'. But somehow I am not able to code it using ODE45, more likely I am not able to write the system to solve using ODE45. It will be really great if someone can help me with the primary syntaxes and function. The equation I am trying to solve is

where c is concentration and is binary diffusion coefficient and z varies from 0 to 2.5e-4 cm.

Thank you in advance.

[darova]

[Anantrao Vijay Shirsath]

Hello, @darova, Sorry for the mistake, I have corrected it. Cheers!

Answer 1

[darova]

Here is the solution:

f = @(t,y) 1/c/D*(y*(NH2+NH2O)-NH20);
[t,y] = ode45(f,[0 2.5e-4],1);

[Anantrao Vijay Shirsath]

The only problem in this anwer is that the value of y should decrease over z, but here it increases!

[darova]

impossible. can you show?