I have been trying this system of differential equations in Mathematica but it doesnt seem to work. Any help would be much appreciated!
ParametricNDSolve[{D[a
