How to solve a system of differential equations?

I'm trying to solve the following set of differential equations. I followed a mathworks examples, however, I receive a solution of "empty sym". I'm ultimately trying to plot X and y as functions of V.

Cao = 10*101325/8.314/400/1000;
k = 10^(-4);
Fao = 2.5;
alpha = .001;

syms X(V) y(V)

Solution = dsolve(diff(X) == k*Cao*(1-X)*y/Fao/(1+2*X), diff(y) == -alpha*(1+2*X)/2/y)

Answer 1

[Walter Roberson]

You do not have much of a chance of solving that for a closed form solution. I tried, and I am still trying to interpret what I got back.

The boundary conditions are important for this; what is known? You have a singularity at y = 0 so a boundary there is not useful.

But if you want to experiment further then you need to be more explicit in your equations:

Solution = dsolve(diff(X(V),V) == k*Cao*(1-X(V))*y(V)/Fao/(1+2*X(V)), diff(y(V),V) == -alpha*(1+2*X(V))/2/y(V))