# I have a simple, two object thermodynamic model with radiation and advection. This model consists of

I have a simple, two object thermodynamic model with radiation and advection. This model consists of two first order quadratic differential equations, what I would like to solve analytically. The equations can be simplified to
$\frac{d}{dt}x={a}_{0}{x}^{4}+{a}_{1}{y}^{4}+{a}_{2}x+{a}_{3}y+{a}_{4}$
$\frac{d}{dt}y={a}_{5}{x}^{4}+{a}_{6}{y}^{4}+{a}_{7}x+{a}_{8}y+{a}_{9}$
I'm trying to find the analytic solution because I would like to use the model in an embedded system with limited resources, where solving it numerically is not feasible. So far I have tried to solve this with Maxima without any success. I have also consulted with my college textbooks, but they cover only linear differential equation systems.
You can still ask an expert for help

• Questions are typically answered in as fast as 30 minutes

Solve your problem for the price of one coffee

• Math expert for every subject
• Pay only if we can solve it

Jordin Church
I am afraid you are out of luck. Even finding the equilibrium is a very complicated endeavor. With general coefficients what you have for the equilibrium is a fourth order polynomial.
My suggestion is you try numerically, there is no general simple analytic solution. You don't actually need all that much resources to solve the problem numerically, since you can always try Euler:
Take $h>0$ a time step. The approximations for ${x}_{i}=x\left(ih\right),{y}_{i}=y\left(ih\right)$ are
${x}_{i}=x\left(ih\right),{y}_{i}=y\left(ih\right)$
The computational cost should be minimal and the approximation for small h is fairly accurate. If you find this method insufficiently accurate you can always try better numerical methods.