Hello,

I have a system of differential equations that can be used to model orbits of stars in binary systems in cartesian coordinates. Normally, I specify the masses of the stars, the period and eccentricity of the orbit to use for initial conditions, I use NDSolveValue and things work great. However, I want to solve the functions symbolically where I do not give initial conditions for the masses of the two stars and, instead get a solution that is mass dependent. I switched to DSolveValue, keeping all of the diff eqs the same and just didn't specify what the masses were. However, I keep getting an error message saying that things are not the same shape. Is it possible to solve these diff eqs without specifying the mass of the two stars?