How can I solve ODE in matrix form?

syms y(x)

DE1=Dy==x*(y-1)/(x^2-1);

DE2=Dy+(x+1)*y^2==0;

DE3=Dy+y/x==0;

DE4=Dy-y*exp(x)==0;

DE5=exp(x)*Dy==exp(-y)*(1+exp(x));

DE6=Dy+2*y-x*y==0;

eqns=[DE1 DE2 DE3 DE4 DE5 DE6];

T=feval(symengine,'symobj::dsolve',eqns,x)

I want to find What are DE_n exactly.

but result is 'T =dsolve(matrix([[diff(y(x), x) == (x*(y(x) - 1))/(x^2 - 1), diff(y(x), x) + (x + 1)*y(x)^2 == 0, y(x)/x + diff(y(x), x) == 0, diff(y(x), x) - exp(x)*y(x) == 0, exp(x)*diff(y(x), x) == exp(-y(x))*(exp(x) + 1), 2*y(x) - x*y(x) + diff(y(x), x) == 0]]), x)'

how can I solve ODEs simultaneously by using feval and get each solution?