I've got a system of ODEs I need to solve for dependant variables phi(t) and a(t) with independant variable t.(appologies, it's a bit messy; the maple tag isn't making things render nicely)

dsys := {(D(D(phi)))(t)+3*a(t)*(D(phi))(t)/(D(a))(t)+dOm(phi(t), a(t)) = 0, a(1000) = sqrt(2*sqrt((8*Pi^3*(1/45))*To)*1000), phi(1000) = 12*M^(alpha+4)*alpha*g^2*sqrt(8*Pi^3*(1/45))*To^((1/2)*alpha-1)*10^(1/(alpha+2)), (D(a))(t) = sqrt(RhoDE(phi(t), a(t))+RhoM(a(t))+RhoR(a(t)))*a(t...