Hello eveyone,

I'm experiencing a frustrating problem, trying to get maple to solve a system of differential equations. I have defined my system as:

f:=1-2*(n+1)*b*xi^2*nu(xi):

sys:={diff(nu(xi),xi)=(theta(xi)^n - 3*nu(xi)), diff(psi(xi),xi)=(n+1)*b*xi*theta(xi)^n*(1+b*theta(xi))/f, diff(theta(xi),xi)=-(nu(xi)+b*theta(xi)^(n+1))*xi*(1+b*theta(xi))/f, theta(xi=0)=1, nu(xi=0)=1/3};
 

The only part worth concentrating on is the part I've made bold. I then can solve the system near xi=0 witha  series solution: