Maple 2019.1 with "Physics" package version 419:
sol := dsolve([diff(x0(t), t) = 3*x1(t), diff(x1(t), t) = 0, diff(x2(t), t) = (2*x1(t)*x2(t) + x3(t))/x0(t), diff(x3(t), t) = (3*x1(t)*x3(t) - 6*x2(t)^2)/x0(t)])
I have:
[{x1(t) = _C2}, {x0(t) = Int(3*x1(t), t) + _C1}, {diff(x2(t), t, t) = (-6*x1(t)^2*x2(t) + 2*x1(t)*diff(x2(t), t)*x0(t) - 6*x2(t)^2)/x0(t)^2}, {x3(t) = -2*x1(t)*x2(t) + diff(x2(t), t)*x0(t)}]
If I add build
command to dsolve
I have a solution:
{x0(t) = 3*_C2*t + _C1, x1(t) = _C2, x2(t) = -9*(_C2*t + _C1/3)^2*WeierstrassP(((3*_C2*t + _C1)^(1/3) + _C3*_C2)/_C2, 0, _C4)/(3*_C2*t + _C1)^(4/3), x3(t) = -(3*_C2*t + _C1)*WeierstrassPPrime(((3*_C2*t + _C1)^(1/3) + _C3*_C2)/_C2, 0, _C4)}
In LaTex code:
$\left\{ {\it x0} \left( t \right) =3\,{\it c2}\,t+{\it c1},{\it x1}
\left( t \right) ={\it c2},{\it x2} \left( t \right) =-9\,{\frac {
\left( {\it c2}\,t+{\it c1}/3 \right) ^{2}}{ \left( 3\,{\it c2}\,t+{
\it c1} \right) ^{4/3}}{\it WeierstrassP} \left( {\frac {\sqrt [3]{3\,
{\it c2}\,t+{\it c1}}+{\it c3}\,{\it c2}}{{\it c2}}},0,{\it c4}
\right) },{\it x3} \left( t \right) =-{\it WeierstrassPPrime} \left(
{\frac {\sqrt [3]{3\,{\it c2}\,t+{\it c1}}+{\it c3}\,{\it c2}}{{\it c2
}}},0,{\it c4} \right) \left( 3\,{\it c2}\,t+{\it c1} \right)
\right\}
$
Check by odetest
:
odetest(sol, [diff(x0(t), t) = 3*x1(t), diff(x1(t), t) = 0, diff(x2(t), t) = (2*x1(t)*x2(t) + x3(t))/x0(t), diff(x3(t), t) = (3*x1(t)*x3(t) - 6*x2(t)^2)/x0(t)])
[0, 0, 0, 0]